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The Council regrets that, as a consequence of the stoppage of work following 
the recent printing dispute, publication of the present issue of Monthly Notices 
has been so long delayed. 

It is inevitable that the appearance of the large number of papers now awaiting 
publication will continue to be retarded, but it is hoped in forthcoming issues 


progressively to reduce the delay and to restore the normal rate of publication by 
the end of the current volume. 
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NUTATION: COMPARISON OF THEORY AND OBSERVATIONS 
Harold Feffreys 
(Received 1958 December 11) 
Summary 


‘heoretical values are given for the principal, semiannual, and fortnightly 
nutations, based on the central particle and Roche Earth models. ‘These are 
compared with Federov’s observed values. Agreement is better than in any 
previous comparison. ‘The fortnightly terms are almost the same for both 
models and agree with observation. Both components of the semiannual 
motion agree best with the central particle model. ‘The observed 19-yearly 
motion in obliquity agrees best with the central particle model, that in 
longitude best with the Roche model, and some discrepancy would remain 


for anv mode i. 


1. ‘The rate of precession and the amplitude of the 19-yearly nutation, on 
the hypothesis that the Earth is a rigid body, are given theoretically by a pair of 
equations containing the Earth’s dynamical ellipticity (C— A)/C and the Moon’s 
mass. ‘Their observed values are the basis of the estimates of the two latter 
quantities given in the older lists of astronomical constants. ‘The position 
was changed by Hinks’s determination of the lunar inequality from the Eros 
observations of 1go1. ‘This inequality is a monthly motion of the Earth about 
the centre of mass of the Earth and Moon together. It was found to give the 
mass of the Moon with much smaller uncertainty than the previous method, 
and more recent procedure has been to use the lunar inequality to estimate the 
mass of the Moon and then the rate of precession to estimate the dynamical 
ellipticity. However Jackson (1930) pointed out that the results of the two 
methods were actually inconsistent. His conclusion was emphasized by Spencer 
Jones (1941) who estimated from the lunar inequality and the precession that 
the amplitude of the nutation in obliquity should be 9":2272 +0”-0008, while 


/ 
the observed values according to several reductions lay between 9’:207 and 


g”:213. ‘The standard value of Newcomb is g”:21 
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\ccording to the work of several writers, especially Poincaré (1910), allowance 


for elasticity of the Earth has no appreciable effect on the theoretical value. But 
sondi and Lyttleton called my attention to a passage in Lamb’s Hydrodynamics, 
1 on work of Poincaré, concerning the effects of a liquid core. No data 

il core had of course been available to Poincaré and to previous writers 

ect, including Kelvin, Greenhill and Hough, but the qualitative 

that a liquid core would shorten the free (Eulerian) period and 

of the 1g-yearly nutation. The latter effect depends 

ellipticity of the core boundary. Assuming that this is 

hydrostatic theory of the Earth’s interior, which is likely to be 

at such depths, I found (Jeffreys 1948) that with a rigid shell and a 

ous core the free period would be reduced from 305 days, as for a 
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Summary 


Theoretical values are given for the principal, semiannual, and fortnightly 
nutations, based on the central particle and Roche Earth models. These are 
compared with Federov’s observed values. Agreement is better than in any 
previous comparison. ‘The fortnightly terms are almost the same for both 
models and agree with observation. Both components of the semiannual 
motion agree best with the central particle model. The observed 19-yearly 
motion in obliquity agrees best with the central particle model, that in 
longitude best with the Roche model, and some discrepancy would remain 
for any model. 





1. The rate of precession and the amplitude of the 19-yearly nutation, on 
the hypothesis that the Earth is a rigid body, are given theoretically by a pair of 
equations containing the Earth’s dynamical ellipticity (C — A)/C and the Moon’s 
mass. ‘Their observed values are the basis of the estimates of the two latter 
quantities given in the older lists of astronomical constants. The position 


was changed by Hinks’s determination of the lunar inequality from the Eros 
observations of 1g01. ‘This inequality is a monthly motion of the Earth about 
th centre of mass of the Earth and Moon together. It was found to give the 
mass of the Moon with much smaller uncertainty than the previous method, 
and more recent procedure has been to use the lunar inequality to estimate the 
mass of the Moon and then the rate of precession to estimate the dynamical 
ellipticity. However Jackson (1930) pointed out that the results of the two 
methods were actually inconsistent. His conclusion was emphasized by Spencer 
Jones (1941) who estimated from the lunar inequality and the precession that 
the amplitude of the nutation in obliquity should be 9’:2272 + 0"-0008, while 
the observed values according to several reductions lay between 9”-207 and 
g":213. ‘The standard value of Newcomb is g’-210. 

According to the work of several writers, especially Poincaré (1910), allowance 
for elasticity of the Earth has no appreciable effect on the theoretical value. But 
Bondi and Lyttleton called my attention to a passage in Lamb’s Hydrodynamics, 
based on work of Poincaré, concerning the effects of a liquid core. No data 
on the actual core had of course been available to Poincaré and to previous writers 
on the subject, including Kelvin, Greenhill and Hough, but the qualitative 
results were that a liquid core would shorten the free (Eulerian) period and 
reduce the amplitude of the 19-yearly nutation. The latter effect depends 
considerably on the ellipticity of the core boundary. Assuming that this is 
given by the hydrostatic theory of the Earth’s interior, which is likely to be 
accurate at such depths, I found (Jeffreys 1948) that with a rigid shell and a 
homogeneous core the free period would be reduced from 305 days, as for a 
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rigid Earth, to 270 days, and the nutation amplitude by about 1 part in 160, 
about three times too much to account for the observed discrepancy, but in the 
right direction. 

Later work showed that allowance for elasticity of the Earth’s shell would 
again lengthen the free period and make the reduction of the nutation amplitude 
smaller, but a detailed quantitative comparison, based on an Earth model in 
accordance with modern geophysical information, was rendered possible by 
work of Takeuchi (1950) on the elastic yielding of the Earth under tidal forces. 

There is a complication due to the fact that the core is not uniform. The 
radius is about 0-55 of that of the Earth. Within it, however, there is an inner 
core, with radius about 0-36 of that of the main core, which may be solid. Its 
density is doubtful; the mean density of the core is about 10, but there may be 
a discontinuous increase of density at the inner core boundary by as much as Io. 
The effect of this is small for most observable phenomena since the volume 
concerned is small. 

Theoretical values of the corrections to the coefficients in the 19-yearly, 
semiannual and fortnightly nutations have been calculated by Jeffreys and 
Vicente (1958 a,b). ‘These are based on ‘Takeuchi’s solution for the shell, 
itself based on one of Bullen’s distributions of density, and on two models for 
the core. In one of these the core is replaced by a homogeneous incompressible 
fluid with an extra particle at the centre. In the other, which we call the Roche 
model, it is taken as having a quadratic density distribution, the variation of 
density being taken as wholly due to compression. In both the constants are 
adjusted to make the mass and moment of inertia of the core agree with a com- 
promise between Bullen’s extremes. In fact about half the departure of the 
core from uniformity of density (judged by the departure of the ratio J/Ma? 
from 0-4) is due to compression and half to the probable rise of density in the 
inner core, and interpolation between the two extremes should be at least a 
useful guide. 

In all discussions of observations until recently the ratio of the amplitudes 
of the nutations in obliquity and longitude has been assumed to be the same as 
for a rigid Earth. ‘This is not true for an Earth with a liquid core. The data 
have however recently been rediscussed by Federov (1958 a,b) who has used 

35000 observations for the 19-yearly term and 230000 for the fortnightly 
term; he quotes also a solution by Popov for the semiannual term based on 
observations at Poltava. ‘The coefficients in obliquity and longitude are estimated 
separately. ‘Thus there is now a basis for a detailed comparison between theory 
and observation. 


> 


Before proceeding to such a comparison I offer some remarks on the theory 


for a rigid Earth. The published theory is rendered more difficult by the facts 
that (1) the Euler angle ys (longitude of the pole) is sometimes measured in the 
direct and sometimes in the retrograde sense, (2) some theories refer only to the 
axis of figure, while others use an intermediary, namely the instantaneous axis 
(Oppolzer and Woolard (1953)), or the axis of angular momentum (Federov), 
and give displacements of the intermediary axis as the ‘‘nutation’’, while the 
difference between the axis of figure and the intermediary axis is called the 
‘diurnal nutation’’. I think that this separation makes for confusion. Either 
intermediary axis involves derivatives of the angular coordinates, and it is impos- 
sible to measure a time-derivative, since the error does not tend to o when the 


’ 
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interval between observations does. It is said that the instantaneous axis is 
determinable by observing through a sidereal day, but some of the nutations 
vary within a day by quantities larger than the error tolerated in the calculations. 
What is observed is the displacement of the vertical at a place relative to the 
stars, and for a rigid Earth the whole effect is determined by the displacement 
in space of the axis of figure, and there is no need for any intermediary. 

If 6, %, fare Euler’s angles, so that @ is the celestial colatitude of the pole of 
figure and y¢ its longitude, measured from west to east, we have as usual 


é + ¢cos 0 =w=constant (1) 
and the equations of motion, to the first order, for a periodic disturbance, are 
Ab + Cw sin 0 = «cos ot (2) 
Cw6 — A sin 0% = Bsinot. (3) 
For the fortnightly and semiannual nutations « is nearly —Scose, where ¢ is 
the mean obliquity. Put 
6=e+6', W=Acosot, J=%,+pt+y’, wv’ =psinot. (4) 
Then 
pas CwB — Aca _ Cwa—AcfB (5) 
Ctwtc— A2os’? C88 — A2a? 
In the elementary treatment (for instance in Tisserand’s book) 8 and # in 
(2) (3) are neglected, giving the approximations 
(A,n)=(—B/Cwo, «/Cwo). (6) 
But for the fortnightly term this introduces errors of about 1 part in 14. 
The displacements of the instantaneous axis relative to the axis of figure are 


a ae ro. 
(—sinOb/w, O/w)= (- ycosot, — “sin ot). (7) 
Ww Ww 

Thus if A’, »’ refer to the instantaneous axis 

(C*w?a — A*o)\’ = — CwB —(C—A)ou + ABo*/w (8) 

(C?w*a — A®o*)u’ = Cwa—(C— A)oB — Aao*/w. (9) 
On account of the smallness of C— A the terms in o on the right are small, and 
even the terms in o? approximately cancel the term in o? in the first factor on 
the left. Thus the neglect of second derivatives in the equations of motion 
gives the motion of the instantaneous axis with much greater accuracy than that 


of the axis of figure. This is pointed out by Woolard. 
The displacements of the axis of angular momentum relative to that of figure 


are APO nog ot, — Ado sin ot, and we find 
Cw c 


(Casto — Ato8)( ~ 7) = —CwB+ ol (10) 
Cw Cw 
so that the elementary theory actually gives the motion of the axis of angular 
momentum exactly, so far as a is concerned. This result is given by Federov. 
Woolard’s method of calculating A—A’ and 4» —p’ is to rewrite by transferring 
second derivatives to the right side, substitute the first approximations to @ 
and ys on the right, and integrate. This has the effect, in comparison with (5), 
of removing the term A’o* in the denominator. The resulting errors of order 
o?/w* of the motion are usually negligible, but not quite so for the fortnightly 


F sie 
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motion when the calculations are done to 0”-o0o01, but the difference will be 
neglected here. 

If o is positive, A, » are always greater than the first approximations and have 
the same sign, so that the motion is direct. If negative, these statements are 
reversed. As Federov points out, the axis of angular momentum has one definite 
advantage over the instantaneous axis, since it is totally unaffected by the 
14-monthly and annual variations of latitude. 

Since the dynamical theory as developed by Vicente and me takes a circular 
motion of the pole as standard («= — 8) and actually « is not quite equal to — 8, 
it is necessary to resolve each nutation into a pair of circular movements in 
opposite directions. 

3. Woolard’s values for the axis of figure, with his two parts stated separately, 
are as follows. 

Obliquity Longitude 
(cosines) (sines) 
Principal (§2) “+2100 —0” 0010 17"*2327 —0”*0034 
Fortnightly (2)) "0884 +0”-0059 0”*2037 +0”-0162 
Semiannual (20) "+5522 +0"-0028 1°*2729 + 0"*0075 
These need an adjustment since the latest work on the solar parallax and the 
mass of the Moon indicates that 9”:210 should be increased to 9”-2272 or 9”:2242 
according as Spencer Jones’s or Rabe’s parallax is adopted (Jeffreys and Vicente 
1957 b, p. 172). I take the former; it requires lunar terms to be multiplied by 
1-00187 and solar ones by 0-99598. I also multiply the coefficients in longitude 
by sine =0-39738 to give displacements. ‘Then altogether the relevant terms 
for a rigid body are: 
Obliquity 1,=9":2262 cos 9, +.0”-0945 cos 2) +0"*5528 cos2O. 
Longitude m,=6"-8594sin {. +0”:0875 sin 2) +.0":5068 sin2©. 
We have to resolve these into pairs of circular motions; we have 
1, +-im, =6' +1 s1n eds’ 
=8"-0428 e'? +17-1834e7% 
+0”-0910 e*) +.0”-0035 e~ 7! 
+0"+5298 e#© +0”-0230e-70, 
This corresponds to , of previous papers (Jeffreys and Vicente 1957 a, b). 

‘To derive ¢, the value with allowance for elasticity of the shell and fluidity 
of the core, we must multiply by the factors C/¢, given previously, the symbol n 
of the previous papers being o of the present. ‘Then for the central particle 
model (factors in Jeffreys and Vicente, p. 160) 

l+im=8"-0138e'2 + 1°-1877e 

+ 0”-0934 e?) + 0”-0038 e~? 
0”+5483 e” o casie™ 
g”:2015 cos §) + 0°:0972 cos 2 


i(6”"-826sin 2 +0":08g6 sin 2 


For the Roche model (factors in Jettreys and Vicente, p. 169) 


l+im 8”-0339 e' 1”-1848 e 
+ 0”-0934 e2!) f ‘ 
0”-5143 e7 
9":2187 cos , + 00971 cos 2) +0"*5 403 cos 2O© 
+ 1(6"-8491 sin Q + 0"-0897 sin 2) +0"-4883 sin2C 


’ 
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We have suggested (p. 169) that for the actual core the and 2) terms might 
be interpolated, with weight 0-4 for the central particle model and 0-6 for the 
Roche model. Such interpolation would not be justified for the 2© terms 
because there are resonances in the neighbourhoods, but it may as well be done 
as an experiment. Then the interpolated values would be 


1+im=g"-2118 cos Q, +0":0971 Cos 2) +.0"°5§35 cos2© 
+1(6":8399 sin 2. + 0":0897 sin 2) + 0"-5023 sin2Q©). 
Federov’s values (1958 b) are* 


0’ = (9-198 + 0-002) cos 2 + (0"-0949 + 0010) Cos 2) 
+ (0”-578 + 0”-004) cos2© 
sin € ys’ = (6"-853 + 0"-002) sin . + (0”-0918 + 0”-0010) sin 2) 
+ (0"-0533 +0”-0004) sin2©. 


There is some check on the uncertainties because Federov also estimates terms 
giving phase shifts, 
6’ =(—0"-001 + 0”-002)sin 2, + (0”-0019 + 0” 0010) sin 2) 
sin 4 ys’ = (— 0”-008 + 0-002) cos 9, + (0”-0019 + 0”:0010) cos 2). 


For reasons already given (Jeffreys 1950, p. 466; Jeffreys and Vicente 1957 b, 
p. 171), it is extremely unlikely that there are any measurable phase shifts. ‘The 
equilibrium amplitudes of the larger diurnal tides K and O are of the order of 
1ocm, and even if the phase shift was 47 they could only produce couples of 
the order of 10-5 of the nutational couples; and the evidence of tidal friction 
shows that the phase shift cannot be more than a fraction of a degree. Then it 
is reasonable to use the estimates of the damping terms to give an additional 
estimate of uncertainty. With weight 1 for the {) terms and 2 for the 2) terms, 
on account of the different numbers of observations, the standard error would 
be 0”:0037 for the {) terms and 0”:0027 for the 2) terms. These are on only 
four degrees of freedom but at any rate suggest that Federov’s uncertainties may 
be a little too low on account of correlation of errors. 

The cos 9 term in the obliquity is in agreement with the central particle 
model, but not with the Roche model or the interpolated model. On the other 
hand the sin {) term in the longitude agrees better with the Roche model. 

With the increased estimate of uncertainty the sin 2) terms agree with 
either model. 

The observations are affected by the tidal deflexion of the vertical. The 
value for a rigid Earth is multiplied by the factor 1+k—J/ (where / will not be 
confused with the direction cosine). Its theoretical values corresponding to 
components considered here are as follows: 


C.P. model Roche model 
— 2) *162 1°188 
—20 134 +226 
R 119 “161 
—-2 ‘121 163 
20 "149 “180 
2) ‘160 -183 





* A mistake in copying has affected the 2) terms. I give here revised values communicated 
to me by Federov in a letter. 
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Federov has applied a correction using the factor 1-1, which is a little too small, 
but as the whole effect is only in the third decimal (1958 a, p. 86) the difference 
is unlikely to be important. 

In any case the agreement between theory and observation has been much 
improved, though it is still not quite satisfactory. It should be noticed that the 
different behaviour of the obliquity and longitude terms is quite fundamental. 
It depends on the fact that the steady state solution is exact for speed o, for 
any model; and for speeds +o, with o very small, the correcting factors must 
be of the form 1+p. No change of model could account for the fact that each 
component of the §) terms agrees with one model and not with the other. It is 
in any case curious that three comparisons agree better with the central particle 
model, and one better with the Roche model, since the former makes no allowance 
for the compressibility of the Earth’s core, which is very well known, whereas 
the density of the inner core is still in much doubt. Compressibility alone 
would produce an increase of density towards the centre comparable with that 
assumed in the Roche model, and its effects cannot be negligible. 

Most correlations between errors would be likely to affect both the components, 
and also both phases for each component. It seems likely therefore that the use 
of the anomalies of phase to estimate an uncertainty is valid. In addition there 
is an unexplained systematic variation during the night, which was discussed 
at the Moscow I.A.U. meeting (I was afterwards reminded of it by Dr Atkinson). 
This might conceivably contribute to the semiannual terms but hardly to the 
19 yearly ones. 

It may be recalled that both models, with a correction for the effect of the 
ocean, give almost the same value for the free period, and this agrees with observa- 
tion. ‘The main object of the work has been to test whether Earth models 
consistent with the geophysical data are also consistent with the astronomical 
ones. Agreement would give a comprehensive check on the Earth models. 
The actual agreement is not quite satisfactory, but much better than in any 
previous comparison; and the discrepancies that remain are of a sort that could 
not, apparently, be removed by any change of the Earth model at all. 


160 Huntingdon Road, 
Cambridge: 
1958 December 10. 
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RADIATIVE RECOMBINATION OF HYDROGENIC IONS 
M. #. Seaton 
(Received 1959 February 3) 
Summary 
Using the first three terms in the asymptotic expansion of the Kramers- 


Gaunt factor, calculations are made for the rate of recombination and for the 
mean kinetic energy of the recombining electrons. 





For the study of the recombination spectra of gaseous nebulae one requires 
the rate coefficients for radiative recombination to individual quantum states of 
hydrogenic ions and for the study of the ionization equilibria one requires the 
radiative recombination rates summed over all quantum states (1). Further, 
for the study of the thermal balance in nebulae, it is necessary to calculate the 
rate of loss of electron kinetic energy due to recombination processes (2, 3). 
For hydrogenic ions all these quantities may be calculated exactly but the 
calculations become increasingly difficult with increasing m, where n is the 
principal quantum number. If errors as great as 20 per cent are of no importance 
one may put the Kramers—Gaunt g-factor equal to unity and hence obtain 
relatively simple expressions. A considerable improvement is obtained on using 
the asymptotic expansion of the g-factor, as derived by Menzel and Pekeris (4) 
and corrected by Burgess (§). From the expansion of g we obtain asymptotic 
expansions which enable the rate coefficients to be calculated with errors which 
will not exceed 2 per cent for temperatures of order 1o*°K or less but which 
may be greater for temperatures of order 10°°K. We present a systematic 
tabulation of the various functions which occur in these asymptotic expansions. 

1. Coefficients for recombination to level n.—Let «,,(Z, T) be the rate coefficient 
for radiative recombination to level of a hydrogenic ion, Z being the nuclear 
charge and T the kinetic temperature; the number of recombinations on n per 
unit volume per unit time is then «,N,N,. The recombination coefficient is 
given by 


>\12 
a,(Z, T)= 5 () (mkT)-32 2nte!n kt 


[ ” (hvPe-*"kPq (Z, hv)d(hv). (1) 


. In 
where 


hRceZ? 
I,= ) (2) 


is the threshold ionization energy for level m and R is the Rydberg wave number. 
The cross-section for photo-ionization from level n is 
: 2°aa,” n 
a, (Z, hv) = ——= =F, (1 + n*e)~*g1;(n, €) (3) 
v3 2 
3V 3 


where « is the fine-structure constant, a, the Bohr radius and 


hows hRet( =, + ‘) 
n* 
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Note that the energy of the ejected electron is Z%e in Rydberg units (13-60eV). 
The Kramers—Gaunt factor g,; in (3) is of order unity. 

From equations (1) to (4) we obtain 


. Xr} 2 
2,(Z, T) = FZ—x,S,(2) 


where 


28 (/n\12 ‘ 
Q= — (=) atca,? 
3 \3 


= 5197 x 10~'* cm® sec™ 


hReZ? 
kT 
72 
= 1578905, (T in °K), 


1 
’ 


A= 


x, =A/n? 


n, € e~tnu 
(A)= ” £11(%, €)e7** du 
— (1+u) 
with u=n’e. 
The asymptotic expansion for gy; is (4, 5) 
(n,€)=1+0°1728n-*3(u + 1)-*3(u—1) 
—0-0496n-*3(u + 1)-*3(u? + $u+1)+.... 
Substitution of (10) in (g) gives 
(A) = S(x,) +AM3,SM(x,) +A-23S@(x,) +... 


where 
S(x)=e7Ei(x), 


S(x) = 0°1728x!3X, (x), 
S® (x)= — 0-:0496x?°X,(x), 


6i(x) = ( ms 


+ 1)~53(4 —1)e-™ du, 


+ 1)~*3(u? + 4u+ 1)e—™ du, 


Extensive tables of the exponential integral, &i(x), are available (6). 
all finite x 


SO(x)=e74) —Inx—y + = 1 eee (18) 


where y=0°5772... is the Euler—-Mascheroni constant. In the limit of large x, 
S® has asymptotic form 


S(x) ‘ia 3, 1! ! 5 Se. cooPe (19) 


2 | oe 


z-< 
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The integrals X, and X, may be expressed in terms of confluent hypergeometric 
functions (7). We obtain 
X, (x)= —3+(1+3x)'¥ (1,4; x) (20) 


and 
X,(x)= —3(1+x%)+(1+2x4 $x?)V(1, 3; x) (21) 
where 
(c—1)¥(1,¢; x)= (c)x!~e? — O(1,€; x) 
and where the series 


2 

x x 
M(1,c;x)=1+- + —— +... 
c e(c+1) 


converges for all finite x. In the limit of large x, 





0 2 
seit x x 


' (3,¢: 2} « fy ad : enys— 9 —..}. 


(24) 
The functions X, and X, have been evaluated using (23) for «<2, numerical 
integration for 2<x< 10 and the asymptotic expansion (24) for x>10. Table I 
gives xS®, xS® and xS® for o-o2<x<20. The intervals in x are such that 
linear interpolation will introduce errors no greater than 2 in the last significant 
figure. ‘lo the accuracy of Table I, for small x, 


vSO =x e7{ —Inx—0°5772+x!} 


v¥ SM = 0°4629x(1 + 4x) -0368x49(1 1 


> 


“© 0°0672x(1 + 3x) + 0°14885°9(1 + 


and for large x, 


x S“ 
S 325 — n - ae 11440] x2>20. (26) 
x) (3x)? (3x)® = (3x)*) | 
rS2 0:0496x23 I 3 \ 32 448 ) | 
J 


aa = 2 \3 
L (3%) (3x)? (3) 
2. Total recombination coefficients.—The total rate coefficient for radiative 
recombination is 
@ 

Ztot i) wf 
n° =] 
For large n’ the sum may be replaced by an integral. 


ny” 
’ 


From 
(f° oO 


a,,-dn’ 


we obtain 


. >; = : 

a, DZ"? ¢ ¥n S,(A)+¢,,(a ) 
oes ln if 
where 


and, using (11), 
a, (A) =o (x,,) AM (x, ) +A*B (x) +... 
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where 
o?)(x) = [ S (x) dx’. 
Using (18) it may be shown that ‘ 
o (x)= SO(x)+Inx+y. 


The functions o, o® and o® are given in Table II for 0-02 < 20. 


a =xd(1 4 *)(—Inx— me 4 


For small x, 


\ 


‘1) i] 
ol Ps 


J 


of) = “4 — 0°0672 (: + :) + 0:0893x* (: + x) 


and for large x, 


f - 
"Pa 

- x4 0°4629(1 + 2x)—0°7776x13 (14+ =. 
a”) 


£0) 


o = SO +1nx+0°5772, 


7 
a) = 0:926 —0°5184x'3 <1 4+ on 
L (3) 
2 my | 6 
a = 0-134 —0°0744x73 414+ — 
\ 
- =. - 
lable IIL gives values of x,S,, x,S,/2 and }x,,S,,/n for various values of (T/Z?). 
1 
As compared with the results obtained using the asymptotic expansion, the 
expression 


ttot = ZZA'*{0-4288 + } InA+0-469A-} 3} (36) 


is in error by less than o-5 per cent for (7T/Z*) < 10°, by 3 per cent for (7T/Z*) = 10° 
and by 31 per cent for (7/Z?)=5 x 10° °K. 


Taste III 
(T/Z*) x,S, Xoo ad Cae 
K 2 ~~ + 
10° ‘775 0°427 3°046 
5 x 10° 0°763 0°399 2°305 
10 0'748 0°371 1°998 
gx r0* 0660 0260 1°324 
10° 0'585 0°202 1°055 
5 x 10° 0350 0089 0°529 
10° 0'252 0:056 0°364 
5 x 10° 0094 0018 0°127 


3. The energy of electrons recombining on level n.—Let kT8,N,N, be the 
kinetic energy, per unit volume per unit time, of the electrons recombining on 
level n. We obtain 


AuNz 
B,(Z, T)= GZ — x,"R,,(A) (37) 
where 


R,(A)= [ ceseny udu. (38) 


Using (10) we have 


R,,(A) = R(x,,) +A*3RO(x, ) +A-23 R(x, ) 
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xR® 


044 
043 
"043 
"042 
"042 
"042 
‘O41 
‘O41 
‘O41 


"040 
"040 
*039 
039 
°039 


038 
*037 
036 
035 


"034 
"034 
"033 
"032 
"031 
"031 


where 


p™ 


‘0029 
"0043 
‘0056 
‘0069 
‘0081 
"0093 
*O104 
‘O1IS 
0126 


"0146 
"0165 
0184 
‘0201 
‘0217 


"0255 
"0288 
0317 
0343 


0377 
‘O41 
"0445 
°0477 
"0509 
"0540 


TaBLe IV 


p 
— -0009 
‘0013 
‘0018 
‘0022 
*0026 
‘0031 
"0035 
“0039 
"0043 


‘OO51 
"0059 
‘0067 
"0075 
0083 


‘OI02 
“OI2I 
*0O139 
*O157 


‘OIQI 
"0225 
"0259 
‘0291 
"0323 
‘0354 


xR = 1—xS, 


ate 
xRO = — o-1728x8 — X,(x), 


xR® = 0-0496x° cs X,(x). 


(40) 
(41) 


(42) 


The functions R® and R® are given in Table IV for 0-02 <x<20. For small x, 


xR® = 0°1543(1 — 8x) + 1°94424 


¢R® = —oo448(1 - 3 x) —0°268x° 


and for large x, 


xR® = —0°1728x | 


xR® = —0-0496x us y 


r 


( 


re 16 an 
(3«) (3x)? 


(3x) | (3x)? 


* £002, (43) 


3200 | 


Gx)” 


1792 
(3x)* 


57200 
(3«)* } 


x2> 20. 


For the total kinetic energy loss due to recombination we require 





ha 
Brot = > By 
n’=1 
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For large n, 


2 ) 
dn’ =~ GZX +f = R,,(A) + pa(A) (46) 


where 
pa (A) = p(x) + A291 (x,,) + A-29(x,) 40. (47) 
and 
p(x) =| x" R(x’) dx’, (48) 
0 
We obtain 
p(x) = (1 — x) SO (x) + Inv + y = 0 (x) — xS(x). (49) 
The functions p™ and p® are given in Table 1V for o-o2<x<20. For small x, 
p™ = 0°1543x(1 — 4x) +0°833x78 ] 
P (50) 
p~ 


and for large x, 


(3x)? 
48 


— ae + > 


(3x) (3x)?  (3x)¥J 


x 
Table V gives values of x,?R, and > x,?R,,/n for various values of (T/Z?). As 
l 
compared with the results obtained using the asymptotic expansion, the 
expression 


Bro. = ZZA'24 — 0-0713 + —InA+0-640A-13} (52) 
is in error by less than 1 per cent for (7/Z*)<5 x 10* and by 4 per cent for 
(7/Z*) = 10° °K. 

TABLE V 
(T/Z?) © x,2R, 
K 4 anid n 
10° wf i 2°578 
5 x 10° m 1°857 
10! , 1°564 

5 x 108 , 0926 

10° . 0°679 


4. Accuracy of the calculations.—The asymptotic expansions will be accurate 
if the third term is much smaller than the sum of the first two terms. For a 
given temperature the accuracy will improve with increasing n and, for a given n, 
the accuracy will improve with decreasing temperature. 

In Table VI some values of x,,S,, calculated from the asymptotic expansion 
are compared with the results of more accurate calculations made by Burgess (5) ; 
the results of Burgess are exact in the limit of To. For (7/Z?)<2x 10*°K 
the maximum error in «, will be 2-5 per cent and the errors in «,, for m>1 will 
be smaller than 1 per cent. 
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TasBle VI 


% yS''(x,,)  A~tx, SO%%e,) A~4x,S'(x,) XS (A) x,S,(A) Difference 
(Asymptotic) (Burgess) 
‘1728 — 0496 7776 7973 +0197 
‘1089 —*O197 8714 8763 +:°0048 
‘0831 —"O1Is "9054 "9075 +°0021 


"1480 "0470 748 763 +°OI§ 
0645 “0169 "742 “748 + -006 
"0305 “0089 661 *660 —‘Ool 


"1293 "0451 "723 726 +°003 
0436 "0156 661 663 +002 
0148 0081 548 “549 +°0o1 


From the magnitudes of the various terms in the expansions, it is estimated 
that our results for aot may be in error by about 10 per cent for (7/Z*) = 10° 
and by about 25 per cent for (7/Z*)=5 x 10° °K. 


University College, 
London, W.C.1: 
1959 February 2. 
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THE SOLUTION OF CAPTURE-CASCADE EQUATIONS 
FOR HYDROGEN 


M. F. Seaton 
(Received 1959 February 3) 


Summary 


Calculations are made for T'=2:'5 x 10°, 5 x 10°, 1x 104 and 2x 104 °K 
taking an infinite number of levels into account but assuming that 
b,,= ((2l+1)/n?] 6,. The results obtained by Baker and Menzel are 
shown to be in error by about 5 per cent for Case A and 20 per cent for Case B. 

An improved method of calculation results from the introduction of the 
cascade matrix, the elements C,,,, giving the total probability of the n’+n 
transition when all cascade routes are taken into account. 





1. Introduction.—In Case A of Baker and Menzel (1), valid for an optically 
thin nebula, it is assumed that the excited states of hydrogen are populated by 
radiative capture and by cascade from higher states and that they are depopulated 
by radiative transitions to all lower states. In Case B, valid for a nebula which 
is optically thick in the Lyman lines, it is assumed that the rate of depopulation 
of excited states by emission of Lyman lines is exactly equal to the rate of 
population by absorption of Lyman lines. In Case B, for all m>2, one may 
equate the rate of population by capture and cascade to the rate of radiative 
transitions to all lower states n” with n”>2. For both cases the equilibrium 
equations may be written 


N, N +e be | S N vA, Y.nal= N A 1 (1) 


ll ia wal nl 


n'l’ 
where n’ is greater than m. ‘The number of electrons per unit volume is N,, 
the number of protons N. and the number of hydrogen atoms in quantum 
states ml is N,,. ‘The radiative recombination coefficient for n/ is x,, and the 
radiative transition probability for n'l’->nlis A, ,,.. The total radiative transition 
probability for n/ 1s 


> Ant n'1 (2) 
"wht 


where* n,=1 for Case A and n,=2 for Case B. 


l 


In conditions of thermodynamic equilibrium all quantum states of the same 
energy have equal populations. We then have 


Now +) (3) 


n* 
where 


ni 


N= > Noy 
l 


* It is, of course, to be understood that Any, »*;" is identically zero for n” < (/”+1); hence it is 
only for /=1 that the 4, are different for Cases A and B. 
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Although the solutions of (1) will not satisfy (3) we consider the approximation 
of simplifying (1) by assuming (3). We then have, in place of (1), 


NN io ‘uit 2 NwAw, a N,An (5) 


n’=n+1 


a, = Loup (6) 


I ’ 
Ay. ha a = (2/ +s )A wr, nl (7) 
nil 


‘Aon i Ve wat for Case 3): (8) 


Pr Ny = 2 for Case B 


Baker and Menzel (1) have considered equations (5) with an infinite number 
of levels. Burgess (2) has obtained solutions of the equations (1) in the 
approximation of considering only 12 levels and has then applied a correction 
obtained on comparing solutions of (5) for 12 levels and for an infinite number 
of levels. Searle (3) has obtained solutions of (1) for 10 levels. 

Burgess has drawn attention to some errors in the work of Baker and Menzel. 
An error in the third term of the asymptotic expansion of the Kramers—Gaunt 
factor, as obtained by Menzel and Pekeris (4), gives rise to errors of a few per cent 
in the calculated hydrogen level populations. There are more serious errors in 
the Case B calculations of Baker and Menzel. Burgess notes that for T= 2 x 104 °K, 
some of the N,, given by Baker and Menzel in the approximation of considering 
only 14 levels are greater than the corresponding N,, given for an infinite number 
of levels. This cannot be correct since the inclusion of more levels should increase 
the rate of population due to cascade and hence increase all of the N,,. Burgess 
checked the solutions for 14 levels and was able to conclude that there were errors 
in the Baker and Menzel Case B calculations for an infinite number of levels. 

It may be shown that the solutions for hydrogenic ions of nuclear charge Z 
at temperature T may be obtained from the solutions for hydrogen at a temperature 
of (T/Z*). In planetary nebulae T is usually between 1 x 104 °K and 2 x 10* °K. 
To obtain results for H1 and Het! we therefore require calculations for hydrogen 
in the range 2-5 x 10°°K to 2x 10*°K. The lowest temperature considered 
by Baker and Menzel is 5x10*°K. The need for calculations at lower 
temperatures has been emphasized by Wurm (5). 

In the present paper we consider equations (5) for an infinite number of 
levels and obtain solutions for T=2-5 x 108, 5 x 10%, 1x 10* and 2x 104 °K. 
These calculations correct the errors in the work of Baker and Menzel and extend 
the solutions to lower temperatures. Our main interest in this work has been 
in the development of convenient methods for the solution of the infinite set 
of simultaneous linear equations for all values of n. Similar methods will later 
be used for the solution of exact equations of type (1), both for hydrogen and for 
complex atoms. 

2. Solution of the equilibrium equations.—Let P,,,, be the probability that 
population of n’ is followed by a direct radiative transition tom. Then 


A, n 
Py n= a: (9) 
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Let C,,. ,, be the probability that population of n’ is followed by a transition to n 
via all possible cascade routes. We then have 


— 
n=l 


a . 
+2. 0 Poe vil n rn tf Pavan ’ (10) 


n+1,n 


a Y Y A > 
14+3,.n= Pass n+al n+a,n bt Pnse, n+1Cn+.0 + Pass, n 
With the convention C,, , =1 we have, quite generally, 


n 
om = > | ee om n* (11) 


p=n'—1 
We shall refer to C,,.,, as the cascade matrix. Using (11) the cascade matrix 
may be generated once the matrix P,,, , is known. 
The solutions of equations (5) are given by 


N,A,=N.N4 > tC ns (12) 


n’ = 
The form of the solutions used by Baker and Menzel may be obtained by the 
following argument. We consider level m to be populated by: direct capture 
on m; captures on higher levels followed by direct radiative transitions to n; 
captures on higher levels followed by cascade transitions to n via one intermediate 
state; cascade transitions via two intermediate states, .... The solutions are 
then 


( ca oo n’—1 


rT T , ) * ) > 
NN 04% t+ > ayPyrnt 2 a ae A ee 
q n=n+1 n’=n+2 p=n+1 


(13) 
Re-arrangement of terms shows (12) and (13) to be identical. The particular 
advantage of (12) is that one may calculate the cascade matrix once and for all, 
this matrix being independent of the particular processes populating the excited 
states. It may then be used for radiative recombination at all temperatures and 
also for non-radiative processes such as collisional excitation from the ground 
state. The only assumption made* is that the excited states are depopulated by 
radiative transitions alone. ‘The cascade matrices are, of course, different for 
Cases A and B, due to the difference in the definition of A.,,. 
It is usual to solve for 5,, defined as the ratio of N,, to the corresponding 
population for thermodynamic equilibrium at the same density and temperature. 


We then have 
N \ \ nb “n oa a ( I 4) 
4 =i d One 
n ¢ n ; LT 


where x, =/,,/kT, I,, being the threshold ionization energy for level n. From 


(12) and (14), 
T\ 32 20 


3. Formulae for hydrogenic ions.—F¥ or hydrogenic ions with nuclear charge Z 
the radiative transition probabilities are given byt 


4 4 
A. = (=) - (16) 
37Agy/ 3j/ n°? 


* We here consider only the essential assumption made in our definition of the cascade matrix 
and disregard the incidental assumption implied by (3) ; one may easily define the cascade matrix 
when (3) is not assumed. 


+ Where possible we use a notation similar to that of Baker and Menzel. 


a 
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where « is the fine structure constant, a) the Bohr radius, c the speed of light and 
2n'g,,, 
ty rt de! = I 
n, n(n’? —n?)’ ( 7) 
Z,’., being the Kramers—Gaunt factor of order unity. Defining 


= ban” (18) 


n’ = ny 
Pyin= saan (19) 
The ionization threshold for level n is F 
Ls hRcZ* (20) 


nz 


we have 


where R is the Rydberg wave number. We put 
ere hRcZ? _ X (21) 
 —. # 
where hReZ* 15-789 x 10*Z? +: ed 
A = = rm . 
;T F (T in °K) (22) 
The radiative recombination coefficients are given by 


‘aii (*t (2) , aca) zysa*0Sa(%) (23) 


3 
” g11(n, €)e~7n™ 


S.0)= | 2 woke (24) 


and where ARcZ*e is the energy of a free electron, g,;; the bound-free 
Kramers—Gaunt factor and u=n*e. From (15), (16) and (23) we obtain 





where 


3 oc 
etnh = ll _* Swlw,n “4 2 
n y™ 3 
nn=n n 


In the notation of Baker and Menzel, 
(S,+V,,). 


e7nb,, = ; 


(26) 


therefore 


(27) 


4. Methods of calculation and extrapolation procedures.—Three-figure tables 
of the factors g,.,, are given by Baker and Menzel for n, n’<30. Accurate 
six-figure calculations have been made by Green, Rush and Chandler (6) for 
n,n’ <20. Spot checks show that most of the values of Baker and Menzel are 
correct to three figures. Errors which are larger, but not serious, occur for 
n,n’ large and (n’—n) small: thus goo ,,=0°782 according to Baker and Menzel 
compared with the exact value of 0-7782.... In all of our calculations the Baker 
and Menzel g,,, , were used. 

Table I gives calculated values of t, for Cases A and B and Tables IIa and 
IIB give the cascade matrices C,,,,, , for n<10, m<10. For 5<m<to it is 
found that the calculated C,,,,,, , can be fitted accurately to 


A 
as =C, ns 28 
nena Cont (28) 


Tables II include the limiting. values, C,,, ,. 
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The functions S,(A) may be calculated using the first three terms in the 
asymptotic expansion of the factor g;;._ Tables for the calculation of S,, are given 
in the preceding paper (7). For n’ large the sum in (25) may be replaced by an 


TABLE I 


tn 
A. 





~ 


A B 
I°gII — 
3°084 1°363 
3°929 2°265 
4°589 2°951 
5*132 3°507 
5°590 3°973 
5°989 4°378 
6°339 4°731 
6-657 5°O51 
6°943 5°340 
7°205 5*604 
7°442 5°841 
7°668 6-068 
7°876 6°277 


ont AM b&w DN 


8-740 7°144 
25 9°410 7°814 
30 9°956 8-360 


integral and in this integral C,,, ,, may be replaced by a mean value independent 
of n’. We then have 





TABLE 
Values of C 
(Case 

5 


a 


3 
oe 


"4418 
"4105 
"3991 
"3934 
"3903 
"3883 
*3870 
*3860 
"3853 
*3848 


Qu WSN 


~ 
oo Os! 


8 


*3796 
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Tas_e IIB 
Values of Casm, n 
(Case B) 


5 6 7 8 9 


+2890 "2454 2163 *1961 +1804 
2458 "2015 +1726 1531 *1382 
2282 -1836 "1550 “1351 *1206 
-2187 1737 "1451 "1253 *1107 
2128 -1676 *1389 ‘1190 *1045 
*2090 *1635 *1346 -1148 *1003 
+2062 *1605 "1315 “1117 *OQg7I 
"2041 +1584 "1292 “1094 0948 
"2025 *1567 "1274 *1076 "0930 
"2012 "1553 *1260 *1062 ‘0916 


of° on aunt wW N 


& 


‘2721 “1895 1431 *1130 "0934 0787 0680 


Tables for the calculation of o,(A) are given in the preceding paper. For n<10 
it was found that good accuracy was obtained with g=n+ Io. 

For n>10, (S,,+V,,) was obtained by extrapolation. This may be done 
quite accurately. Taking all Kramers-Gaunt factors equal to unity, Menzel 
and Baker (8) show that, in the limit of large n, t, ~3Inn, S, ~2lnm and 
V,,~Inn. Using the methods of Menzel and Baker it may be shown that the 
same results are obtained with the asymptotic expansions of the g-factors. It 
follows that 6, ~1 in the limit of large n. The values of t,, given in Table I may 
be fitted to ¢, =3Inn+ const. for n> 20. 

For 5<n<10 it was found that the calculated values of (S,,+V,,) could be 
fitted accurately to 


’ ee 
(S,+V,)=3lnn+e+ = + 5. (31) 


These expressions were used for extrapolation to larger values of n. 


5. Results—Table III gives values of e*»b, for the four temperatures 
considered and Table IV the Balmer decrements (relative emission intensities), 
for 7 =1 x 10 and 2 x 10%, calculated from 


I, = r00($) "pene (32) 
n} byeZ4 9 
Our relative intensities are similar to those of Baker and Menzel but there is 
some change in the (H,/H,) ratio for Case B. In Table V we compare our values 
of 6, with the values of b, obtained by Baker and Menzel and the values of b, 
obtained by Burgess, 5,, being the effective value of b,, for the calculation of Balmer 
intensities. The errors of order 5 per cent in Baker and Menzel’s 6,, for Case A 
may be explained in terms of the error in g,;. The larger errors, of order 
20 per cent, in their 5, for Case B may be explained by assuming that the values 
of (S,,+V,,) used in (26) were in error due to the error in g,; and in addition 
that the values of t,, used in (26) were calculated with g, ,,.=1. This explanation 


seems plausible in view of the fact that Baker and Menzel made some calculations 
with g,, ,”=1. 








M. F. Seaton 
TABLE III 


Values of e®nb,, 
B 


= 





m 
I x 10! 2x 104 


= vf 
x 104 axi1e agxi0" §x 10 


3 


‘193 "315 
‘213 "332 
244 "364 
‘273 "394 
"299 "421 
"322 "443 
‘341 “463 
*360 *480 
‘376 "494 
"390 ‘507 
"402 "519 
“414 "529 
"424 "538 
‘434 ‘547 


wu UU 
+ NK > 
ov 08 O 


wm 


>= NNW WN 


“om 
HOO CON Audbwrn ~~ 


wu wn 
ont mS 


~ 
N 


~~ + 
ub WwW 


‘472 "580 
"499 "603 
*520 ‘621 
TaBLe IV 
Balmer decrements (values of F ,,) 
A 


A. 





. 
1x 104 2x 104 


191 199 
100 100 
58-9 56°9 
37°8 35°6 
25'°8 23°8 
18°4 16°7 
13°7 12°2 
10° g'2 
8: “: 
° 


NNW? Us 
WS Of NP eH UN N 


oO ™ 
wu odo 


TABLE V 
A B 


A — % 
— —~ 


= 1x 104 2x10 
b,, Baker and Menzel . -166 "404 
d,, present results , ‘201 "484 
b,, Burgess ° "200 "494 














No. 2, 1959 The solution of capture-cascade equations for hydrogen 97 


For Case B the calculations of Burgess (2) and of Searle (3), based on 
equation (1), are in fairly good agreement with the results of the present 
calculations. For Case A the agreement is less satisfactory. This is due to the 
fact that, in calculations using (1), transitions to the ground state take place 
only from the np states; equations (1) therefore give differences between 
Cases A and B which are smaller than the corresponding differences obtained 
using (5). 


Department of Physics, 
University College, 
London, W.C.1: 
1959 February 2. 
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Summary 


There appears to be a practical limit to the resolution attainable in solar 
photography from a ground station. In an attempt to obtain greatly improved 
resolution, photographs have been made with a 29 cm refracting telescope 
mounted underneath the nacelle of a manned balloon at a height of 18 oo0 ft. 
Two flights were made; during the second flight on 1957 April 1, 480 frames 
were exposed. Although the ground seeing during the flight was very poor, 
the photographs taken with the balloon-borne telescope during this second 
flight are superior to the best that have been secured at any ground station. 

The complex structure of the photosphere is illustrated by photographs. 
Autocorrelation curves for two scans across the best photograph are given. 
The contrast transmission functions for the complete telescope (objective, 
eyepiece, photographic emulsion) have been measured using sinusoidal 
gratings, and the results used to correct microphotometric measurements of 
contrast. The corrected mean contrast between granules and intergranular 
regions is 40 per cent at a wave-length of 5300 A. 





1. Introduction The fine structure of the solar photosphere, known as 
granulation, has received a great deal of attention during the last hundred years, 
but progress towards an understanding of its nature has been slow. Observations 
are greatly hindered by the daytime seeing, which is usually of poor quality: 
while some of this poor seeing has its origin in the upper atmosphere, much 
of the trouble is due to local heating near the ground or around the observing 
slit of a dome. Also, unless the telescope is carefully designed, instrumental 
imperfections arise, chiefly from the disturbing effects of the large flux of solar 
radiation through the apparatus. ‘This is especially true of local heating near 
the prime focus. 

In the past, good photographs have been taken by paying particular attention 
to the design of the telescope and to the selection of an observing site. We 
should mention here the early pioneering work of Janssen (1) and of Hansky (2). 
The well-known photograph taken by Janssen on 1885 July 5 has scarcely 
ever been improved upon, but it is not often mentioned that it has been selected 
as the best of an extensive collection of photographs taken at Meudon over a 
period of 20 years with a telescope of only 13-5 cm aperture. Better photographs, 
taken by Lyot at the Pic du Midi, were investigated by Macris (3); Lyot’s 
techniques have been still further improved by Résch (4,5) using a refractor 
of 25cm aperture. A smaller telescope, of 13 cm aperture only, has recently 
been designed and used by Bray, Loughhead and Burgess (6,7) with good 
results. Other investigations of solar granulation will be referred to later. 
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Even under the best observing conditions at a good site the proportion 
of first-rate photographs has always been disappointingly small. In all of 
the work described above the theoretical resolution of a 25 cm objective 
has never been realized. As some of the bad seeing originates close to the 
ground we may reasonably expect that if photographs could be taken from 
a balloon at a height of a few thousand feet above the ground, better 
definition would result. The present experiments using a manned balloon 
were made to test this supposition and obtain photographs of still higher 
resolution. More recently, Schwarzschild and Rogerson have taken photographs 
from an unmanned balloon at heights of up to 80000 ft (8, 9). 

The expected scale of the fine structure is very close to the resolution limit 
of a telescope of about 25cm aperture. If quantitative deductions about 
granulation are to be made from the photographs, the least that we can do is 
to test the laboratory performance of the telescope and use this information to 
correct the observed data as well as possible. This also has been done in the 
present investigation. ‘The photometric measurement of granulation is closely 
analogous to the investigation of Fraunhofer line profiles with a spectrograph 
of rather low resolving power. It would be unwise to attempt to measure a 
line profile knowing almost nothing of the instrumental profile of the 
spectrograph; and almost as unwise to assume that the instrumental profile 
is that corresponding to simple diffraction, without aberration or scattering in 
the photographic emulsion. Similarly, it is unwise to deduce the photometric 
properties of granulation knowing that the resolution is only just sufficient to 
show the granulation, but not knowing the precise properties of the telescope. 

2. Apparatus.—Two flights have been made. Accounts of the apparatus 
used on the first flight from the Meudon Observatory on 1956 November 22 
have already been published and only a few notes will be given here. The 
telescope and balloon are described in ref. (10). The telescope was a refractor, 
with a visual doublet of 29 cm aperture and 294 cm focal length which had 
been figured by Schaer. The primary image was projected by an eyepiece of 
2°5 cm focal length to give an image of diameter 46cm in the second focal 
plane. This final image was photographed on Kodak Microfile Pan film through 
a Wratten No. 58 filter with a peak transmission of 5300 A and a bandwidth 
at half-peak transmission of 600 A. The shutter of a Contax camera was 
placed in the prime focus and set to a speed of 0°8x10-% sec. It was later 
discovered to have given an exposure of about 2 x 10~* sec because of the low 
temperature at altitude. The telescope was so supported on an altazimuth 
mounting underneath the balloon basket that it could be directed at the Sun 
by one of the two observers, whilst the second could control the mechanisms 
of the telescope (primary shutter, exposing shutter, camera wind, focusing, 
alteration of shutter speed) electrically from the basket. The automatic 
mechanism for operating the Contax is described in ref. (rr). When the balloon 
is floating at equilibrium altitude and the observers in the basket are still, the 
telescope does not oscillate appreciably. The principal effect is an occasional 
slow rotation about the vertical axis. Under reasonably calm conditions it 
proved possible to keep the telescope directed at the Sun with an error not 
exceeding 20’—30’ for a few minutes atatime. On this first flight 390 photographs 
were taken at heights between 18 000 ft and 20 000 ft, but the photographs 
of granulation had a sheared structure caused by the movement of the telescope 
during the unexpectedly long exposure. 
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On 1957 April 1 a second flight was made from Meudon (12) with similar 
apparatus, but including a much improved shutter. This publication shows 
a photograph of the balloon and apparatus. ‘The Contax shutter was replaced 
by a flying slit shutter of original design which gave a speed of about 2 x 10-4 sec. 
This was sufficiently fast to prevent the appearance of shear on practically every 
photograph even though the balloon was less stable than during the first flight. 
An outline drawing of the shutter mechanism is shown in Fig. 1. An essential 
requirement of the shutter was that it should be as nearly inertia-free as possible, 
to minimize vibrations transverse to the optical axis of the telescope just at the 
moment when the photograph is taken. The shutter proper consists of the 
framework A pivoted and balanced about a horizontal axis and containing in 
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Fic. 1.—Diagrammatic representation of components of flying slit shutter. For clarity the 
bearing that supports the nearer end of the horizontal axle of the shutter is not shown. The shutter 
is here loaded, ready to be released by raising the latch G. 


its upper part the adjustable slit, S, which makes the exposure. ‘The slit is 
set to a width of about 0-4 mm. In operation, we require this slit frame A 
to be accelerated to the required angular speed and then to be allowed to continue 
its motion at constant speed whilst the slit which it contains moves across that 
part of the primary image which is being photographed. At the end of its run 
it must be stopped positively without any rebound. These operations are 
performed by the mechanism shown. ‘The spring loaded arm C can rotate 
freely about the same axis as the slit frame A. The sequence is initiated by the 
solenoid B which lifts the latch G and so releases the slit frame. The arm C 
is then accelerated by the spring and carries the slit frame A with it, until the 
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slit is just before the area of prime focus to be photographed. Part C is then 
brought to rest by a cushioned stop, D, but the slit framework A continues 
in free motion at constant speed until it comes against the stop E. This stop 
is a small permanent magnet which, attracting a steel shoe on the slit frame A, 
prevents any rebound. 

The shutter is reloaded by a solenoid at F, which pushes the whole shutter 
framework until it is caught by the latch G. The shutter is enclosed in a light- 
tight box which contains, however, an entrance aperture which is normally 
covered by a preliminary shutter. This shutter is moved out of the way by an 
electromagnet just before each 0-2 millisecond exposure is made. 

In flight the speed of the shutter could be changed by altering the tension 
of the spring by means of the motor-driven wheel H. The speed was measured 
in the laboratory by illuminating the shutter with a neon lamp fed from a beat 
frequency oscillator. The light after passing through the shutter then shows 
a pattern of light and dark stripes, and the speed of the shutter can be calculated 
from a knowledge of the scale of this pattern and the width of the slit. There 
is a limit to the increase of speed that can be effected by diminishing the slit 
width. The reason for this is that during the flight the position of the shutter 
with respect to the primary focal plane is changed. The change in focus is 
very difficult to predict theoretically. It depends partly on the known change 
in refractive index of the air at altitude. But it depends also on thermal changes 
in the dimensions of the metal framework and the two components of the 
objective, and the temperature at which these attain thermal equilibrium with 
their surroundings is uncertain. It proved impossible to predict the position of 
focus with an accuracy better than +2mm. We must evidently be careful that 
the shutter slit is not displaced so much from the focal plane that the whole of 
the objective cannot be seen through it at once, and, as the focal ratio of the 
objective is f/10, the slit must not be less than o-2 mm in width. 

For the second flight a smaller magnification was used, the image diameter 
being 14:4.cm and the scale 6"°-7=1 mm. During this flight 480 frames were 
exposed at an altitude of 18000 ft. For the reasons just explained the 
position of focus is uncertain and the projection eyepiece was moved through 
©-I mm every ten photographs over a total range of +2mm. During both | 
flights a photoelectric photometer was used to estimate the exposure. The 
film was processed in strips and developed in Kodak D.19b for 5 minutes at 
a temperature of 70 °F. Photographs of a standardized step wedge, previously 
made on portions of the same batch of film, were developed with each strip 
to provide a calibration. The emulsion gave a contrast of y=6-1. 


3. Measurement of contrast transmission functions of complete apparatus 


3.1. Introduction.—The photography of granulation is extremely difficult 
but the deduction of the photometric properties of granulation from calibrated 
photographs is even more difficult. Not only are we observing an extended 
object of low contrast, but we are also working near the resolution limit of a 
medium-sized objective. In visual tests with an eyepiece on a point object the 
objective may appear to perform well, showing clear and regular diffraction 
rings, but a test of this kind does not easily show how much spurious light 
there is in the outer parts of the image; such spurious light diminishes the 
contrast in the image of a complex object if the scale of its structure is of the 
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order of magnitude of the resolution. In all cameras there is the extra 
complication of scattering in the photographic emulsion and, in our apparatus 
particularly, the possible additional image defects introduced by an enlarging 
eyepiece. Although in our apparatus we should not expect these latter effects 
to be serious it seems desirable that such an optical system should be tested 
in its entirety on an extended low contrast object which has the same order of 


angular scale as solar granulation. To effect this we have measured the contrast 
transmission functions of the comple 





eyepiece and 
photographic emulsion—and applied the results to the interpretation of the 
photographs. 

3.2. The contrast transmission functions.—The contrast transmission function 
of an optical system, also known as the sine wave response factor, describes 
its response to various spatial frequencies in the object plane. It is a concept 
first suggested by Duffieux (13) and developed by Schade (14) and is now used 
increasingly to describe quantitatively the performance of an optical system 
in preference to the concept of resolution measured as lines/cm or as the 
** Rayleigh limit’’. Analytically, it is the Fourier transform of the distribution 
of light in the image of a line object. An equivalent definition is that it is the 
ratio of contrast in the image plane to that in the object plane for each spatial 
frequency separately. ‘The numerical value of the transmission function varies 
with spatial frequency, reaching zero at the “cut-off” frequency of the 
system. ‘The significance of the cut-off frequency is that no greater spatial 
frequencies are passed by the optical system; if such spatial frequencies exist 
in the object they are completely suppressed in its image and cannot be restored 
in any way. At wave-length A the cut-off frequency corresponds to an angular 
resolution limit of Ad radians, where d is the linear diameter of the objective. 
If the optical system includes a photographic emulsion, it is necessary to include 
this as equivalent to an additional component. 

3.3. Measurement of the contrast transmission functions.—The measurement 
of the transmission functions is noteasy. It may be done directly by scanning 
the image of a line object and taking the Fourier transform afterwards, or by 
arranging a mechanical scanning device which will automatically measure the 
Fourier transform, as has been done by Ingelstam, Djurle and Sjogrén (15). 
This method has the advantage of measuring the phase as well as the amplitude. 
Alternatively, an extended test object may be used. For many years a square 
wave test object has been popular, but although this has the advantage of ease 
of construction it has the disadvantage of introducing uncertainty in going 
from the square wave response factor to the sine wave response factor (16). 
The use of a sinusoidal test object has been developed by Schade; for this work 
he uses film sound track (17). 

In our tests we have used sinusoidal test objects which have been prepared 
photographically, starting with an accurate transmission grating of spacing 
12 lines'mm and of square wave form. This grating has been reduced in size 
by projecting it with a high quality camera lens and photographing the slightly 
defocused image on a Kodak Maximum Resolution plate. For each of several 
reduction ratios this procedure is repeated using various exposures and different 
amounts of blurring, and the most suitable test object selected. Suitability is 
judged by two factors; contrast, or percentage modulation, and purity of 
waveform. These two characteristics are easily tested by viewing a distant 
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lamp through the grating; the amplitude of each modulation frequency is then 
proportional to the integrated brightness of the corresponding spectral order. 
Using this technique it is relatively easy to produce gratings of at least 
170 lines/mm which attain such perfection of waveform that they show no 
trace of spectra above the first order. The brightness of the first order spectrum 
relative to the zero order depends upon the depth of modulation; this property 
provides a useful guide for the rapid selection of gratings suitable for further 
tests. This spacing of 170 lines/mm is near the cut-off frequency of an f/10 
optical system. 

The modulation of each grating was measured by scanning in a micro- 
photometer at very slow speed, using a microscope objective of numerical 
aperture 0-28. Such an objective has a Rayleigh resolution of 1250 lines/mm 
and may be expected to measure the modulation of a 68 lines/mm grating to 
within 5 per cent. The modulation used increased with the grating frequency, 
being of the order of 15 per cent at the lower frequencies and up to 50 per cent 
at the higher frequencies. 

In our tests of the telescope we have been obliged to use a collimator, as a 
sufficiently distant source would have introduced poor seeing in a long light 
path near the ground. 

The collimator used was a chemically silvered paraboloid of diameter 46 cm 
and focal length 640 cm. As there is a reduction of image size by a factor x 2°3 
it was necessary to use gratings of spacing up to 68 lines/mm only. The figure 
of the paraboloid was known to be of excellent quality. During the tests careful 
observations were made of the local seeing and a test was not made until the 
steadiness of the first diffraction ring of an artificial star source showed that the 
seeing was sufficiently good. Otherwise the identical optical system of the 
flight was used and the film was developed under exactly the same conditions 
as those of the flight film. 

3.4. Results of measurements.—The results for the complete telescope are 
shown in Fig. 2. The curve for an aberration-free system in monochromatic 
light is also shown (18); the last point on this curve is the calculated cut-off 
limit of the telescope. 

This is the first time that the contrast transmission functions have been 
measured for a complete photographic telescope. It is clear from the diagram 
that in spite of the apparently satisfactory results of conventional visual tests 
made on the objective before the flight, the performance of the whole telescope 
seems markedly inferior to that of an aberration-free system; the contrast 
transmission function at a spatial frequency corresponding to 1”, for example, 
being only about one-third of its value for an ideal telescope. 

Although it is known from visual tests that its resolving power is good, 
as defined by the Rayleigh criterion, later examination in the light of the 
measured contrast transmission functions revealed spurious light in the outer 
parts of the image that had previously been thought unimportant. The only 
additional optical component is the chemically silvered paraboloid used as a 
collimator. It has been found by Lyot that such a mirror may form an image 
of a point object which contains excess light in its outer parts, the effect being 
ascribed to the variable thickness of the chemically deposited silver coating (19). 

We attribute most of the fault to the objective. Other tests suggest that its 
defects are due partly to a rapid variation of spherical aberration with wave-length 
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and partly to disturbances in the wavefront produced by local figuring of the 
surfaces. Separate tests of the eyepiece show that it is almost free from defects. 





Contrast transmission functions 
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Fic. 2.—Contrast transmission functions of complete telescope (objective, eyepiece, filter, 
photographic emulsion) as used on balloon flight (lower curve). The upper curve is the computed curve 
for an aberration-free optical system of the same size. 


4. Preliminary discussion of the appearance of the photographs 

4.1. Image quality.—During the flight the focus setting was changed by 
0-1 mm after every ten photographs. ‘Two points now arise: first, that the 
ten photographs taken at the position of “‘ best focus” are not all of equal 
quality, but one photograph in particular is somewhat superior to the others. 
Second, there is not a regular deterioration of image quality with change of 
focus: it sometimes happens that an isolated photograph taken as far as 0-3 mm 
from the “ best focus ” is of better quality than some taken nearer to this position. 
Even 2 mm from true focus there is an occasional photograph that shows a trace 
of fine structure of the order of size of the granulation. These observations 
show that the seeing is still not uniformly good even at 18 000 ft. 

This kind of effect is not easily discerned at ground level with a large 
objective because the wavefront passing through the lower atmosphere is 
deformed by disturbed regions whose lateral dimensions are smaller than the 
telescope objective, and the resulting image is always confused. On the other 
hand, at high altitude the refractive index of air is smaller and the thickness 
of a disturbed region must be greater to produce the same phase change. We 
may expect these regions to be of greater lateral extent, and their principal 
effect to be a random change in the apparent focal length of the objective, 
without a necessary deterioration of image quality. 

Whilst the flight was in progress, careful observations of the Sun were made 
at Meudon, the place of departure, by Mme d’Azambuja: the seeing throughout 
was of very poor quality. That we have obtained under these conditions 
photographs which are superior to the best previously secured from the ground 
confirms the supposition we set out to test, that there is undoubtedly a great 
improvement to be obtained by observing from a balloon at an altitude of 
about 18 ooo ft. 
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4.2. Shapes and dimensions of granules.—Plate 1 shows enlargements of 
two parts of a good photograph taken near the centre of the solar disk. The 
resolution appears to be about 0”+5; it is clearly greater than the mean size 
of the granules so that the granules appear little altered in reproduction. The 
photograph clearly shows the individuality of each granule. Their mean size 
is about 1”-4. Few granules are larger than 2” and few are smaller than 0”-8. 
When examined individually, the granules show a much less regular structure 
than that suggested by previous photographs. The granules are rarely circular. 
Their shape is often polygonal, elongated, or complex: many of them, but not 
all, are of uniform brightness with sharply defined edges. 

4.3. Large-scale variations of structure—The two photographs of Plate 2 
are of regions where the granular structure is particularly irregular. Further 
comments are given in the legend. 

4.4. Associations of granules—It turns out that some granules are not 
completely separated but touch, or merge into each other, forming amorphous 
groups of greater size. In associations they can delineate complex shapes. 
Plate 3 (first row) gives three examples; left, the grouping of the granules traces 
a ring; middle, a crescent; right, a complex filament. 

4.5. Intergranular regions.—The darker material between the granules is not 
of constant brightness. In some places granules seem to be missing, their place 
being taken by small areas whose material is generally darker than the 
neighbouring intergranular space. These are the “ pores”. Plate 3 (second 
row) shows three examples of such formations. 

4.6. Striations.—Curious alignments of the structure are sometimes revealed ; 
the granules are drawn out in lines, isolated or in parallel arrangements over a 
distance of as much as 45”. ‘These linear arrangements may be evident either 
in the granules themselves or in the dark intergranular material. Plate 3 
(bottom) shows an example, and they may also be seen in parts of Plate 1 
(top photograph). Similar structures which have been seen previously on 
other photographs have been interpreted as effects of atmospheric turbulence. 
Our photographs seem to confirm that these striae, although rather rare, are 
probably real. 

4.7. Isophotometric contours.—We have drawn curves of equal intensity for 
different values of brightness, using the isophotometer of the Sacramento Peak 
Observatory, New Mexico, placed at our disposal by the Director, Dr J. Evans. 

Three areas have been selected each of size 15” x20” and are shown on 
Plate 4. In each of these areas five tracings have been made corresponding to 
five equal steps of brightness. The left column is a region of normal granulation. 
The region in the centre column contains a dark “ pore” on the left and a 
particularly bright granule at the bottom. The right-hand column is selected 
from one of the regions referred to above which is darker and contains only 
granules smaller than the mean diameter. 

Examination of the isophotes confirms the complexity and irregularity of 
the granules. As an example the photograph at the top left shows a particularly 
large granule right of centre. The isophotes in the bottom two photographs 
show its complex outline and its internal structure of three nuclei. 


5. Measurement of autocorrelation functions of brightness—We wish to 
obtain finally the power spectrum of the granulation. To derive this we have 
first measured the one-dimensional brightness distribution by scanning across 
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the negative with a very small spot of light. The scans were each made along 
a length of 24 mm with a square spot of o-o2 mm side ( =0”-14), using a recording 
microphotometer. ‘The resulting traces were measured at intervals of 0-021 mm 
(=0":142). This value is within the interval required by the sampling theorem 
(20), which for a cut-off of o”-37 demands an interval not greater than 0”-18 
for the complete description of the microphotometer curve. With the Cambridge 
electronic computer (EDSAC II) the 1100 readings from each scan were 
converted into intensities using the calibration curve of the emulsion, and the 
autocorrelation function found for a shift of up to 100 measures, or 14”. If 
the N successive brightness values are J,, J,...J,, the autocorrelation function 
¢(j) for a shift 7 is defined as 


k=N-j 
SMM) 
¢(j )= 





k=N-j 


2 (,—-M? 
k=1 


where M is the mean brightness defined as 
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In all, six scans were made across different parts of the best negative (frame 221). 
The autocorrelation functions for each scan have been computed, and two of 
them are shown in Fig. 3 (a) and (4), plotted as a function of angular distance 
on the Sun. These curves show an initially steep descent followed by small 


fluctuations about the values 6=0-2 to 0-3. Part of the long period fluctuation 
is probably instrumental in origin, but the fluctuations do not repeat from one 
scan to another and so are chiefly of a random nature. We emphasize that we 
have not used a moving average to eliminate the effects of slow changes of 
brightness across the photograph, real or spurious; had we done so the 
autocorrelation functions would have diminished to zero. 

The “ half-width” of the initial descent is a measure of the resolution of 
the photograph. The half-widths of the curves shown in Fig. 3 are 0:50 
and o0”-72. ‘The half-width of the autocorrelation curves for the other four 
scans are 0"-75, 0-79, 0”-7g and 0”-g7. These results show that the average 
quality varies considerably over the area of the photograph. The best definition 
is along the line of scan of Fig. 3 (a) (half-width 0”-50) but this autocorrelation 
function curve is a measure only of the average definition along this line. As 
the measurements were made with a small spot no appreciable contribution 
to these curves from photographic grain may be expected. Were there a 
contribution from grain a peak at the origin of the autocorrelation curve would 
be expected. Such a peak cannot be discerned on these curves. 

Autocorrelation functions have been published by other investigators. 
Wlérick (21) has studied a photograph taken by W. A. Miller at the R.C.A. 
Laboratories with an 11-3. cm objective. ‘The half-width of Wlérick’s curve 
is 1":23. One photograph taken at Mount Wilson, “ of exceptionally fine 
definition’, has been analysed by Frenkiel and Schwarzschild (22). This 
photograph was taken with the 60 ft solar tower diaphragmed to 10 cm aperture, 
the effective wave-length being 5800 A. The half-width of their curve is 0”-g1, 
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Fic. 3.—Autocorrelation functions computed from two scans across the best negative. 
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but as the calculated cut-off of the telescope is 1”-11, it seems likely that much 
of the first part of their autocorrelation curve is instrumental in origin. Uberoi 
(23, 24) has also analysed a high quality Mount Wilson ‘photograph. This 
photograph, taken with the same 10 cm aperture and also at a wave-length of 
5800 A, yields an autocorrelation curve of half-width 0’-76. For the same 
reason, much of this curve is instrumental in origin, as the author shows. 

The significance of the autocorrelation curve found in this investigation 
will be discussed in a later paper. 


6. Measurement of contrast 


6.1. Introduction —In the measurement of contrast a large instrumental 
correction is inevitable. Even if the optical system is free from aberration 
the correction is a considerable one, while aberrations too small to be noticed 
in ordinary visual tests will increase the correction greatly. In the absence of 
information about the optical quality of a telescope and the photographic 
emulsion the magnitude of the correction is very uncertain. ‘The method usually 
used to obtain the correction is that of Wanders (25) which involves scanning 
radially across the limb. ‘This must be done on separate photographs taken 
with longer exposures and almost inevitably gives erroneous results. In 
addition, the curve given in Fig. 2 for an ideal telescope shows that a great 
loss of contrast may be expected at frequencies near the cut-off. Probably all 
the observations made so far have been undercorrected. 

6.2. Previous observations.—Two quantities characterize the brightness 
properties of the granular structure. One is the root mean square brightness 
variation, and the other is the mean contrast between granular and inter- 
granular regions. In Table I we list previous determinations of these quantities. 
In the last two columns the upper figure for each determination is the uncorrected 
direct measurement, and the lower figure, in brackets, is the value as corrected 
for instrumental effects by the individual investigators. 

We may reconsider these data in the light of the known contrast transmission 
functions for an ideal telescope, assuming a correction appropriate to the mean 
granule size. ‘This would appear to be a minimum correction if only because 
no allowance has been made for the smaller transmission function corresponding 
to the first harmonic, which may be supposed to exist if there is light-dark 
asymmetry in granulation. The correction to the observed contrast must be 
at least that indicated for an ideal telescope and in practice will be larger by 
an unknown amount which depends upon the properties of the instrument 
and the seeing. 

As an example, we take the observations of Keenan. In obtaining his 
correction he assumed a mean granule size of 1”-1. But the transfer function 
for an aberration-free objective of this size is only 0-37. This means that the 
correction must be at least a factor of x2-7, giving a contrast of at least 
28 per cent. The correction may be expected to be considerably bigger than 
this because of instrumental imperfections and seeing, and especially as Keenan 
mentions that there are no granules on his photograph smaller than 1” in 
diameter. 

The observations of ten Bruggencate and Miiller—a continuation of earlier 
work (29, 30) with the Potsdam tower telescope—gave a value that is super- 
ficially in good agreement with that of Keenan. We have, however, already 
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Above: A part of the solar surface ncar the centre of the disk. Dimensions of field 
go” X75”. 

Below: Enlargement of field in which the granulation is particularly uniform. 
Dimensions of field 48” * 40”. 
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Enlargements of regions in which the granulation is noticeably irregular. 
: he centre is a region where the granules are particularly bright. To the 
ht a dark region where the granules are of unusually low contrast. 
top right a region where the granules are small; on the left, dark 
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Top row: Examples of granules of complex shapes. 
Second row: Examples of small dark areas and “ pores’. 
Third row: Left, a region of high contrast. 

Centre, a region of normal contrast. 

Right, a region of low contrast. 
Bottom: Example of structures aligned in parallel streaks. 
(All the same scale; dimension of the squares: 12” X 12".) 


D. E. Blackwell, D. W. Dewhirst and A. Dollfus, Solar granulation observed from a manned balloon, 
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Tas_e I 


holies Refs Aperture Wave-length R.m.s. observed Contrast observed 
q : (cm) (A) (r.m.s. corrected) (contrast corrected) 


ten Bruggencate 28 12 4350 = 9 
and Miiller (18) 
Keenan 19 5550 10°3 


(15°0) 


25 
(35) 


Thiessen 60 5650 


Richardson and 31 10 5800 
Schwarzschild 

Frenkiel and 22 10 5800 
Schwarzschild 

Leighton 32 30 5850 14 

(—) 

Richardson and 31 30 6300 , - 
Schwarzschild 

Waldmeier 35 (=) (~) 


(30 to 40) 
Krat 36 24 (—) ’ _ 
(13) 
pointed out that Keenan’s value is undercorrected, and in addition the correction 
for the large wave-length difference is quite unknown. 

The Mount Wilson photograph used by Frenkiel and Schwarzschild can 
tell us little about the contrast of the small scale granulation because the cut-off 
of the telescope is at 1”*1 and even at 1"°5 the contrast transmission function 
of an ideal telescope at 5800 A is only o-15. 

The r.m.s. value of 3:2 per cent obtained by Richardson and Schwarzschild 
was measured on a spectrogram and no instrumental correction can be estimated 
because of the possibility, for example, of astigmatism in the grating. The 
r.m.s. value of 5-4 per cent derived by these authors is for an average element 
size near 2”, and as the contrast transmission function for structure of this 
size is 0-36 for an aberration-free system, this means that the r.m.s. brightness 
fluctuation is probably at least 15 per cent. 

The contrast given by Leighton is substantially a peak to peak brightness 
measurement over a scan of length 120”. 

Thiessen’s observations were made visually with the 60 cm refractor of the 
Hamburg Observatory and his value of contrast agrees well with that found 
by Waldmeier, who does not, however, give details of the method of observation. 

The work of Plaskett (37) gave a total amplitude variation of 8 per cent on 
his best plates; however, this work refers to brightness variations in a large 
scale structure with elements about 5” in diameter. 

6.3. Present results—On our photographs, taken with a 29 cm objective 
at a wave-length of 5300 A, the r.m.s. fluctuation of brightness measured on 
a microphotometer tracing is 4°6 per cent. Taking an effective transmission 
function of 0-25 appropriate to a mean granule size of 1”-4 (see Section 4.2) 
we obtain a corrected r.m.s. fluctuation of 18 per cent, without any additional 
correction for seeing. 

The same microphotometer scans show that areas with an observed contrast 
of up to 7 per cent occur frequently, while a contrast of 9 per cent is not 
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uncommon (excluding those areas where there are known sunspots or pores). 
The isophotometer tracings show that the mean contrast between granular and 
intergranular regions is 11 per cent. This value is slightly higher than that 
indicated by the microphotometer scans probably because it refers to small 
regions selected for their good definition, while the microphotometer scans 
were made at random along the larger dimensions of the photograph. For 
further discussion we take the observed mean contrast between granule and 
intergranular space to be 10 per cent. This gives the ratio of mean contrast 
to r.m.s. brightness fluctuation of 2:2, which seems reasonable. Applying the 
same correction for a contrast transmission function of 0-25 we get a corrected 
value of 40 per cent. iowever, this correcting factor is a lower limit not only 
because the seeing will reduce the effective contrast transmission function 
(although because of the nature of our observational material we suppose that 
any further correction will be small), but also because the effective width of the 
intergranular regions is probably much less than 1”+4. 

Our final result of a contrast of 40 per cent is not inconsistent with most 
of the values previously obtained by other investigators, provided that the 
observed values quoted in Table I are corrected anew by the larger factors 
indicated by the discussion in Section 6.2. However, the corrections that we 
have made to our own observed contrast are minimal. Simply interpreted 
this contrast implies a difference of temperature between mean granule and 
mean intergranular space of 520°K. 
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MAIN-SEQUENCE S'FARS 
C. B. Haselgrove and F. Hoyle 
(Received 1959 March 9) 


Summary 


Integrations have been carried out for uniform stars of Types I and II. 
It appears likely that the separation at the zero-age main-sequence of the 
different populations does not exceed 0-1 magnitude, provided a correct 
effective temperature scale is used. 





1. Introduction.—Models for main-sequence stars—i.e. stars of uniform 
chemical composition in which the outflow of energy is determined by nuclear 
processes—have been computed for the following cases: 


Type I stars, X=0°75, Y=0-24, A=8x 107, O [MZ 1250. 

Type II stars: 

Case Ila, X=0-99, Y=0-009, A=4x 10-*, O2T M=4O. 
Case IIb, X=0°75, Y=0-249, A=4x 10-*, O[> M230 

M=mass of the star, X=hydrogen concentration by mass, Y=helium 
concentration by mass, AX-!=abundance of metal atoms by number relative 
to hydrogen. 

The integrations were performed on the IBM 704. The method of 
integration was the same in principle as that described in a former paper (1), 
but with the increased speed and storage capacity available to us it was possible 
to use formulae of a more accurate kind for the opacity and for the nuclear energy 
generation. Account was taken of the appropriate surface boundary conditions. 


2. The opacity.—The opacity per gram « may be written in the form 


K=(k,1+K,)7, (1) 


where x, is the radiative opacity and «, is the conduction opacity. Formulae 
for x,, K, were constructed to give representations of the tabular values worked 
out by Keller and Meyerott for the radiative opacity (2) and by Mestel (3) for 
the conduction opacity. The following quite simple expression was found to 
represent Mestel’s values with adequate accuracy (5 per cent in the important 


range of p, T) 


T : 1-685 
oan [ os478( m X) 7a | } " 


where p is the density in gmcm~', and the temperature T is in units of 10° °K. 
The situation is much more complicated for the Keller-Meyerott tables. Because 
of the absorption edges of the various elements, the values do not run in any 
very simple way, so that empirical formulae must of necessity appear awkward. 
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The following were found to represent the tabular values within an accuracy 
of ~ 10 per cent in nearly all cases: 
o1<T<10, «,=o0-19(1+X)+Ep(1 +X), 
10<T<20, x,=0-19(1+X)+£p(1+X)+p(1-o-'1T)(Er+X—L), (3) 
20<T, x,=0'19(1+X)+p, 
the units for p, T being the same as in (2). The quantities €, ¢ are given by 


bat (1-—X-—Y)(1-o0°05T) | _.- P(t +X) ) 
£=6°5 x 108 T? 42-57% ne GU Kose, ae 





X Y , 
+415 x 0! AR 7+ 2507"! + sexp(—|$-28737))) | 





225—155X —190Y 
T35 P 
These formulae enable «x to be computed when X, Y, p, and T are specified. 

Considerable simplifications could be effected if no attempt were made to 
represent the tabular values over the whole of the ranges of p and T. If 
separate formulae were used in different parts of these ranges no doubt each 
formula could be given a straightforward structure. It must be emphasized, 
however, that unless care were taken to preserve strict continuity at the joins 
of the separate formulae severe difficulties in the iterative procedure for the 
numerical solutions would be encountered sooner or later (at any rate in the 
computational method we have used). The present formulae do preserve 
continuity. 

The question arises as to how far the formulae can safely be extrapolated 
outside the range of the tables from which they were constructed. At low 
densities or high temperatures such an extrapolation can be made, since in these 
cases the formulae have a correct asymptotic form. But extrapolation at high 
density and low temperature would be most unsafe. The limitation is important, 
since the latter case arises for dwarf states with M< ~o-5@. The present series 
of integrations were in fact carried down to M~o-2©. Results are not given 
below at this low mass, however, since doubts about the opacity make them of 
uncertain value. 

Nor can the above formulae be used for the low values of T occurring at and 
near the photosphere. Here a quite different procedure must be followed. 
When p, T, X, Y, A, a;, b;, ¢; (i=1, 2, ...) are specified, the electron pressure P, 
and the quantities 0, x, y, 2, u; (i=1, 2,...) are determined by the following 
system of equations: 


P,=8pTO, (5) 
log x = — 13540 + § log T—0-48 — log P,, 


y 
4= 








log y= —24'486 + $ log T + 0-12 —log P,, 
logz = — 54'170+ $ log T— 0-48 — log P., 
logu;= —a,9 +, 3 log T+ 5;,—log P,, 


_ 5040 
O= >? 
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Xx Yy(1 +22) 
r+x 4(t+y+y3) — 
where ® is the gas constant, the logarithms are to base 10, and T is now in °K. 
The physical meanings of the symbols appearing in these equations are: 

x= (ionized H density)/(neutral H density), 

y = (singly ionized He density)/(neutral He density), 

z= (doubly ionized He density)/(singly ionized He density), 

u; = (singly ionized density of constituent?)/(neutral density of constituentz), 


a;= ionization potential of constituent 7 (eV), 


Q= ASS (7) 


ym , 
Ti+; 


b;= weight factor for constituent 1, 
¢; is proportional to the number density of constituent 7, the constant of 
proportionality being chosen so that }c;=1, 
1 
4X-1- >; (number density of constituent 7) 
in (number density of hydrogen) 





Multi-stage ionization of constituent i has been neglected for all 7 in 
equation (7), © having the simple interpretation 


©/m,,= number of electrons per gram of material, 


where m, is the mass of the hydrogen atom. This approximation can be made 
because the contribution to © of all elements except hydrogen and helium is 
always very small unless x, y, <1, and these inequalities require T to be low, 
in which case multi-stage ionization can be neglected. Indeed at low T we can 
neglect all elements except the metals, and of these only magnesium, silicon, 
iron, aluminium, sodium, and potassium are important. The above summation 
over i can therefore be taken as referring to these six metals only. The appropriate 
ionization potentials a,, weight factors b;, and relative abundances c¢; are 


a; (eV) b; Cj 
Magnesium 7°61 0°04 0°3750 
Silicon 8-11 —0°42 0°3450 
Iron 7°86 —o0'08 0'2250 
Aluminium 5°96 —0'98 0°0370 
Sodium 5°12 —0°64 0'0167 
Potassium 4°32 —0'92 00013 


The relative abundances have been taken from Suess and Urey (4). The 
quantity A in equation (7) represents the overall abundance of the six metals. 
The tables of Suess and Urey give A~8 x 10~°, and this can be regarded as 
typical for Type I stars. The value for Type II stars was taken to be twenty 
times less, viz. A= 4x 107°, 

A satisfactory approximation for the opacity is given by 


. XP, 0°870 X exp [15-94 — 1161/7] 
= 0280 xp | —6-278 704), Pits 
erin aegis aT? | iis ] (1+x)(1+0°057*) 


+ arial exp| 13°81 ~ | +y exp| 13°81 - =|) (8) 
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T in units of 10*°K. The first term on the right-hand side of (8) gives the effect 
of electron scattering, the second term gives the contribution of the negative 
hydrogen ion, the third is the photoelectric contribution of neutral hydrogen, 
and the remaining term in Y represents the combined effect of neutral and singly 
ionized helium. Formula (8) was used at all temperatures below 10°°K. No 
attempt was made to preserve strict continuity between (3) and (8). Special 
provision was made to ensure that no computational difficulty could arise from 
a lack of strict continuity in this particular case. 


3. The nuclear energy generation.—The following formula was used for the 
energy production 


8 \2 

10 °*pX?T* + pX(1-X— (=) erg per gm, (9) 
T in 10°°K. The first term represents the contribution of the P—P chain, and 
the second that of the C—N cycle. 

The numerical coefficient used for the P-P chain is greater by a factor 1°5 
than the rate given by Burbidge, Burbidge, Fowler and Hoyle (BBFH, §). 
This increase takes account of the recent radical change in ideas concerning the 
mode of completion of the P-P chain, viz. 


He3(Het, y)Be7(e-, v)Li?(p, y)Het 


|-(p,y)B*(B+, v)Be**(x)Het, 

competing strongly with He*(He*,2p)He* whenever the concentration of Het 
is appreciable and the temperature is above 13m.°K. According to Fowler (6) 
the present increase of the rate by 1-5 should give a good approximation to the 
effect of the somewhat complicated competition between these different reactions, 
provided that we are not concerned with dwarf stars appreciably fainter than the Sun. 
The rate falls off at temperatures below 13 m. °K, at first slowly, and then rather 
steeply for temperatures below 8 m. °K, because the P-P chain then ends at He’®. 
The effect is to reduce the energy production rate by a factor ~3 in faint dwarfs 
with M~ 0/3. 

The coefficient chosen for the C—N cycle rate takes into consideration the 
enhanced rate of the O'*(p,y)F! reaction, as discussed by BBFH. This has 
the effect of augmenting the nitrogen abundance through F!"(8+, v)O"(p,a)N™. 
The carbon and nitrogen abundances were estimated from the tables of Suess 
and Urey (4) and the prescription of adding half the O' abundance to the 
N!* abundance was used in order to take reasonably satisfactory account of these 
conversion reactions. It is true that this procedure exaggerates the coefficient 
of the C-N cycle for stars with central temperatures < 20m. °K, but in such stars 
the main contribution to the energy generation comes in any case from the P-P 
chain. The situation is that the coefficient must be increased in all stars in 
which the C-N cycle is important. 

In accordance with this prescription the concentration by mass Xpy of 
the group of elements that promote the C-N cycle is given by 

X71Xoy = 12C + 14(N+0°50), 


where C, N, O are abundances of C!, N'4, O1* by number relative to hydrogen. 
We can rewrite: 





, _[ 12C+14(N+0°50)], » 4, 
Xou=| XG —X-¥) (1-X-Y). 
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The point of expressing X,y in this way is that we expect the absolute concentration 
of carbon, nitrogen, and oxygen in primaeval stellar material to be proportional 
to 1—X— Y, so that the concentrations relative to hydrogen are proportional to 
X~1(1—X-—Y). The value of the quantity in the square brackets should therefore 
be the same in all primaeval samples of stellar material. The value of this 
constant quantity was estimated from the tables of Suess and Urey, so that 


Xoy = (known quantity)(1—X— Y). (10) 


Assuming that the non-resonant reaction N'4(p,y)O™ is the determining 
step of the cycle, the formula for the rate of C-N energy production is (cf. (5)) 


7°86 x 107" pXXoyfoyT~** exp (— 152°3T-**), (11) 


T in 10° °K, where f,y is an electronic shielding factor ~ 1-5 under conditions 
where the C-N cycle is important. For 20<7'<45, the range of J under 
consideration in this paper, (11) can be well represented by 


(known constant)pXX,,7"*. (12) 


Combining (12) with (10) yields the second term in (9). 

It is particularly to be noticed that when written in this way the contribution 
of the C—N cycle needs no change when we pass from Type I to Type II stars. 
The proportionality on 1—X— Y ensures that the coefficient of the C—N cycle 
changes in a correct way when X + Y changes. For Type I stars 1—-X— Y2o-o1, 
whereas for Type II stars 1—-X—Y is probably <o-oor. 


4. The surface condition.—The method described by Hoyle and Schwarzschild 
was used (7), except that the differential equations were integrated numerically 

i.e. the approximate integrals of Hoyle and Schwarzschild were not employed. 
A slight change was also made in the formulae used for the convective energy 
flux. In place of formula (84) of Hoyle and Schwarzschild 


, 282)... en 


was used, again with 8 =0°3. 


5. The homology relations.—Although homology relations cannot be used 
for determining strict numerical results, they are of value in showing the 
approximate dependence of the radius and luminosity on the opacity and 
energy generation coefficients, and on the mass and mean molecular weight. 
The relevant formulae were obtained by Hoyle and Lyttleton (8) and are: 


Ys 1 2n+6 l4n+45 10n+31 ) 
L pt x E 2n+5 K 2n+5 bh 2n+9 M 2n+5 , 
I 
2 2n—15 20-7 ( 4) 
J 


R, cE %+5 «x ~ Int pp 2nts Mints, 


where L,,, R,, u represent the bolometric luminosity, radius, and mean molecular 
weight, where the opacity per unit mass approximates to the Kramer’s form 
xp/T*°, and where the nuclear energy generation has the form Ep7” per unit 
mass per unit time. 
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When on the other hand the opacity is due mainly to electron scattering the 
corresponding expressions are: 


Lng © eM, 


n—4 n—1 (15) 


1 
R, cEa+3 pats Mats, 


6. Results —The relevant quantities for determining a stellar model (with 
composition specified) are L,,, R,, the central pressure P, and the central 
temperature T,. Results for these quantities, and also for the effective surface 


temperature 7,, are given in Tables 1, 2, and 3. The surface temperatures 
were obtained from 


Ly = 7acR?T,', 
where a is Stefan’s constant and c is the velocity of light. 


TABLE 1 
Type I stars; X=0°75, Y=0'24, A=8x 1075 (Mo =1'985, Lo =3°780, Ro =0°6951) 
Mass L pol 5 A TZ; R, 
(10*8 gm) (10* erg sec~') (°K) (10!* dyne cm~?) (10° °K) (101! cm) 
247°78 6-9899 x 10° 57400 2°O195 44°746 9°5184 
165-19 3°5443 10° 54250 2°1175 43°365 7°5875 
110°13 1°6694 x 10° 50320 2°2819 41°806 6°0517 
73°417 771938 x 10% 45670 2°5543 40°046 4°8228 
59°945 4°5490 x 10° 43130 2°7529 39°086 4°3012 
39°963 16792 x 10° 37760 3°3479 37°003 3°4081 
26°642 5°5664 x 104 32290 4°3685 34°709 2°6839 
17°761 1°6594 x 104 26950 6°1275 32°225 2°1043 
11°841 4°4020 x 108 21850 Q*IOI9 29°576 1°6479 
"7384 1°2027 X 10° 18410 15°047 27686 1°2141 
+8342 7°4701 X 10? 16960 17°090 26-698 1°1268 
"7407 38363 x 10° 15090 20°197 25°349 10204 
"9037 7°7211 X10 IIIIO 25°900 21°631 0°84462 
0168 2°3215 X10 8640 26046 18°457 0°76553 
8976 1°8930 X IO 8240 25°054 17°795 0°75986 
7817 1°5332 X10 7840 23°886 17°124 0°75483 
"5592 9°8896 7030 21°283 15°838 0°75482 
"3545 6°3535 6565 18-723 14°667 0°69379 
1699 5°1202 6390 16°462 13°632 0°62352 
9963 2°6958 5880 14°480 12°687 0°56299 


7 
6 
5 
3 
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TABLE 2 
Type Ila; X=0-99, Y=0°009, A=4 10-* 
Lau i P, Yi R, 

(10* erg sec~) (°K) (101* dyne cm~?) (10® °K) 
5°7808 x 10? 16720 28°147 27°836 
*3822 x 10? 15100 30°074 26°154 
"7033 X 10? 13080 30°098 23692 
1*2034 X 107 12095 29°026 22°367 
83758 x 10 11120 27°400 20°992 
3°7543 X 10 9200 23°103 18°143 
16009 X 10 7500 18-613 15°526 
1'0516 X 10 6820 16-603 14°402 
6°8727 6410 14°745 13°362 
4°4923 6060 13°072 12°408 
2*9802 5730 11°631 11°560 
2°4401 5580 10°988 11°169 
1°3547 5095 9°3061 10°095 
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TABLE 3 
Type IIb; X=0-75, Y=0-'249, A=4x 10-* 


Mass L pot as ¥. Za R, 
(10% gm) (10** erg sec~) (°K) (10'* dyne cm~*) (10° °K) (1014 cm) 
60614 5°5458 x 10° 17800 29°818 29°083 088171 
‘8260 *3544 X 107 15220 34°884 26°730 0°78574 
8423 *3541 X10 12580 36°128 23 fo2 


0°72440 
2014 


"1533 X10 7000 32°391 20°043 0°70036 
0170 "1409 X 10 6950 30°563 19°630 0°69539 
8337 2°3205 X10 6910 28°497 18°575 0°69086 
6523 6743 X10 6865 26°268 17°496 0°68653 
3128 8-3721 6860 21°813 15°443 0°68642 
0168 4°1327 6370 17°870 13°640 0°59426 
9579 "5526 6240 17°ra1 13°289 0°57322 
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From the point of view of stellar structure, low metal content and presumably 
low carbon, nitrogen, and oxygen content are the essential characteristics of 
Type II stars. In the Type II models integrated, the metal content was reduced 
to 1/20 of that of the Type I stars and the carbon, nitrogen, and oxygen content 
was reduced to 1/10. It is uncertain whether the helium content should also be 
reduced. ‘I'wo cases were therefore considered: the more likely case where the 
helium content was also very markedly lowered, and a less likely case in which 
the Type II helium concentration was taken to be the same as for Type I. These 
cases are shown in Tables 2 and 3 respectively and are described as Ila and IIb. 


a 


7. Discussion. —The dependence of L,,, on the mass M changes quite 
markedly over the range given in Table 1 for the Type I stars. For M in the 
neighbourhood of © we have L,,0°~M*; for M~2© the dependence is 
already reduced to L,,0°0~M*; for 30<M<30© we have L,,0~M!; 
and for 30© <M < 1256 radiation pressure reduces the dependence still further 
to L,,,« ~M?. Apart from the latter effect of radiation pressure, these changes 
are determined by the opacity which alters from an essentially Kramers type of 
dependence at M~ © to an electron scattering dependence for M>30. 

The mass range of the integrations for Type II stars is not large enough to 
show such marked opacity effects. Indeed Loc M‘ is a tolerable approximation 
for these stars over the whole of the mass ranges given in Tables 2 and 3. The 
reduced effect of the metals weakens the Kramers type of dependence at M~ ©. 

Although there is no precise value of M at which the nuclear energy 
generation may be said to change from the P-P chain to the C—N cycle, as a 
general statement it can be said that the P-P chain dominates in stars with 
T,<21m. °K (Type I), 7,<22m. °K (Type II). Maximum values of P, occur 
at these temperatures. Jt is also at these temperatures that convective cores first 
appear. 


Reference to (13) shows that when the Kramers type of opacity is dominant 


R,x<M Sn+5. Now when Kramers opacity is dominant M is sufficiently small 
for the P-P chain to be more important than the C-N cycle. Thus we must 
put 7~ 4, so that R,oc ~ M"8—1.e. the radius depends only very weakly on the 
mass. This effect is shown very clearly in Table 1, where the values of R, decrease 
with M until a value ~ 0-75 is reached. Then R,~0-75 is maintained very closely 
for a while as M decreases. 
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A further consideration of the dependence of R, on M shows, however, that 
R, does not continue to stay nearly constant as M continues to decrease. Thus R, 
declines markedly below o-75 for the last three entries of Table 1. This behaviour 
cannot be understood in terms of the simple homology relations. It arises from 
the surface boundary condition. As 7, decreases with decreasing M the ionization 
of hydrogen at the photosphere also declines. This causes a marked decrease 
of photospheric opacity as T, falls below 7000 °K. The consequent rise of the 
photospheric density produces an important effect on the structure of the outer 
regions of the star, causing the surface zone of convective instability to penetrate 
downwards to deeper and deeper levels. The decline in R, due to this cause was 
first considered by Osterbrock (9). 


8. Accuracy of the integrations.—For computations performed on a digital 
computer, questions of purely numerical accuracy scarcely arise. Accuracy is 
determined by the physical approximations that have been introduced in setting 
up the problem. 

The Sun provides a valuable check on the physical theory used for stars on 
or near the main-sequence. The mass, luminosity and radius are known. Since 
moreover there is no reason to believe that mixing has ever taken place between 
the central energy producing regions of the Sun and the surface, spectroscopic 
analysis of the outer solar layers presumably determines the values of X, Y 
and of the carbon, nitrogen, and oxygen concentrations as they existed within 
the primaeval solar material. 

Against these advantages, the Sun is not a ‘‘zero-age’’ main-sequence star. 
Geophysical evidence suggests that the Sun is probably about 5 x 10° years old, 
During this time the luminosity, radius, and the composition of the central 
regions must have changed appreciably. To allow for this, an evolutionary run 
was obtained by one of the present authors. This was done as a part of the 
investigation of the ages of Type I and II stars (cf. following paper). The results 
will be given here, however, since they are of more immediate application to the 
present paper than to the following one. The actual present-day solar values 
are given in brackets in the last line of Table 4. The Kelvin-Helmholtz 
contraction time has not been included in the age estimates. 


TABLE 4 
M=1'985 x 10" gm, X=0°75, Y=0'24, A=8x107%. 
Evolution of the Sun 
L poi (10* erg sec) R, (10! cm) Age (10° years) 
2°620 0*5630 0000 
2°733 0°5665 0*700 
2°919 0°5740 1°75° 
3°213 0°5875 3°150 
4°016 0°6310 5°950 
(3°780) (0°6951) (~ 5"000) 


Interpolation between the last two entries yields an age ~5-2 x 10° years 
at L,.=3°78 x 10% erg sec~!, in excellent agreement with the geophysically 
estimated age. It may be noted that no allowance has been made for the loss of 
energy by neutrino emission from B® (cf. following paper). According to 
Fowler (6), this effect is not likely to be important for stars with internal 
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temperatures < ~13m.°K, but might well amount to a 3 or 4 per cent effect 
—in the sense that the calculated age must be reduced to this extent. 

This age agreement means that, subject to the ratio Y/X having been chosen 
correctly, the opacity values used in the present calculations must also be correct. 
In this connection it may be noted that an opacity error can be compensated 
(falsely !) in the luminosity calculation by a corresponding error in Y/X. Thus 
a calculation in which Y/X is taken as a free parameter scarcely has much value 
—one error may be cancelling another. Inthe present case the Y/X ratio indicated 
by the spectroscopic data was used, however. 

The situation concerning the opacity is less complicated in the other stars we 
have considered than it is in the Sun. Hence the opacity values in the other stars 
are expected to be at least as reliable as they are for the Sun—and indeed the latter 
seem to be correct within a few per cent. 

It is of interest that the luminosity of the Sun has increased by about 0-4 magn. 
during geologic time. By hand computation Schwarzschild, Howard and Harm 
obtained an increase of 0-5 magn. (10). 

The position concerning R, is more complicated. If the convective surface 
boundary condition had not been used in the present calculations, the value 
of R, at zero-age would have been close to 0°75, in accordance with the simple 
homology conditions discussed above (i.e. if we had used the old-fashioned 
boundary condition p, T->o at the surface, R, would have been ~o-75 at 
zero-age). Subsequent evolution would have again increased R, by about 
IO per cent—i.e. to ~o0°83. Thus the boundary condition formerly used (by 
Eddington for example) would have given a present-day radius some 20 per cent 
toohigh. By introducing a physically realistic boundary condition we have reduced 


the radius but have apparently overswung the correct value by about 10 per 
cent. 


There are two ways of viewing the situation. Inthe broad sense it is satisfactory 
that the same rather simple analysis of the surface zone, found by Hoyle and 
Schwarzschild to be appropriate for giant stars (7), should also give tolerable 
results for dwarf stars. On this basis it appears that the analysis can probably 
be applied over the whole H-R diagram provided very accurate results are not 
required. 

If really accurate values are required, the matter needs closer attention 
however. One possibility is to reduce the convective efficiency of material in 
the sub-photospheric levels. Thus the factor 8 in (13) could be reduced from 
the value 0-3, used by Hoyle and Schwarzschild for the giants, until R, came 
into agreement with the required solar value. This in fact was the procedure 
used by Schwarzschild, Howard and Harm (10), and more recently by Sears (11). 

A second way of resolving the discrepancy lies in an adjustment of the energy 
generation constant for the P—P chain. Reference to (14) shows that for 7=4, 
R,«xE*3, The radius can therefore be increased by 10 per cent through an 
increase of 2 in E—i.e. by a doubling of the rate of energy production at fixed 
p, T, X. In this connection it may be pointed out that within the last two years 
the nuclear energy estimates have in fact changed by factors of this general order 
on two occasions. In these circumstances it seems somewhat dangerous to take 
up the error in the calculated value of R, by an adjustment of the surface 
condition. Until the nuclear situation becomes certain beyond any possible 
question it has seemed to us that the original value of the convection efficiency 
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should be adhered to, and that the error in R, should be left in abeyance for the 
time being. 

From one important point of view it would be preferable if it should ultimately 
prove to be the P-P energy generation rate that still is in error. The lithium 
abundance observed at the solar photosphere is extremely low. This has been 
interpreted by Greenstein and Richardson (12) to mean that at some stage during 
the history of the Sun the sub-photospheric convection zone penetrated inwards 
to regions where the temperature was as high as 2m.°K. On the other hand 
Greenstein and Tandberg-Hanssen (13) find the solar beryllium abundance'to be 
normal, implying that convection never extended inwards to a temperature of 
3m. °K (lithium is destroyed by nuclear reactions at 2m.°K and beryllium 
at 3m. °K). 

Now convection turns out to be deepest at the beginning of the solar 
evolution. At the last solution of Table 4, convection extends inwards to 
T~1-5m. °K, whereas at the first solution (zero-age) convection extended to 
T~2m. °K, exactly as is required by the lithium—beryllium data. In contrast 
to this, Schwarzschild, Howard and Harm found that after weakening the surface 
convective transport sufficiently to secure agreement in R, the convection zone 
only extended down to 7~1m. °K, even at zero-age. These authors suggested 
that lithium may have been destroyed during the Kelvin-Helmholtz contraction 
that preceded the main-sequence era. 


g. The main-sequence as a calibration standard.—Fitting to the lower part 
of the main-sequence provides a powerful tool for determining the distance 
moduli of distant star groups. The procedure is to secure a standard zero-age 
main-sequence from nearby stars that are either young enough or faint enough 
to be regarded as essentially unchanged by evolution. Faint stars of large 
trigonometrical parallax delineate the lowest part of the main-sequence, as also 
do the stars of nearby clusters when a group parallax can be determined (as in 
the case of the Hyades). Young clusters yield the higher portions of the 
main-sequence simply by fitting at the faint end of the main-sequence to that 
given by the trigonometric stars. The distance modulus of any star group can 
then be determined by a simple superposition of its main-sequence on to the 
standard sequence (provided that main-sequence stars within the group can be 
resolved to a magnitude faint enough for them to be regarded as essentially 
unevolved). 

This method would plainly be seriously upset if the zero-age main-sequence 
were sensitive to the composition of the stars. The degree to which the 
main-sequence is changed by quite drastic alterations of composition is shown 
in Fig. 1. Because of uncertainty in the value of Y for Type II stars both cases 
have been given: Ila for Y=o-009, IIb for Y=o-249. The very satisfactory 
result emerges that irrespective of this uncertainty the main-sequence for Types I 
and II stars cannot differ by more than 0-25 magnitude. Much greater differences 
have been claimed by observers. Indeed sub-dwarf Type II sequences as much 
as 1 magnitude below the Type I sequence have been described in the literature. 
It now seems certain, however, that errors have been made in the estimation 
of 7, due to line-blanketing effects (14). 

The position concerning the coincidence of the main-sequences is probably 
much better than is suggested by Fig. 1. Although observational uncertainty 
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exists concerning the helium abundance of Type II stars it appears very probable 
on theoretical grounds (cf. following paper) that the helium abundance was 
initially very low. ‘This suggests we should adopt case Ila. Now by an 
oversight (g) was used for the energy generation rate in all cases. ‘This means 
that in case Ila the P-P energy production rate was overestimated by a factor 1°5, 
since Be’ is not formed in the absence of He*. The situation is easily corrected, 
however, since the resulting changes are small, and for small changes the homology 
relations (14) can be used. ‘These relations show that for stars in which the 
P—P chain is dominant (y= 4) the radius R, is proportional to the 2/13 power 
of the P-P energy generation rate. Thus the values of R, in Table II should be 


DEPENDENCE OF MAIN SEQUENCE ON 
COMPOSITION 
I X*0-75, Y#0-24, A= 81075 
Ila X = 0:99, ¥20:009.A*4x107° 
Ib X2 0:75, Y= 0'249, A=4x107° 
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decreased by 1-5? in all except the first four solutions (the C-N cycle dominates 
in these). Fig. 2 gives a revised plot incorporating this correction. ‘The 
important result now emerges: 

The separation of the lower main-sequences of Type I and II stars is not likely 
to exceed 0-1 magn. 

On this basis the main-sequence can be used as an excellent calibration 
standard, provided line-blanketing effects are properly allowed for. It is 
particularly to be noted that a change of the surface boundary condition or of 
the overall P-P energy generation rate (cf. Section 8) would not change this 
result since both types of star would be affected in the same way. 
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THE AGES OF TYPE | AND TYPE II SUBGIANTS 
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Summary 
Evolutionary sequences corresponding closely to the lower portions of 
(i) the globular cluster sequences, 
(ii) the M67 sequence, 
have been worked out with the aid of the IBM 704 computer. These lower 
portions of the sequences determine the ages of the Type I and Type II stars. 
By using observed magnitudes for certain nearby Type I stars of large 
trigonometrical parallax, together with the condition that Type II stars 
be older than Type I stars, it is concluded: 


(1) that the cluster-type variables cannot be brighter than My= +0°6, 

(2) that the age of the Galaxy must be greater than 10!” years. 

Evidence is also given that the initial heltum content of Type II stars was 
low. 





1. Introduction.—On the basis of hand computation Hoyle and Schwarzschild 
(1) obtained an age of 6-2 x 10° years for the Type II stars. This value was 
confirmed by Haselgrove and Hoyle (2), who obtained 6-5 x 10* years from 
calculations carried out on a digital computer (EDSAC 1). These estimates 
were based on C!*(p,y)N™ as the determining reaction of the carbon—nitrogen 
cycle. It now appears, however, that the non-resonant reaction N'4(p, y)O” 
must be taken as the determining step of the cycle (3), and that the overall C-N 
energy production rate (at given XY, X¢y, p, T) must be taken lower than the 
rates assumed in former work by a factor 10 to 10? (cf. preceding paper). The 
question accordingly arises as to whether the former age estimates stand in need 
of appreciable revision. 

The lowering of the C-N energy generation rate enhances the importance 
of the P-P chain. It turns out that during the early stages of evolution of stars 
in the mass range © to 1-5© (the range in question for Type II stars) the energy 
generation is essentially due to the P-P chain. ‘This introduces an important 
increase in the time required for evolution to carry stars away from the main- 
sequence—an increase that comes about for the following reason. 

Because of the much lower temperature sensitivity of the P—P chain (~ T* as 
compared with ~ 7" for the C-N cycle) hydrogen is burned more uniformly 
throughout the central regions of a main-sequence star. Hence dX/dX, is 
appreciably larger for the P-P chain than it is for the C-N cycle, where X is the 
mean hydrogen concentration for the whole star, and X, is the concentration 
at the centre. Now soon after X, declines to zero a star always moves to the 
right in the H-R diagram away from the main-sequence. At this evolutionary 
stage, XY must be appreciably less when the P-P is dominant than it would be 
if the C—N cycle were dominant. At a given luminosity level this greater decline 
of X requires a greater length of time to elapse before the star departs from the 
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main-sequence. Thus on simple grounds it can be seen that stars evolving off 
the main-sequence at a given luminosity level must have a greater age if the energy 
is mainly generated by the P—P chain than they would have if the energy were 
mainly generated by the C-N cycle. 

2. Results.—A sequence of evolutionary models was determined for a star 
of mass 2°6856 x 10% gm, with initial hydrogen concentration 0-99, initial helium 
concentration 0-og, and metal concentration A = 4 x 10~* (for the meaning of A, 
see preceding paper). ‘The formulae used for the opacity, the energy generation, 
and for the surface condition were the same as in the preceding paper. 

The important parameters determining the solutions are given in Table 1. 
The relation of these solutions to the main-sequence is shown in Fig. 1. The 
way in which the hydrogen concentration changes from one solution to another 
is given in Fig. 2. The abscissa in this figure, M,/M, changes from zero at the 
centre of the star (M,=0) to unity at the surface (M,=™M), M, being the mass 
interior to a sphere of radius r (<.R,). Curves of M,/M are given in Fig. 2 only 
for those solutions marked by an asterisk in Table 1. It has not been thought 
necessary to mark which particular X-curve belongs to which solution since the 
solutions run down the table in the sense of decreasing X, so that their respective 
association with the curves is obvious. 


TABLE 1 
Evolution of Type II star 
M=2:6856 x 10% gm; initially X=0-99, Y=o’009, A=4x 107% 


L pot A F. qe Age R, 
(10 erg sec") (K) (10%* dyne cm~*) (10° °K) (10" sec) (10! cm) 


4°4920* 6065 12-986 12°378 
5°0088* 6135 14°646 12°901 
5°6538* 6210 16-870 13°523 
6-5068* 6290¢ 20°117 14°312 
0693 6335 22°453 14°814 
3-1542* 6410 27°468 15°751 
9°0953 6460 32°456 16°540 
10°004 6505 37°844 17°289 
10°721* 6535 42°563 17°875 
12°473 6600 56-149 19°326 0°96186 
14°277* 6650 72°803 21°023 1°0133 
17°193* 6665 103°84 25°024 , I*1077 
19‘O7I 6655 176°52 22°474 ° 1‘1690 
22°468* 6550 325°82 23°782 : 1°3095 
24°525 6425 540°09 24°431 ’ 1°4223 
26°363* 6250 1007°5 25°32! ; 1°5599 
27°720 5975 2437°8 27°128 , 1°7499 
27°657 5630 8438-9 30°700 ‘ 19677 
30°619* 5490 24911 34°059 ‘ 2°2544 


0°68387 
"5189 0°70523 
037 0°73118 
"5567 0°76434 
‘8162 0°78549 
"2054 0°82468 
4648 0°85700 
6594 088676 
*7891 0°90942 
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The curves of Fig. 2 demonstrate quantitatively the point that was mentioned 
at the end of the previous section. It is immediately clear that at the stage where 
X, declines to zero X has declined from the initial 0-99 to less than o-85—i.e. 
some 15 per cent of the original hydrogen has been burned. This may be compared 


10* 
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with the integrations of Haselgrove and Hoyle (2) which gave a burning of about 
11 per cent at the corresponding evolutionary phase. 


EVOLUTION OF TYPE IT STAR 
x =099, Y=0:009,A=4x107° 


ZERO AGE TYPE 


MAIN - SEQUENCE 


a Log Leo (10%erg sec™') 
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The last solution given in Table 1 has a well-marked significance. It lies at 
the second turning point of the ‘Type II evolutionary sequence; the point where 
the sequence mounts parallel to the magnitude axis of the H—R diagram at a 
colour index of about +0-6. It is emphasized that this solution can readily be 
recognized in performing the integrations because at this point the surface 
condition first begins to exercise a dominating influence on the structure of the 
star(1). It might also be added that even if an error were made in the identification 
of this particular point the effect on the calculated age would be very small indeed. 
The values given in the age column of Table 1 show that the main contribution 
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to the age comes overwhelmingly from the early stages when the star lies very 
near the main-sequence. 

According to Sandage (4) this second turning point lies about 2-1 magn. 
below the cluster-type variables (RR Lyrae stars). Thus since the bolometric 
luminosity at the last solution of Table 1 is +24 it follows that the integrations 
are applicable to a Type II sequence with RR Lyrae stars at bolometric absolute 
magnitude +0°3. The relevant age is 1:23x10! years. The bolometric 
correction to be used for the RR Lyrae stars (surface temperature ~ 8000 °K) 
is about 0-2 magn. We can conclude therefore that for a sequence with 

My= +05 
the age ~ 1-2 x 10!° years. 

The absolute visual magnitude of the RR Lyrae stars has in the past been 
set, somewhat conventionally, at o-o. Evidence has been put forward, however, 
that +0°5 is a more correct value (5,6). Indeed recent work by Eggen and 
Sandage requires these stars to be fainter still by about 0-2 magn. (7). 

If we assume that for small changes of the starting point on the main-sequence 
the whole evolutionary sequence is changed homologously, it is easy to estimate 
the age for other values of My at the RR Lyrae ‘‘gap’’. Under homologous 
conditions, age is simply proportional to M/L,;, where L; is the initial main- 
sequence luminosity. For Type II stars on the main-sequence L,oc ~ M* (cf. 
preceding paper) so that the age t varies as L;°'°, Under homologous conditions 
there must also be a constant ratio between the starting luminosity L; and the 
luminosity of the RR Lyrae variables, so that My, for these stars is given by 

My = —2°5 logy) L;+ constant. 
Thus 
logy)t = — 0°75 log,) L;+ constant = 0°3/M,, + constant. 


Combining this equation with the age result. obtained above, viz. 
t= 1:23 x 10" years at M, = +0°5, we have 


lee (sag x70) =0'3(My—o'5), t= 10+0-9%V years, (1) 


The age estimate determined by Haselgrove and Hoyle (2) referred 
to M,=o-0. It is seen that with this same value of My we now obtain 
t=8-7 x 10° years. This gives the extent of the correction introduced by the 
lower efficiency of the C-N cycle. Jt is particularly to be noticed that these age 
estimates assume no energy loss by neutrino emission from the decay of B® (cf. preceding 
paper). This question will be discussed later. 

3. The age of M67.—The significance of the comparison between the 
evolutionary sequences of the globular cluster stars (Type II) and of the stars 
of M67 (Type I) has been discussed by Sandage (4). A set of integrations 
(M =2:1699 gm, X=0-75, Y=0'24, A=8x 107°) giving a theoretical sequence 
close to the observed sequence of M67 was computed. Results are set out in 
Table 2. The relation of these solutions to the Type I main-sequence is shown 
in Fig. 3. 

At the last solution of Table 2 the bolometric magnitude is +2°67. This 
agrees quite closely with Sandage’s value M,,,= + 2:9 taken at the corresponding 
surface temperature of the M67 sequence. The age is 9°23 x 10° years. 

We again assume that Type I stars with slightly different values of the initial 
main-sequence luminosity L; follow homologous evolutionary sequences. Writing 
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M,,, for the bolometric magnitude at any particular stage of the evolution we have 
M gq = — 2°5 logy) L;+ constant. (2) 


The constant in this equation depends of course on the particular evolutionary 
stage that happens to be chosen, but the important point emerges from Table 2 
that the constant is very insensitive to the evolutionary stage provided that it 
is chosen at the almost constant luminosity level at which the evolutionary track 
moves markedly to the right in the H-R diagram. 

We can again argue that the age ¢ is proportional to M/L,;. For Type I stars 
with masses of solar order we have L,oc ~ M® (cf. preceding paper) so that 
txL;°*. Using equation (2) this proportionality becomes 


t= (constant) . 10°32 Bot, (3) 
The constant here is also very insensitive to the precise stage of the evolution 
7 precise stag 
provided that it is chosen at the almost constant luminosity level of Table 2. 
Thus since we have already obtained t=9g-23 x 10° years for M,.= + 2°67 at 
this almost constant level, we have 
t=109 +0-32M Rol, (4) 


where M,,, must be taken as the bolometric absolute magnitude of subgiants at 
the almost constant luminosity level in the H—R diagram at which the evolutionary 
track first moves markedly to the right. 

4. 5 Ertd. and » Herc. A.—Ina survey of nearby subgiants Eggen (8) obtained 
M,=+3°67, B-V=+0-94 (T,~4650°K) for the star S5Erid. Since the 
trigonometrical parallax is very large (0-109”) the value of M,, cannot be seriously 
doubted. The bolometric correction for such a Ko IV star is probably 
about 0-37 magn. (9), yielding M,, = =~ +33. Substitution in (2) gives 
t=1'47 x 10! years. 

The same result follows for the G5 IV star »Herc.A. The appropriate 
values are My= + 3°58, B—V=+0°75 (T,=5100°K), bolometric correction 
~o0-28, M,,~ +3°3- This star has parallax 0”-108. 

A fair number of Type I subgiants are also known at comparable values of 
M,,» but the determinations of M, are somewhat less certain than they are for 
6 Erid. and pw Herc. A. 

5. The B® problem.—The question now arises as to how far the above age 
estimates must be modified to take account of recent changes in our understanding 
of the termination of the P-P chain, viz. 

He*(x, y)Be’(e~, v)Li’(p, y)He* 
L, (p,y)B*(B*, »)Be** (a)He! 
competing with He*(He’,2p)He*. If the P-P chain were always completed 
through B® all estimates of lifetimes would be reduced by the factor 0-73. The 
reason for this is simply that the effective energy output from hydrogen burning 
is then reduced by 0-73, so that the time required for a star to burn a given quantity 
of hydrogen at a specified luminosity is reduced by the same factor. 

According to an analysis by Fowler (xo), the situation is much more complex 
than this, however. The energy loss by neutrino emission from B* depends on 
the temperature and on the concentration of He*. The numerical values given 
by Fowler show that little correction is needed up to, and including, the fourth 
solution of Table 1 for the Type II stars—the central temperature and the 
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helium concentration are both too low for the correction to be appreciable, even 
for the largest reasonable cross-section for the reaction Be’(p,y)B*. Likewise 
little correction is required for central temperatures above 20m. °K, since the 
main energy generation then comes from the C—-N cycle. This leaves only the 
range of solutions where the central temperature ranges from ~14m. K to 
~ 20m. °K. 

The following prescription should give a reasonably accurate estimate for the 
corrections to be applied for the Type II stars: 

T..< 14m. deg., no correction (Y small), 

14< T,< 20m. deg., 20 per cent reduction, 

T..> 20m. deg., no correction. 


This prescription amounts to an overall reduction of age by ~o-go. Probably 
the reduction is then somewhat exaggerated, but the uncertainty in the final age 
is probably less than 5 per cent—the factor could scarcely be as low as 0-85 or 
as high as 0-95, although it might well be as high as 0-g2 or 0°93. Formula (1) 
then becomes 
t= 109 %+0-3Mp years, (5) 
The reduction factor must be greater for Type I stars, because the initial 
helium concentration is almost certainly higher in Type I stars. Also the central 
temperatures in ‘l'able 2 are systematically higher than in Table 1. Against this 
latter point, however, it must be remembered that the solutions of Table 2 refer 
to an evolutionary sequence more luminous by about o-7 magn. than the 
sequences of 6 Erid. and « Herc. A. Thus the central temperatures along the 
sequences of the latter stars would probably fall closer to Table 1 than to Table 2. 
\s a reasonable estimate for old Type I stars, a reduction factor of 0-85 may 
be used, again probably correct to within 5 per cent. Then formula (4) becomes 
t = 10°-04+0-32M Bot years, (6) 


where 1/,,, is the absolute bolometric magnitude of the old Type I subgiants 
—e.g. 6 Erid. and » Herc. A at M 


4 > 


~ ~s 
Bol = J 3° 


6. Conclusion.—If we accept the inference that the stars of the globular clusters 
are older than the oldest Type | stars (an inference that can be based either on 
dynamical arguments or on considerations of nucleogenesis) the value of t given 
by (6) must be less than the value given by (5). Immediately we have 


o-30(M, er Lyrae > 0°32(M,, )s Erid 0°86. 
Setting M,,,= +3°3 for 5 Erid. (or for » Here. A) we obtain 
(M, er Lyrae stars > +0°6, 
in agreement with the recent work of Eggen and Sandage (7). And with 
M,> +0°6 in (5) we have t>1-2 x 10! years. The age of the Galaxy is in excess 
of 10'° years. 

In view of the cosmogonical significance of these conclusions it is of interest 
to consider what errors might possibly influence the above arguments. A change 
in the calculated age of the Galaxy could come about from: 

(a) errors in the observed apparent magnitudes of 5 Erid. and pw Herc. A, 

(6) errors in the parallaxes of these stars, 


(c) errors in the bolometric corrections, 
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(d) errors in the assumed initial helium—hydrogen ratio, 
(e) an error in estimating the effect of neutrino emission from B®. 


Before considering these possibilities it may be noticed that an error in the 
opacities used in the calculations would not affect the age calculation in an 
appreciable measure. An error of opacity would change the theoretical estimate 
of the masses of the stars. In this connection it may be noted that the present 
calculations give a mass 1-093© at subgiant bolometric magnitude + 2-67. 
For other values of M,,,, the homology relation (with L,oc ~ M® on the main- 
sequence—WM is the mass of the star) gives 


login = =0'2523 —0°08M,,,. 


‘Thus for 6 Erid. and » Herc. A at Mg, = +3°3 we obtain Mzo-970. 

Likewise a moderate error in the overall energy generation coefficient for the 
P-P chain would not affect the age estimates. Such an error would change the 
radius R, slightly (cf. preceding paper) and would hence alter 7, somewhat. 
But M,,, for the subgiants is substantially independent of T,, as can be seen 
from Table 2. Hence such an error would be of no real consequence in the 
present considerations. 

Of the possibilities listed above, (a) and (b) are too improbable to be seriously 
considered. A bolometric correction of 0-28 magn. was used for uw Herc. A 
at B—V=+0°75. It seems most unlikely that this correction is significantly 
too small—which it would have to be to change the above estimate for the age of 
the Galaxy at all appreciably. Nor can the age estimate be much reduced by 
(e). We have already noticed that 0-73 gives the greatest possible factor of 
reduction in the ages, corresponding to the extreme case where the P-P chain 
always proceeds through B*. In our estimates for Type I stars the factor used 
was 0-85. If instead the factor 0-73 had been used the age inequality for the 
Galaxy would have been t > 1-07 x 10!° years. This estimate is certainly too low 
since the P-P chain does not always pass through B*. 

The initial helium—hydrogen mass ratio of the old Type I stars was taken as }. 
This falls close to the value given by Suess and Urey (11), and is close to the 
value associated with the primaeval solar material. Since the Sun is apparently 
much younger than the old Type I subgiants it would indeed be surprising if the 
solar helium—hydrogen ratio were less than the value appropriate to the latter 
stars (as it would have to be to reduce the calculated age of the Galaxy). Rather 
should the reverse be the case. 

It may be remarked that an increase in the helium—hydrogen ratio would 
decrease the calculated masses of 6 Erid. and uw Herc. A below the value 0-97 
obtained above. Accurate mass determinations would provide a test of this 
question. 

Turning now to the value of M, appropriate to the RR Lyrae stars, we ask: 
how would the above conclusion (M, > +0-6) be affected if the Type II stars 
were taken to have initial composition X =0-75, Y =0-249, instead of X =o-g9, 
Y=o-00g9 (the metal abundance remaining unchanged at A=4x10-*)? A 
qualitative answer can be given at once. Since the amount of hydrogen is thus 
reduced, the corresponding evolutionary development at a given luminosity 
level proceeds more quickly. Thus to obtain the same age for two evolutionary 
sequences, one with X =o-gg (IIa) and the other with X =0-75 (IIb), the general 
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luminosity level for IIb must be lower than for Ila. Hence the condition that the 
ages of Type II stars be greater than the ages of Type I stars demands a lower 
luminosity level in case IIb than in case Ila. Hence for case IIb the RR Lyrae 
stars must have M, fainter than for case Ia. 

To examine this effect quantitatively an evolutionary sequence for case IIb 
was computed at mass 2°3128 x 10% gm. Results are given in Table 3. Un- 
fortunately the value of M used in the calculations was too large, causing the 


TABLE 3 
Evolution of Type II star 
M=2:'3128 x 10" gm: initially X=0°75, Y=0°249, d=4x 107° 
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general luminosity level to be too high. By applying the homology condition 
the relevant results can readily be extracted, however. The luminosity level in 
Table 3 1s higher than the level in Table 1 by about 0-7 magn. Remembering 
that the level for the RR Lyrae stars fell at My = +0°5 for the sequence of Table 1, 
the level for the sequence of Table 3 must fall at M,~ —o-2 according to the 
homology assumption. The corresponding age is 4-78 x 10° years. 

\n argument exactly similar to that leading to (1) now gives the following 
general relation between the age ¢t and the luminosity level M, of the RR Lyrae 
Stars: 

t= 10°°74+0-3M 7 vears. (7) 


Once again a correction must be made for neutrino emission by B*. In the 
present case the correction factor must be closely the same as for the Type I stars, 
since the helium—hydrogen ratio is closely the same. Thus we now apply the 
correction factor 0-85 instead of the factor o-go used formerly for case Ila. This 
gives 

t = 10°67 +0-3My vears. (8) 
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The requirement that the age (8) for the case IIb be greater than the age (6) 
of old Type I stars (with M,,= +3°3, as for » Herc. A and 6 Erid.) leads 
immediately to 
My, > +1°3. 

This is the condition imposed on the visual magnitude of the RR Lyrae stars. 
Such a low value of M, makes it appear unlikely that Type II stars correspond 
to case IIb, the upper limit +1-3 being significantly fainter than the value 
suggested by observation (7). 
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A NOTE ON THE RRLYRAE VARIABLES 
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Summary 

The pulsation criterion P+/p=Q is used to derive the absolute magni- 
tudes of the RR Lyrae variables as a function of their colour and period; an 
error in a similar derivation by Sandage is corrected. According to this 
model the observed difference in mean periods of the variables in M3 and 
w Cen is due to a difference of 0-3 absolute visual magnitude between the stars 
in the two clusters. 

Rough agreement is found between the predicted absolute magnitudes 
and those derived from trigonometrical parallaxes for four field variables with 
P<o:2 days. 

The pulsation model predicts a difference of o™-7 between the abnormal 
short period ab types found by Gaposchkin at the galactic centre and the 
normal ab and ¢ types in the same field. Since this is not observed, either 
the model or the periods of the abnormal variables are incorrect. 

A list is given of 16 stars of ab type and P <0-36 days in the general field. 
These stars are shown to lie closer to the galactic plane than a sample of 
normal ¢ type variables of the same period, in agreement with Detre’s 
observation of a relatively smaller solar motion for these stars. Further 
observations are required for such stars both at the galactic centre and in the 
general field. 





In 1940, Schwarzschild (25) showed that the RR Lyrae stars in M3 were 
restricted to a definite zone in the colour magnitude diagram and that the 
pulsations of those variables of Bailey type ab obeyed the period—density 
pulsation criterion Py p=QO*. Recently, Roberts and Sandage (21) determined 
colours and magnitudes for 47 RR Lyrae stars in M3 on a photoelectric scale. 
They showed that both Bailey types ab and ¢ obeyed the pulsation criterion with 
O,,/Q. ~1°5, so that type ab probably vibrate in the fundamental mode and 
type c in the first overtone. Sandage (24) has further suggested that the 
pulsation model may explain the distribution of mean period (P) of these variables 
in globular clusters which was found by Oosterhoff (18). For one group of 
clusters (e.g. w Cen, M15, M53 and Mgz2) P,, ~o-63 and P, ~o-37, while for 
the other group (e.g. M3, M4 and Ms) P 


» ~O'55 and P.~o-31. Sandage’s 
calculations are repeated here because, in obtaining his relation between absolute 
magnitude and period (viz. AM,=—2-8AlogP), he used an incorrect sign 
in the expression for the mass-luminosity relation. His corresponding 
calculation for galactic cepheids (23) is however correct. 

* This relation holds for homologous stellar models. It was first found by A. Ritter (Wied. Ann., 
8, 172, 1879) for a homogeneous star, and H. Shapley (Ap. 7., 40, 448, 1914) first applied it quanti- 
tatively to the cepheid variables. Later, Eddington (M.N., 79, 2, 177, 1918) found that the relation 
also holds approximately for stars of a non-uniform density distribution. 
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The pulsation criterion may be written as 


fal rm z 
log P+ 4 log —— +0°3(Mpoi—J l £ =| ‘ 
og 4 8 AIG 0°3(Mpo1— Moi ©) + 3 log TO og O (1) 
From Morgan and Keenan (17) and Johnson and Morgan (9), we have: 


log T, = 3-983 — (0°332 +0°030)((B—V) for o17<(B-—V)<o-40 (2) 
and from Kuiper (13) and Eggen (2), we have: 
Mooi = My — 0°23 + (0°28 +0°16)((B—V) for o<(B—V)<o-7 (3) 
while assuming that the slope of the main-sequence mass-luminosity relation 
(Kuiper (14)) may be applied as in (23), we have: 


( 


fH 
log ~ =const. — (0-122 + 0°014) Myo. 
og MoO ons ( 4)Mpvo1 (4) 


Substituting (2), (3) and (4) in (1), we have: 
My, =(+3°89 + 0°41)(B— V)— (4°18 + 0-12) log P+ (4°18 + 0-12) log O + const. (5) 


The quoted errors are merely the uncertainties in obtaining linear relations from 
the data; the actual errors of the coefficients in equation (5) may be somewhat 
larger because the relations (2), (3) and (4) do not apply exactly to RR Lyrae 
stars. ‘This error of application is likely to be quite small for the bolometric 
correction in equation (3). Iwanowska (8) has deduced from a mass—luminosity 
relation of Parenago and Masevitch (19) that the coefficient of Mpoi in (4) varies 
from —o-100 for pulsation stars of Population I to — 0-086 for all Population II 
stars. ‘These values lead to an error in the coefficient of log P and logQ in 
equation (5) of +0-32. The colour-temperature relation in equation (2) is 
derived from main-sequence stars. According to Sandage (22), however, 
RR Lyrae stars have a surface gravity 0°03 that of main-sequence stars of the 
same 7',, and the electron pressure is reduced by a factor 5. The application 
of the slope of the plot of (B— V) against T, for main-sequence stars to RR Lyrae 
stars (see Fig. 2 of (22)) will therefore introduce errors which will show up as 
an error in the coefficient of (B—V) in equation (5) which it is difficult to 
estimate. The uncertainty in the coefficient of (B—V) is therefore likely to be 
somewhat larger than that quoted above. 

The zone of instability in M3 is 0-17<(B—V)<0-23 for c type variables 
and 0-23<(B—V)<o-39 for type ab. Sandage (24) has assumed that this is 
true for all globular clusters, i.e. the zone is vertical on a c-m diagram. There 
is no observational evidence for this because of the unknown amounts of reddening 
in different clusters. For W Virginis stars and classical cepheids it is known 
that the zone has a slope of A(My)/A(B—V) ~11. The adoption here of a vertical 
zone does not materially affect the arguments in this paper. 

Fig. 1 shows a colour—-magnitude diagram for RR Lyrae stars in the form 
given by Sandage (24). Lines of constant period are drawn calculated from 
equation (5) assuming Q,,/O,=1°5 (Roberts and Sandage, Joc. cit.). The M3 
variables with known photoelectric colours and periods (further than 180” 
from the cluster centre) are plotted in the figure, using their colours and periods 
alone. Plotted in this way, a range of 0-4 in absolute visual magnitude is predicted ; 
it must be noted that the ordinate in the figure (My) has an arbitrary zero. The 
predicted range must be diminished when allowance is made for errors in colour ; 
it could in fact be produced entirely by a maximum error of + 0-05 in the observed 
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colour. The directly observed spread in visual magnitude for the selected 
ab stars in M3 is 0-23 magnitude. 

The ratio P,,/P, of the periods of ab type and c type variables with the same 
My, in any cluster is found from equation (5) to be: 


log? =0'93[(B—V),, —(B—V),]+log ee ‘ (6) 
It would be desirable to see if Q,,/O, is the same for the w Cen group of 
clusters (P,,=04-63 and P,=0137) as for the M3 group (P,,=0%55 and 
P,=0%-31), but sufficient variables of known colour are not available for the 


Bailey type c Bailey type ab 
Y ee eee SY 


(arbitrary zero) a 
aoe (A) TEN 


M3 


Galactic Centre 
(Groups II ond I) 





Galactic Centre 
(Group) 


a | ome A 


————-- 


/ 








j i. . rr we 
+ 0-20 + 0.30 (B- Vv) +0.40 





Fic. 1.—The colour-magnitude diagram for RR-Lyrae variables in the form given by Sandage 
(24). The ordinate My is computed from the colour and period by equation (5) with Qgp/Qe= 1°5 
and is adjusted so that the variables in M3 have zero mean magnitude. 
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w Cen group. If we assume that the mean colours of the ab type and the c type 
variables are the same for all clusters and further if Q,,/Q, is constant, we 
should expect P.,/P, to be constant. The values of P_,/P, for the w Cen group 
and the M3 group are 0-170 and 0-177 respectively, suggesting that Q,,/Q, is 
in fact constant. In this case, equation (5) predicts that the variables of the 
w Cen group are ~o™-3 brighter than those in M3. Also, it can be seen from 
the figure that the w Cen variables should be bluer than those of the same period 
in M3, when absorption effects have been removed ; it is therefore very necessary 
to obtain colours for the variables in a cluster of the w Cen group to check this 
prediction. At present we can only conclude that the pulsation model is consistent 
with the available observations of the RR Lyrae variables in globular clusters. 

Additional checks on the model may be obtained by considering its implications 
for the short period variables in the general field and this is discussed in the 
following sections. 

Variables with 0-05 <P <o-2 days.—About 30 of these stars are known and 
they are often distinguished from RR Lyrae stars for which 0-2 < P<1-o days. 
Absolute magnitudes from trigonometric parallaxes are known for four of them 
and Woltjer (27) has attempted to estimate absolute magnitudes for SX Phe 
and AI Vel from radii determined from phase shifts in the light curves. The 
absolute photographic magnitude (Mpg) has been calculated for these stars by 
equation (5) using Eggen’s colours (3) for SX Phe, p Pup and 4 Sct, while the 
colour of AI Vel has been deduced from the mean spectral type. Mpg was assumed 
to be zero for the variables in M3 in order to determine the constant in equation (5). 
The values found are compared with those found from parallaxes and by Woltjer 
in Table I. It may be noted that the values derived from parallaxes and from 
the pulsation model will be similarly affected by a small unknown amount of 
reddening e: the parallax magnitude will be too bright by 4 and the pulsation 
model magnitude by 4-9. The assigned Bailey type can be no more than a guess, 
in view of the variable nature of the light curves of these very short period stars ; 
this considerably strains the applicability of a simple model with one pulsation 
period. Considering this and the likely errors in the parallaxes, the discrepancies 
found for individual stars might have been expected. Both model and trigono- 
metrical parallaxes however give a mean Mpg for the group of +2°5 which 
suggests that the coefficient of log P in equation (5) is roughly correct. 


TABLE I 
The absolute photographic magnitude (My,) for stars with o-o5 <P <0-20 days 
Assumed Mpg 
Star Bailey Period B-V_  Pulsation Trigonometrical Woltjer 
type (days) model parallaxes (27) 

SX Phe ab 0056 +°20 +3°7 +4°'1t0'°9 +4°5 
(3) (27) 

Al Vel ab O'1It +°26 +2°7 +4°3t+0°4 +2°2 
(27) (26) 

p Pup ; O'14! +42 +2°4 —0°4; +03; +1°6 
(3) (3) 

8 Sct "194 +°34 +1°4 +1°7; +16 
(3) (3) 


RR Lyrae variables at the galactic centre.—It is found that in globular clusters 
the variables with periods ~o-3 days always have symmetrical light curves 
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(Bailey type c). Greenstein (6) noticed, however, that there are a few stars 
of type ab of this period in the general field. Recently Gaposchkin (4), in studying 
a field near the galactic centre, has found that about one-half of the stars of this 
period are of type ab. Pavlovskaya (15), however, claims that some of 
Gaposchkin’s short periods should be o-5 instead of 0-3 days. 

A comparison of Gaposchkin’s stars with a sample from the general field 
(on a plot of the period (P) against the asymmetry (/—m)) shows that the 
greatest difference for the two samples is for stars with 0-20 < P< 0-36 days and 
(M—m)<0-22; ~30 per cent of Gaposchkin’s stars lie in this region compared 
with ~1 per cent in the general field. Ninety-three of Gaposchkin’s stars 
taken from the General Catalogue (5) are divided into three groups according 
to period and asymmetry in Table II. 


Tas.e II 
Stars at the galactic centre (Gaposchkin) 
: Mean 
Bailey Period Mean Mean 


M—m apparent pg No. 


GC 
3roup type range period (M—m) magnitude 


| ab 0'20<P<0°36 0°32 <0'22 O17 17°62 +0°08 36 
II ‘ 020 < P<0°36 0:28 > 0°22 0°34 17°57 +0°06 32 
III ab P> 0-36 0°43 o'1g 17°63 0°15 25 
Asymmetry (/—m) is the time of brightness increase for the variable expressed as a 
fraction of the period. 


Each group of Gaposchkin’s stars shows the same mean apparent magnitude 
with only a small scatter; it may be inferred therefore, following Baade, that 
they all occupy the same volume in space and have the same absolute photographic 
magnitude within o-1 magnitude. If we assume that the mean colour of the 
ab variables is about 0-1 magnitude redder than for the type c, then the My for 
all groups is practically the same. ‘The horizontal dashed lines in the figure 
labelled Galactic Centre (Groups II and III) and Galactic Centre (Group 1) 
indicate the approximate absolute magnitude (/,) which the pulsation model 
predicts for these groups as judged from their mean periods and Bailey types alone 
since no colours are available. It may be seen that the stars of Groups II and III 
have the same absolute visual magnitude, which is ~o™-25 fainter than for the 
variables in M3, and which agrees with the observation that the mean apparent 
magnitude of the two groups is the same. ‘The asymmetrical short period stars 
of Group I are however anomalous on the pulsation model, since this predicts 
that they should be about 0-7 magnitudes fainter than the other two groups and 
this is contrary to Gaposchkin’s observations. In the light of this discrepancy 
and also the results of Pavlovskaya’s re-analysis, it seems very desirable that 
new observations should be made to check the periods, colours and light-curves 
of these galactic centre stars. 

RRLyrae stars with o20< P<0-36 days.—As Greenstein noted (loc. cit.) 
some RR Lyrae stars of short period ( ~o-3 days) and type ab do occur in the 
general field. Baade (1) has suggested the possibility that RR Lyrae stars of this 
type may be present in the Galaxy with a disk-like distribution. A list of 69 stars 
was therefore made from the General Catalogue (5) of all RR Lyrae stars with 
0:20 < P<0-36 days and with known asymmetries (M—m). Sixteen of these for 
which (/— m) < 0-22 may be considered to be of type ab and are listed in Table III. 
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Tasie III 
Stars in the general field with 0-20 < P<0°36 and (M—m) < 022 
Star M-—m P 1 b m, orm, R.A. (1900) Dec, 


v672 Aql 0°20 0°346 14° —10° 12°6-13°5p 1oh4oM™4q8 +07°46"9 
v753 Aql 007 «0'299 «18 12s g8-15°5P =—19 55 21 = F959 


IU Car o'18 , 10°4-I1'5p 06 51 40 —59 281 
QR Cas 0°20 ° 14°O-15°4p 23 46 12 +55 08 


SW Cru 018 328 1I'Q-12°7p II 54 03 —59 44'8 
IV Cyg O'10 , 13°6-14'9p 19 46 28 + 34 04°2 


RV Del , , 13 12°9-13°8v 20 16 04 +14 56°3 
BK Del , , 17 13'0-14'6p 20 24 38 +11 33°7 


CP Del , "345 2 14 14°4-15°6p 20 21 18 +14 28 
DQ Del 17 , 28 17 13°9-15'3P 20 36 30 + 14 23°4 


XX Hya “I 337 +22 10°5-12p 09 05 08 —15 118 
TV Lib ° , +38 10°8-12°3p 15 13 00 —08 05°8 


AQ Lyr ‘16 , 23 +14 12°0-13°3V 18 30 48 +26 28-8 
DD Lyr org ; 25 +10 12°7-13°8v 18 49 07 +26 56°4 


v413 Oph 0°16 2 332 +25 11°8-13°6p 16 19 42 —10 18°4 
v1176 Ser 0°16 "35 338 — 6 12°0-12'9p 18 18 24 —21 38 
For some of these stars the periods are doubtful on the basis of large residuals 
between the observations and the quoted light curves; on the other hand, SW Cru 
for example has a period which has been well determined over a number of years. 
A plot of asymmetry against galactic latitude (b) for all 69 stars appears to 
show that the asymmetrical stars do occur at lower latitudes than the remainder. 
If the sky is divided into two zones of equal solid angle—an equatorial zone 
(|5]< 30°) and a polar zone (|b|> 30°), then the distribution into these zones 
according to (M-—~m) is as given in Table IV. Application of Fisher’s exact 


TaBLe IV 
Distribution in galactic latitude of short period variables according to asymmetry (M—m) 
(M—m) < 0°22 (M—m) > 0°22 
Equatorial zone (|>| < 30°) 15 35 
Polar zone (|o|> 30°) I 18 


method for a 2x2 table gives these two populations as significantly different 
at the 5 per cent level. The mean perpendicular distance (|Z|) from the galactic 
plane was found for each group assuming zero absolute magnitude for all stars. 
The interstellar absorption for a group of 270 B stars at a mean distance of 2 kpc 
has a mean value of 1-9 mpg and shows no correlation with distance (cf. Hiltner 
(7)). This value was therefore adopted for the 13 stars which have |b|<8° 
and the cosecant law of absorption, Apg=o-28cosecb (12), was used for the 
remainder. The following values were found for the mean height from the 
plane (|Z]), the r.m.s. dispersion in Z about the plane (a7) and the mean distance 
from the Sun projected on to the plane (R) in kiloparsecs: 

(M—m) <0-22 |\Z|\=0°65, o,=0-8+0-2, R=2-7, 

(M—m)>o0-22 |Z2|=1-49, o7=2°6+03, R=3-2. 
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The mean period of the 16 stars with (M—m)<o-22 is 0-326 days. On the 
pulsation model, an absolute magnitude of +0-g relative to the M3 stars would 
be predicted. On this hypothesis and assuming the M3 stars to have zero absolute 
magnitude, the group of stars with (M—m) <o-22 will have: 


|\Z|\=0°43, o,=o6+01, R=18. 


In either case the asymmetrical group has a more flattened distribution with 
respect to the plane than the normal c type variables. The approximate equality 
of (R) when the same absolute magnitude is assumed for each group suggests 
that observational selection has not produced this result—as it might if the 
asymmetrical type stars were much nearer to the Sun than the others. 

Joy (10, 11) has given radial velocities for two of the stars with (M—m) <o0-22 
(XX Hya and TV Lib); both have the same low velocity (— 10km/s) which is 
in accord with a disk-like distribution. Pavlovskaya (20) reports that Detre* 
has investigated RR Lyrae stars with 0-20 <P<o-40 days of Bailey type a and 
has found a solar motion of 48km/s which is appreciably less than for c type 
variables of the same period. ‘This is to be expected if the asymmetrical group 
(type a) has a more flattened distribution than type c stars of the same period. 

[wanowska (8) has separated RR Lyrae stars into two sequences according 
to her observed mean spectral types and period, and shown that stars with 
P<o:s5 days form a sequence of later spectral type than those with P>o-5 days. 
Repeating this with additional stars for which a reliable spectral type is available, 
the effect, although small, is confirmed for type ab stars, as is shown in Table V. 


TABLE V 
Mean spectral type and mean period for type ab variables 
Period range Mean period Mean spectral type 
020<P<o'50 0°442 A7‘4 
P>0'50 0°589 A5‘0 

On the pulsation model the group of larger mean period should be nearly 
o™-5 brighter in absolute visual magnitude. Lubienska-Iwaniszewska (16) has. 
shown that the group of shorter mean period has a smaller height above the 
plane and lower mean velocity than the other group. This suggests that the 
properties which we have found for the group of short period asymmetric stars 
in the present paper are an extreme example of a general tendency for the ab 
variables in which their space distribution, space motion and spectral type are 
a function of their period; presumably a similar relationship holds for the type c 
variables. In this case the sixteen stars listed in the present paper would be 
expected to have a later mean spectral type than type ab variables of longer period 
or type c variables of the same period. 

There is no indication for the stars with 0-20 < P< 0-36 days that the type ab 
variables exist preferentially in the galactic centre hemisphere compared with 
those of type c. A preponderance of type abd stars of short period would not 
therefore be expected at the galactic centre on this meagre evidence. The 
well-known increase of RR Lyrae stars of all types towards the galactic centre is 
however shown in Table VI. 

To summarize, the pulsation model has some success in relating the different 
mean periods of the type ab and type c variables found in globular clusters and in 


* Paper unavailable to author. 
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Tasie VI 


Concentration to galactic centre as function of asymmetry (M—m) 
(M—m) <0-22 (M—m)> 0°22 


Galactic centre hemisphere 13 42 
Anti-centre hemisphere 3 9 


predicting the absolute magnitudes of very short period variables within the 
observational errors and likely theoretical limitations. The disagreement with 
Gaposchkin’s observations together with the results of Pavlovskaya’s re-analysis 
make further observations of RR Lyrae stars at the galactic centre desirable. 
The short period variables of type ab in the general field lie closer to the plane 
than the type c variables of the same period and have a smaller solar motion. 

It would seem worth while to check on the periods of the 16 stars listed in 
Table III and also discover if they differ from the type c variables of the same 
period or type ab variables of longer period in physical properties (such as colour 
variation through the cycle or spectrum). To this end, spectra are now being 
secured at the Radcliffe Observatory of SW Cru and IU Car; the remaining stars 
are within reach of Northern Observatories. 


I should like to thank the Radcliffe Observer (Dr A. D. Thackeray) for 
suggesting some of these problems to me and for his advice in writing this paper. 
I am also grateful to Drs A. R. Sandage and M. W. Feast for several most 
helpful discussions. 

This work was carried out partly with the aid of a Radcliffe Travelling 
Fellowship and partly with the aid of a grant from the Department of Scientific 


and Industrial Research to the Radcliffe Trustees. 


Radcliffe Observatory, 
Pretoria: 


1958 December. 
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THE SYSTEM OF ¢ PHOENICIS 
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Summary 


Radial velocities have been measured from 71 new spectrograms of 
¢ Phoenicis, taken at Mt Stromlo. 

The spectroscopic orbit of the primary star has been computed, and the 
semi-amplitude, K,, of the secondary star has been determined from differ- 
ential measures. Absolute dimensions of the two components have been 
derived by combining the spectroscopic figures with the photometric results 
of A. R. Hogg, and the spectral type of the secondary star has been computed. 
Certain difficulties of the system seem to be solved. 





1. Introduction.—The system of { Phoenicis («=01"06™-3, 5= —55°31° 
(1950), HD 6882) belongs to the class of double-line spectroscopic binaries that 
are eclipsing variables as well. It was reported as a spectroscopic binary by 
Wilson from four plates obtained at Santiago (1). Altogether 16 plates were 
secured at Santiago, from which Colacevich (2) determined provisional orbital 
elements. The lines of the fainter component were found on but five plates and 
then only with difficulty. The system was reported as an eclipsing variable by 


Hogg from photoelectric measurements taken at Mt Stromlo between the years 
1947-49. His results and certain discrepancies between the photometric and 
spectroscopic solutions were discussed in this journal (3). Spectroscopic 
observations with the Cassegrain spectrograph were made at Mt Stromlo during 
the years 1954-57 at the suggestion of W. Buscombe. 

2. Observations, spectral lines, measurements.—Altogether 80 plates of the 
spectrum of ¢ Phoenicis were taken at Mt Stromlo. During the years 1954-56 
mostly Kodak 103a-O were used, but during the observing period of 1957 
Kodak IIa—O were used at the author’s suggestion, as the finer grain of these 
plates enhanced the very faint lines of the weaker component. The linear 
dispersion was about 35 A/mm at Hy. ‘The slit-width used from 1954-56 was 
0-04 mm, and was changed later to 0-06 mm. This corresponded to a projected 
slit-width of o-o1 resp. o°o15 mm. The spectra show the following lines between 
A 3922 to A4528, the range used for radial velocity measurements: He 4026, 
4388, 4471; H 3970(He), 4102 (H8), 4340 (Hy); Mg 4481. From the 
available 80 plates, 71 were measured for radial velocities by the author; 9 plates 
were rejected as too weak. The quality of the lines for radial velocity determina- 
tions was poor. Lines due to He and Mg were faint and diffuse, and not visible 
on some plates because of variations of the density of the plates. The H-lines, 
although strong, were broad; H 3970, moreover, showed as a heavy blend on all 
spectra, and was unusable. Accordingly the lines H 4102 and H 4340 were 
mostly the only ones to be used for the determination of radial velocities. The 
measurements were made with a Zeiss-Abbe Comparator, the advantages of 
which will be described later. 
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TAaBLe I 
Radial velocity observations of © Phoenicis 


Date (U.T.) Phase Obs. Velocity1 Velocity 2 
d 
1954 Sept. 6-69 26 WH + 38 
12°67 56 WB — 65 
14°57 "79 WH — 106 
58 ‘80 WH 104 
57 "09 WB 116 
S*51 +36 AP 64 
‘61 40 AP + 85 
ef 56 AP 126 
"50 68 WB 103 
60 78 WB 
‘68 “86 WB 
‘Ol HG 
‘10 HG 
‘19 HG 
‘00 WH 
‘II WH 
"34 AP 
AP 
“42 AP 
"52 AP 
‘34 WB 
44 WB 
“51 WB 
“a3 AP 
AP 
HG 
HG 
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TaBLe I (continued) 


Date (U.T.) Phase Obs. 
d 

1957 Oct. 1°46 
. 1°56 

04 

“II 

"22 

1‘20 

1°25 


Fe PS a Pe a Pe ee By 


Observers : W. Buscombe (WB), H. Gollnow (HG), W. D. Heintz (WH), 
A. Przybylski (AP). 
3. Period.—Colacevich determined a period P=1%-66958 from the early 
spectroscopic observations; Hogg from the photometric observations found 
P=14-66990. A new period P=14-6697597 was derived, using the old and 





| "ed 


T=J.D. 2435660:500 








| 1 i j i 
° : ‘ 2 ; t2 a 6 8 


PHASE [days] 
Fic. 1.—All individual radial velocities of the primary star of ¢ Phoenicis plotted against the 
radial-velocity curve drawn according to the final elements. Lick observations are indicated by 


crosses. 





new spectroscopic observations. ‘This period represented all observations, 
spectroscopic and photometric, with the exception of the first observation listed 
by Colacevich in his Table I (2, p. 85). According to a private communication 
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from G. H. Herbig, the published J.D. should read ... 39-93 instead of ... 31°93. 
This correction removed the discrepancy. Between the first and last observation 
the system has completed 10 800 revolutions. Fig. 1 shows the radial velocity 
curve with all observations for the main star. The epoch of maximum radial 
velocity is given by J.D. 24356600 + 19-6697597E. 
Tase II 
Normal places 
Observation 


N Mean phase Mean R.V. Residuals 
INOS. 


12, 57, 27 "04 + 133 
44,45, 5 "08 +130 
53, $2, 33 *I2 +115 
31, 14, 59 
32, 34, 4! 
7°, 17, 42 
18, 7 

35; 
46, 
64, 
66, 
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Fic. 2.—The 24 normal points from Mt Stromlo observations of the primary star of 
€ Phoenicis shown against the radial velocity curve from the final elements. 











No. 2, 1959 The system of € Phoenicis 


4. Orbital elements.—Orbital elements for the main star were derived from 
Mt Stromlo observations only, as this material is homogeneous and the Lick 
observations are not numerous. In such a case Schlesinger suggests the 
omission of early observations, because the residuals of these observations depend 
largely upon the adopted period (4). The period was adopted as final and 
omitted from the solution. Preliminary elements found graphically confirmed 
Hogg’s result that the eccentricity was very small, and hence in a first solution 
a circular orbit was assumed and corrections to the elements deduced by Sterne’s 
method (5). The observations were formed into 24 normal places given in 
Table II and Fig. 2. A second, and final, solution was made for an elliptic orbit, 
using also Sterne’s method. The results are shown in Table III. 


Tase III 


Orbital elements of £ Phoenicis 


Element Circular d Elliptical (adopted) 
1°6697597 
"44 +12 


033 "034 + ‘O15 
94° 97°+25° 
121°5 km/s 121°4+1°9 
15°2 km/s 1s‘4t1°1 
Note.—The time T is the date of periastron passage, as usual. Sterne uses the “ epoch of 
the mean longitude ”’ instead, but the dates can easily be converted into each other. 
The errors are all m.e. 


5. Absolute dimensions of the components.—The orbits and absolute dimensions 


of the components can be derived by combining the results of the spectroscopic 
and photometric observations, if the lines of both components are measurable 
for radial velocities. Besides the elements for the primary star, it is sufficient to 
know the semi-amplitude, K,, of the secondary star. In the case of { Phe the 
lines of the secondary star were so weak and diffuse, and visible only near the 
nodal passages, that it was not thought advisable to try to determine very uncertain 
orbital elements. The value K, was found from differential measures of primary 
and secondary spectra on a number of plates obtained near the nodal passages. 
Distances between the lines were measured in microns, which, converted into 
velocities and combined with the results for the primary star, gave the velocities 
for the secondary star. The results are also shown in Table I. The Zeiss-Abbe 
comparator was of particular advantage for measuring the distances between the 
lines. The comparator has no screw, and so the plate carriage can be moved 
rapidly in both directions under the adjusting microscope. During this movement 
the eye can distinguish the very faint line of the secondary star more easily and 
gets a more definite impression of it. When the impression is sufficiently definite, 
the carriage is stopped with the index mark on the line, and the reading is made. 
From all available plates the value of K, was found to be 247+ 10 km/sec. The 
combination of the spectroscopic results with the photometric figures of Hogg (3) 
gave the quantities shown in Table IV. 

6. Spectral class, absolute magnitude, parallax.—From Mt Stromlo plates 
A. de Vaucouleurs has classified the primary star on the MK system as B6 V (6). 
Due to the poor quality of the lines, it was not possible to classify the secondary 
star and it was thus computed according to the method described by L. Plaut (7).. 





G. Hagemann 


TaBLe IV 


Absolute values from spectroscopic orbit and photometric results 
Star 1 Star 2 
Semi-major axis of the orbit km x 10 2‘80 5°70 
Mass Solar Units 6-1 ‘Oo 
Radius 3°4 zo) 
Mean density 17 37 
Using Hogg’s value of the ratio of the surface brightnesses, J,:J,=1-8, the 
computed spectral class of the secondary star becomes AoV. ‘The absolute 
magnitude as given by Morgan and Keenan (8), is —1-1 for the primary, and 
+ 0°3 forthe secondary star. The spectroscopic parallax, using A. de Vaucouleurs’ 
modulus m,—M=5-51, is 0”-008. The difference of 1™-4 between the two 
components is in good agreement with Hogg’s value of 1™-5, from the ratio 
L,/L-. 

Hogg lists a number of discrepancies between the values derived from 
photometric and spectroscopic observations. ‘These difficulties are based on 
the assumption of the same spectral class for both components as B8 + B8, as 
given in the H.D. Catalogue, and on the values for the spectroscopic orbit of 
the secondary star by Colacevich. With regard to the spectral class of the 
secondary star of spectroscopic binaries, Plaut remarks generally: ‘‘ The 
computed spectral type of the secondary component is to be preferred to the 
observed one’’. Already Hogg had come to the conclusion that the spectral 
type of the secondary is definitely later than B8. ‘The orbital elements for the 
secondary star by Colacevich, however, must be considered as very uncertain. 
His radial velocity curve for the secondary star rests on five observations only, 
and these are not favourably placed for the only important element, Ks. 

Assuming now the primary and secondary spectra to be of types B6 V and 
\o V respectively, the present observations lead to results for masses, radii and 
densities which, as shown in ‘Table V, agree well with the normal figures for 
stars of these types given by Allen (9). 


TABLE V 


Absolute values from spectroscopic and photometric observations compared with values 


from Allen’s tables for sturs of the same spectral type. (Units are solar units) 


Observed Tabulated 
Mass : primary star 6°1 I 
secondary star 3°0 ‘O 
Mean radius : primary star 
secondary star ‘O 
Mean density : primary star 16 13 
secondary star 38 28 


5 
3 
3°4 3°4 


2 


7. Conclusion.—The system of ¢ Phe seems to be a normal one under the 
assumption that the spectral type of the secondary star is Ao V. As the quality 
of the spectral lines is poor for the determination of radial velocities, it does 
not seem worthwhile to accumulate more spectroscopic observations. For the 
same reason, non-orbital effects in the observed radial velocities were not allowed 
for. Besides, the non-orbital effects influence the values of e and w, which have 
a doubtful accuracy in nearly circular orbits, as that of ¢ Phoenicis. 
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NOTES ON FIVE SOUTHERN GASEOUS NEBULAE 
David S. Evans 
(Communicated by H.M. Astronomer at the Cape) 
(Received 1958 November 27) 
Summary 
The results of observations of the spectra, radial velocities, internal 
motions and form of five southern objects, NGC 3918, 6302, 6326, IC 4776, 


and Anon 165 13™-5, —51° 52’, made with the Radcliffe 74-inch reflector, 
are reported. 





These notes aim at providing details of structure, radial velocity, internal 
motion and spectrum for these five nebulae, which have been observed in the 
course of a more general programme on planetary nebulae. The second of 
them, NGC 6302, does not, strictly, rank as a planetary nebula, although it 
has been so classified in the past, and, for the last-named, the velocity is not 
given. 

Planetary nebulae show a very great range in form and surface brightness, 
and a considerable amount of experiment has had to be devoted to discovering 
the best methods of investigation. This is partly due to the fact that the best 
methods for obtaining different kinds of data are, to an extent, in conflict with 
each other. 

For the determination of radial velocities and the description of the spectra, 
there are available four different dispersions at the Cassegrain spectrograph 
(known as a, b, c, and d) (1), together with the small spectrograph for use 
at the Newtonian focus (2). The c-dispersion of the former instrument 
(49 Aymm) has proved by far the most useful, and is now employed almost 
exclusively for Cassegrain work. ‘This camera lens has the enormous advantage 
of producing a spectrum in excellent focus to beyond 5000 A, while as far 
as 6700 A the quality remains good, and perfectly adequate for line identification 
and the determination of line strengths. 

For the reduction of these plates new tables were computed using two 
separate Hartmann formulae overlapping at 4500 A. The wave-lengths of the 
forbidden nebular lines were taken, as far as possible, from Ira S. Bowen (3), 
of permitted and other forbidden lines from Charlotte E. Moore (4), and other 
wave-lengths, e.g. of unidentified lines and of blends, from A. B. Wyse (5). 
In deriving velocities, the lines used were selected so as to exclude faint lines, 
blends, lines of uncertain wave-length, and all lines to the red of N,. An iron 
arc comparison spectrum was used, with the exposure tripled as compared 
with the normal for stellar radial velocity work, in order to give a satisfactory 
number of comparison lines. The results from the Cassegrain spectra (Table I), 
show good agreement with those of Campbell and Moore (6). It now seems 
clear that an adequate coverage for determination of velocity and spectral 
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characteristics, for those objects observable with the Cassegrain instrument, 
can be secured with three plates of each: one fairly short exposure (e.g. 
1-2 hours), on fast blue-sensitive plates (Kodak 103aQ), and two long exposures, 
each pushed as far as is practicable, one blue (103aO), and one panchromatic 
(103aF). ‘These long exposures are limited by the fact that, so far, it has not 
proved practicable to extend an exposure over more than one night, and, in 
consequence, the “depth” reached is usually not as great as that attained 
by Wyse (5). Only if an examination of the lines in the violet reveals structure 
is it worthwhile to attempt to secure spectra on higher dispersion for a study of 
internal motion. 


TABLE I 


Nebula Maximum exposures (minutes) for various 
dispersions and emulsions (103a series) 


a:0 c3O c:3B 2:7? -438 Neo Vc+m Total 
plates 
NGC 3918 99 or 260 vem 135 —16%4 
6302 < 63 320 240 oe 150 — 36 
6326 on 270 120 deb 180 155 +11 
IC 4776 ae 60 re 180 aie 120 +188 
Anon 165 13™-5 i... cas AK. sia i 180 coe 


— 51° 52° 
* Hg sensitized. 


Direct photographs at the Newtonian focus have usually been available 
for the study of structure. At the same focus the small spectrograph giving 
a dispersion of 400 A mm has been used to obtain ‘‘ deep” spectra of bright 
objects, or to reach very faint ones. It might be hoped that, by duplicating 
observations of the same object at Cassegrain and Newtonian foci, sufficient 
information might be gained concerning the systematic and accidental errors 
of Newtonian radial velocities to open up the whole field of the velocities of 
very faint southern planetaries. Unfortunately, further work will be necessary 
before this hope can be realized. For objects with large velocity dispersion, 
such as extragalactic nebulae or globular clusters, a probable error of 20 km/sec 
from a single spectrum could be tolerated. For useful results on the planetaries, 
however, this figure must be reduced by a factor of at least four. Other things 
being equal, the reduced dispersion alone may be expected to put up the 
probable error of a Newtonian result by eight times as compared with the 
Cassegrain c-dispersion. The lines given by the Newtonian spectrograph 
are curved: for all extragalactic systems yet observed, the spectrum fills the 
whole dekker, the lines may be measured at mid-dekker, and a fairly well- 
established curvature correction applied. For the planetaries, the difficulties 
of offset guiding are such that the spectrum is not always centrally placed, 
and this introduces another uncertainty. Again, the best lines for extragalactics 
are H and K, where the prismatic dispersion is twice as great as near N, and 
N,, usually the best lines for planetaries. Finally, the low dispersion of the 
Newtonian spectrograph leads to blending of lines, and the relative strengths 
of the components of a given blend usually vary from planetary nebula to 
planetary nebula. 
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However, even though satisfactory velocity determinations of the velocities 
of planetaries cannot yet be obtained with the Newtonian spectrograph, much 
information can be derived from these spectra. The instrument is extremely 
fast, and permits observations of objects completely out of reach of the 
Cassegrain spectrograph. For the brighter systems the continuum is usually 
well shown, and, although there is much absorption below 3900 A, this is 
still much less than in the Cassegrain spectrograph, so that a far better account 
of the violet and near ultra-violet spectrum is rendered. 

After these preliminary remarks we may now return to a detailed discussion 
of the objects described in these notes. In ‘Table II are listed the lines found 
in the various spectra, these being associated with numbers empirically indicating 
the relative line-strengths within each spectrum. No claim is made that these 
represent true line-strengths: it is only hoped that, in general, a larger number 
corresponds to a stronger line. The figures are based on eye-estimates (these 
have been found to be surprisingly consistent), adjusted for plate-sensitivity 
according to the manufacturer’s specifications. For the extreme violet the 
intensities have been ‘‘ built on’”’, using the observed intensities in the 


Newtonian spectra. 


Tas_e II 


Anon 
Ion NGC 3918 165 13™-5 NGC 6302 NGC 6326 IC 4776 
—51° 52’ 


1 * 


-_ 
on 


3726°04 [O 11] 
3728°80 (O 11] 
3750°15 H1 
3759°87 Om 
3770°63 H1 
3797°90 H 1 
3819°61 He! 
3835 °39 H1 
3868-76 [Ne 111] 
3888°8 H 1-He 1 
3967°47 [Ne 111] 
3970°07 Ht 
3994 a 
4026°1 He 1-He 11 
4068-60 [S 11] 
4076°35 [S 11] 
4097°31 Ni 
4101°74 H 1(+ N 111) 
4120°8 He 1-O 11 
4177" [Fe 11] 
4199° He t1-N 11 
4207°: —e 
4244" [Fe 11] 
4267-2 Ci 
4287° [Fe 11] 
4340° H1 
4363° [O un} 
4387° Her 
4415° O n-[Fe 11] 
4471°48 He! 


We nN 
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TaB_e [I—(continued) 


Anon 
Ion NGC 3918 16413™-5 NGC630z2 NGC 6326 IC 4776 
—51° 52’ 


N itl aie ? 
He 11 
Si? 


Nl 
N u1-O 
Ciun-O1 
[Fe 111] 
He 11 
[Fe 111]-O 11 
[A tv]-He I 


[Fe 111] 
[A rv] 
[Fe 111] 
H1 
Ol 
He I 
[O 111] 
{O 111] 
Hel 
5197°94 (N 1] 
5200°41 [N 1] 
5207°2 aa 
5754 [N 11] 
5978 ? Sill 
6242 
6253 “oe 
6300 {O 1] 
6312° [S 111] 
6363°8 [O 1] 
6403°7 ~ 
6435" [A vj 
6461 _ 
6493 eee 
6548°06 [N 11] 
6562 °82 H1 
6583°37 (N 11] 
6678-15 Hel 
6716-42 [S 11] 
6730°78 [S 11] 


‘57 

5875°6 Hel 
"4 
BS 


NGC 3918, 11" 47™-8, —56° 54’ (1950), /=262°-6, b= +4°°4. 

This object has been discussed by Evans and Thackeray (7), who give a 
photograph (Fig. 9, Plate g). It was then found to be a simple, round, very 
bright object, of diameter 13”, with no detectable central star. The spectra used 
in the present study show doubling of the lines on high dispersion. Ordinary 
measures of the separation on three a-dispersion plates gave a mean expansion 
velocity of 20-1 km sec. On one plate the whole lengths of AA 3686 and 3967 
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and, on another, of the former line were so well defined that it was possible 
to measure the shapes in detail. It was found that the shapes were consistent 
vith the uniform expansion of a shell of diameter 18”-4 with a velocity of 
19-9 km sec. ‘The diameter so inferred is rather greater than that observed, 
ither on the Newtonian plates or in terms of the line-lengths of the Cassegrain 
spectra, and, indeed, the double lines fade away at their ends before forming 
a complete lens shape. On the Newtonian spectrum there is a continuum 
detectable down to approximately 3550 A; it shows a fine structure—too fine 
to be measurable—and a rapid decrement, interpreted as the head of the Balmer 
series. ‘The numerous lines recorded on this spectrum, and the very small 
scale, make it impossible to be dogmatic, but there is nothing in the continuum 
vhich can definitely be identified as an absorption line originating in the spectrum 
of a central star. The “ line-strengths” of Table II suggest an object of 
fairly high excitation. 


Anon 16" 13™-5, —51° 52’, /=299°:3, b= —2>1. 

This system is also included in Evans and Thackeray’s list ((7), Fig. 15, 
Plate 11), but is far too faint to be within range of the Cassegrain spectrograph. 
Recently, two Newtonian spectra have been obtained. It will be seen from 
the line list of ‘Table II, that this object has the V, line only just detectable, 
and N, less than one-tenth as strong as HB. Most of the lines in the spectrum 
seem to be attributable to forbidden Fett, and the degree of excitation is 
phenomenally low. An example of the trouble given by blends may be instanced 
here. The [Fe 111] line at A 4658 is twice as strong as .\,, and it may be that 
the latter is badly blended with other lines of the same multiplet occurring 
close enough to it to form a blend at the very low dispersion which obtains. 
For any ordinary planetary, it would be fantastic to feel any such concern, but 
not in the case of this candidate for what has been called ‘‘ The Chamber 
of Horrors ” of planetary nebulae. 


Fic. 1.—Anon, 160 13M-5, —51° 52’. Cf. Evans and Thackeray ((7) Fiz. 15, Plate 11). 


As seen in the telescope this grisly creature seems to consist of no more 
than two arcs of nebulosity to north and south of the central star. A longer 
exposure such as that shown in Evans and Thackeray (7) plumps these out 
into two blobs of nebulosity, one of them being pointed. By an accident, 








Montuty Notices of R.A.S. Vor. 119, Prare 5 


1950 September 13, 
103a0, no filter, 7™. 


Scale: 3°-4 mm. 


lugust 8, 103aE +-Wratten filter 


3 "9g mm 


David S. Evans, Notes on five southern gaseous nebulae 
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this system has been included at the edge of two deeper plates taken on 1958 
August 16 by Dr Thackeray and, particularly on the red plate of the pair 
(103aE plus GGrr1 filter), it can be seen that the creature has legs which 
complete the ant-like form (Fig. 1). The dotted lines superposed on the blobs 
indicate the form as shown on short exposures. The dotted area to the north 
indicates a region of doubtful nebulosity. 


NGC 6302, 17" 10™-4, — 37° 03’, /=317°1, b= —0°°3. 

[his is a very remarkable diffuse nebula, practically in the galactic plane, 
and only io from the galactic centre. ‘The system is illustrated in the 
accompanying plates. ‘The small sized reproduction shows a short exposure 
blue plate, from which it can be seen that the system comprises two lobes of 
luminous gas, shaped like the wings of a butterfly, separated by a relatively 
dark lane. There is no central star, but there are a number of stars involved 
in the system, suggesting that it is much larger and more distant than an ordinary 
planetary. The two other plates are of longer, matched exposures, and were 
taken under conditions of superb seeing. Reference to the list of line-strengths 
shows that, although Hz is very strong, A 6583 is stronger still, and must really 
be responsible for the major part of the density on the red plate. This pair 
of plates repays closer study. ‘The two halves of the system incorporate a series 
of opposed shells or fronts, and the regions of greatest luminosity lie within 
these fronts near the point where they are closest to each other. Behind these 
fronts lie a series of wisps of gas which, in the red plate, show fine filamentary 
structure, strongly reminiscent of that seen in solar prominences. Evidently 
the whole region contains a good deal of gas which, in the outer parts, is lighted 
up in a patchy way. For example, two filaments form the letters “ON ” 
near one of the bright stars in the field, and there is another faint wisp Sf 
from this. ‘The impression conveyed is that of dynamic interaction, as 
exemplified by the star immediately to the south of the centre of the system. 
On both plates this shows a patch of luminosity “ banked up ”’ on its preceding 
side. Another striking feature is the large number of extremely red stars 
which appear. The two prints are well-matched as far as the brighter stars 
are concerned, but half a dozen stars can be picked out by inspection which 
are bright on the red plate but almost invisible on the blue. 

All the spectra have been taken by guiding on the brightest part of the 
nebula, with the slit east and west. The lines often resemble skittles showing 
the intensity drop at the centre. On the Newtonian spectrum, where there 
is no dekker for the nebular spectrum (only for the comparison), the extreme 
length of the lines corresponds to about 110”. The overall diameter on the 
prints is 3 to 4 minutes of arc. There is no distortion of the lines on the 
Cassegrain spectra due to differential radial velocity which could profitably 
be measured, and it is thought that deviations from the mean could not exceed 
at most 5 km sec. ‘The lines in the violet were only detected on the Newtonian 
spectrum, where the oxygen pair at A 3726-28 is blended. On this film, 
continuous emission was detected, extending from about 3900 A to the limit 
of film sensitivity (103aQ). 

Since the form suggests dynamic interaction, but there is no marked radial 
velocity change along the lines, it may be inferred that the motions are 
predominantly transverse. A possible interpretation is that of a collision 


12 
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with the production of shock fronts. It may be profitable to search for radio 
emission from this system, which, so far as the author is aware, has not 


previously been noted as a radio source. 


NGC 6326, 17" 16™-8, —51° 42’, /=305°-6, b= —g°°5. 

This object has been described by Evans and Thackeray ((7), Fig. 18, 
Plate 11). It may be a ring nebula with a bright spot on the ring, or a second 
star associated. Otherwise it is not particularly remarkable. ‘The panchromatic 
plates have not shown up any lines in the red except Ha, and it does not seem 
as if there is a real paucity or weakness of lines in the red. The Newtonian 
spectrum shows a continuum starting well below A 3727. 


IC 4776, 18" 42™-6, —33° 24’, /=329 -4, b= —14°'9. 

This object was described by W. H. Wright (8) as “a small bright 
irregularly elliptical nebula”’. Wright identified nine lines in its spectrum. 
The spectrum is very rich: the excitation is low. On the Newtonian spectrum 
the continuum due to the central star can be traced throughout the film 
sensitivity, (103aO), but no detail can be made out. 

Acknowledgments.—The collection of the observational material has not been 
easy, and | am much indebted to the skilled assistance of Mr J. Churms and 
Dr P. A. Wayman in this task. My thanks are also due to Dr A. D. Thackeray 
for drawing my attention to his recent photographs of the anonymous system, 
and to Dr T. D. Kinman, who, both in personal discussion and by correspondence 
has kept me informed of his extremely careful and illuminating studies of the 
performance of the Newtonian spectrograph. 

Note.—A referee has drawn my attention to the paper of L. H. Aller, 
I. S. Bowen and R. Minkowski on the spectrum of NGC 7027 (Ap. 7., 122, 
62, 1955). The nebular wave-lengths adopted there represent a further slight 
revision of those given by Bowen (3). ‘Tests of the reduction of our spectra 
using these wave-lengths show a maximum change in the derived velocity 
of o-6 km sec and almost invariably a slight reduction in the probable erros 
of the result. The changes are never significant in comparison with the 
probable error of the mean. So far as identifications go, there seems no need 
for change except to split the blend at 4712 A into 4711-34, [A iv] and 
4713°17, Hel, in fact assigning the tabulated lines to the former. It will 
probably be preferable to use the new wave-lengths in future work. 


Royal Observatory, 
Cape of Good Hope: 
1958 November. 
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GLOBULAR CLUSTERS. I. THE RADIAL VELOCITIES 
OF SOUTHERN GLOBULAR CLUSTERS 


T. D. Kinman 


(Communicated by the Radcliffe Observer) 
(Received 1959 February 6) 


Summary 


The radial velocities of 30 southern globular clusters are given and their 
errors are discussed. Spectra were obtained for 28 clusters using the 
Radcliffe Cassegrain spectroscope (86 A/mm) comprising spectra in integrated 
light for 15 clusters, spectra of individual stars for 10 clusters and spectra 
of both types for 3 clusters. No previous velocities were known for 16 of 
these clusters. Spectra were obtained with the Newtonian spectroscope 
(380 A/mm) in integrated light for 4 clusters, for 2 of which Cassegrain 
spectra were also available. 

A comparison of the present velocities with-those obtained by Mayall (5) 
and Joy (3) for the 12 clusters in common shows satisfactory agreement and 
a table of the velocities and distances of 70 globular clusters is given. 





Introduction.—A study of the motions of the globular cluster system of the 
Galaxy is of particular interest because these clusters are prototypes of Baade’s 
Population II and are believed to contain the oldest stars in the Galaxy. 


N. U. Mayall (5), in a remarkably fine investigation with the Crossley nebular 
spectroscope (430 A/mm at Hy), obtained velocities for fifty of these clusters, 
which is rather less than half of the known total. The distribution of the 
globular clusters in the sky is unfavourable to the northern observer; although 
Mayall managed to observe as far south as 6= — 40°, he noted that there were 
23 clusters south of this limit, of which 17 might be observable spectroscopically. 
Many of these southern clusters are in the region of the antapex of the Sun’s 
galactic motion and it is therefore important that their velocities should be 
observed in order to make a reliable determination of the solar motion with 
respect to the globular cluster system. An attempt has therefore been made 
to observe as many of these southern clusters as possible using the Cassegrain 
and Newtonian spectroscopes of the Radcliffe 74-inch reflector; in addition, 
twelve clusters for which either Mayall (5) or Joy (3) give velocities have been 
observed in order to obtain the velocities on a uniform system. The Cassegrain 
and Newtonian spectra are described and their reduction and errors are 
discussed in Sections 1 and 2 respectively of this paper. In Section 3, the 
radial velocities and distances are listed for 70 globular clusters. The spectral 
types obtained from the Radcliffe observations will be given and discussed in 
a second paper, while the analysis of the radial velocities of the 70 globular 
clusters is to be’ given in a third paper. 
1. Determination of velocities with the Cassegrain spectroscope 

1.1. Observations.—All spectra were taken using the shortest focal ratio 
camera of the Radcliffe 2-prism Cassegrain spectroscope (86 A/mm at Hy) (1) 

12* 
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with a slit width of 0-20 mm which projects to 27 » on the plate. In general, 
it the core of a cluster is sufficiently bright and condensed to obtain a spectrum 
in integrated light then the individual stars are too faint for this spectroscope; 
conversely when the stars are bright enough to observe individually, the 
cluster is too extended to obtain integrated spectra. Only for four clusters 
were both types of spectra obtained and for one of these (47 Tuc) the stellar 
spectra were obtained by M. W. Feast and A. D. Thackeray and will be 
discussed by them elsewhere. 

Sixty-nine spectra of 18 clusters were obtained in integrated light. The 
cluster cores were in general sufficiently unresolved for it to be unnecessary 
to put the image out of focus on the slit (cf. Mayall (§) and Morgan (7)); 
satisfactorily uniform spectra could be produced by trailing the image along 
the slit and by scanning in a perpendicular direction. Only for four clusters 
(NGC 104, 6397, 6752 and 7078) were the cores sufficiently resolved for it to 
be necessary to have the image out of focus. Projected slit heights of either 
0°50 or 0:25 mm were used, depending upon the concentration and diameter 
of the cluster; exposures ranged from about one to six hours. 

Sixty-four spectra of forty-one stars in 13 clusters were obtained; spectra 
of at least 3 stars per cluster were taken if the velocity of the cluster was 
previously unknown. Stars were selected from the cluster giants which under 
a power of about 300 often may be distinguished by their reddish colour. ‘To 
increase the probability of finding a cluster member, stars close to the cluster 
centre were chosen, rejecting, however, any not well separated from fainter 
stars. Of the spectra taken, only 3 were rejected as being field stars: their 
velocities and spectra are quite different from those of other stars in the cluster 
concerned. ‘The remaining stars, with one exception, are clearly cluster giants; 
the exception is a blue star in NGC 6752(*4) with almost the same velocity 
as the cluster. While its membership remains dubious, its inclusion has not 
affected the final velocity of this cluster. ‘lhe projected slit height of the 
stellar spectra varies from 0-25 mm for the brightest stars observed (m,,, ~ 12-6) 
to o-10 mm for the untrailed spectra of the faintest (m,,~14°5); the exposure 
times range from about one to six hours. Sixty-nine spectra of 48 standard 
stars of spectral types Fo-Ks5 were taken concurrently with the cluster 
programme. Forty-one of these stars were of “‘a’’ quality in the Mount 
Wilson General Catalogue (12) and seven were stars for which the Cape kindly 
supplied well-determined velocities. All of these stars were observed within 
25 of the zenith to reduce atmospheric dispersion errors, and neutral filters 


were used so that the exposure times were not less than 10 minutes. 


are Comparison spectrum was impressed on either side of all spectra 

it the begining and end of each exposure and in the case of the cluster spectra 
at intervals during the exposure. All spectra were taken on Kodak 103aO 
emulsion and developed in a metolquinone developer for 5 minutes at 67 °F. 
Veasurements.—A slight curvature of the field of the F 2 camera lens 
spectroscope limits the spectrum range available for measurement to 

it 400 A in the photographic region. The range AA 4400-4000 (centred 

on the wave-length of minimum deviation and the central ray of the photographic 
plate) does contain, however, almost all the features which are suitable for 
measurement in these cluster spectra. At the correct plate tilt, the departure 
F 1 


from a flat field over this range does not exceed the accuracy of focusing 
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(+15). All spectra were measured once (direct and reverse) on a standard 
Hilger photomeasuring micrometer. A few repeat measures gave almost 
identical results to the initial measures with similar residuals for given lines. 
This agrees with Mayall’s conclusion (5) that the errors of measurement are 
smaller than inherent plate errors so that re-measurement is unnecessary. 


TABLE I 
Adopted wave-lengths for iron arc spectrum (86 A/mm) 

3956°681 4030°755 4132°060 4282-406 4315 °087 4383°547 

3969-261 4045°815 4202°031 4294°128 4325°765 4404°752 

4005 246 4071°740 4260°479 4307°906 4337°049 4415 °125 

Tas_e II 
Measured stellar wave-lengths (86 A/mm; projected slit width 27 py) 
A (Go) (Fs) (Ko) A (Go) (Fs) (Ko) 

4005 °67+0°08 0°00 0°00 4215 °66 +0°08 +0°24 —0°36 
4032°65 0°08 —0'07 +0°21 4226°52 0°08 —0'03 +0°07 
4045°60 0°08 — 0°04 +008 4254°93 0°14 0-00 0°00 
4063°53 0°03 —0'07 +0°05 4260°39 ‘II — 0°04 +0°06 
4067°27 ovll +0°03 +0°26 4272°31 ‘11 —o'18 +0°24 
4071°70 0°03 —o'ol +0°14 4290°02 "14 —0°27 + 0°27 
4077°38 0°08 +O'l4 +0°08 4325°30 "14 +0°'29 +0o°o1 
4101°80 0:08 +0°08 —0'23 4339°89 "12 +0°22 — 0°33 
4133°03 O°14 — 0°04 +0o°19 4352°10 "12 +0:'06 +0°09 
4143°33 Ov1l +0°06 +0°08 4375°06 o-'12 +0°'09 — 0°07 
4172°81+0'11 —0'08 +0'18 438411 +015 +0o°I9 —o'1g 
Neither iron arc nor stellar wave-lengths are available for the present slit 
width and dispersion. In the first place, therefore, all visible features in the 
stellar spectra and apparently suitable unblended comparison lines were 
measured for a number of standard stars. Using I.A.U. iron arc wave-lengths 
as in (1), Hartmann dispersion curves were found for the iron lines in these 
spectra and those lines which consistently stood off the correction curves were 
rejected, leaving the 18 lines listed in Table I. Knowing the stellar velocities, 
mean provisional wave-lengths were then deduced for the features in the stellar 
spectra. The cluster spectra were then reduced following the method of 
J. S. Plaskett (10)*: standard micrometer readings being set up for both arc 
and stellar lines. The differences between the observed and standard readings 
for the stellar lines, when corrected from a curve derived from the arc line 
differences and multiplied by a factor found from the Hartmann formula, give 
the velocities directly. Similarly, velocities were found for the remaining 
spectra of standard stars observed concurrently with the clusters and from these 
a correction curve (in km,s) was plotted for twenty-two stellar lines as a function 
of spectral type. These corrections to individual lines were then applied to 
the cluster spectra. This procedure enabled almost all the features in the 
cluster spectra to be used, even though several features showed a strong 
wave-length dependence on spectral class; only the G-band components were 
invariably rejected from the features clearly visible in the range AA 4400-4000. 
Col. 1 of Table II lists the final wave-lengths of the 22 features used for class Go 
together with an assessment of their internal consistency. Cols. 2 and 3 give 


* The controlling temperature of the Radcliffe spectroscope is varied over a range of about 
10 °C between summer and winter; the resultant change in the spectrum scale (~1°5 «% per mm) 
is small enough for only one standard Hartmann curve to be necessary. 
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the corrections to these wave-lengths for classes F5 to Ko and may be linearly 
interpolated. 
The heliocentric cluster velocities are given in Col. 5 of Table III. 


/ 


TABLE III 
Radial velocities of globular clusters from Cassegrain spectra 


NGC number of cluster; J=integrated spectrum, *=spectrum of individual star 
(spectra of the same star are bracketed) with identification number. The positions 
of these stars are shown on the finding charts. 

Date. 

Exposure time in minutes (t). 

Number of lines measured (nm). 

Heliocentric radial velocity in km/s (V). 

R.m.s. error of this velocity in km/s (derived from line residuals in the case of 
individual spectra and from plate residuals in the case of mean velocities) (c). 
Assigned weight on Mayall’s system (W’). 


NGC Date t ; W 


104 


1957 7 57 
1957 5 113 
1957 9 60 
1957 10 6 97 


1957 


4°4 
4°4 
4°4 
4°4 
4°4 


272 


Cou 


ww 


1956 
1957 
1956 
1957 


1957 


1950 7 2 Il 
1956 I 
1950 °¢ 5 I 
1956 : ; 
1956 2: 13 
1957 15 


Mean 


1950 
19560 


1956 


1956 
1950 
1956 


1957 


1956 
1957 
1957 


1957 
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TaBLe III (continued) 


t n 
180 12 
210 15 
240 8 
217 8 
300 II 
183 16 


ut ouMr hb OO 


265 
270 
270 
180 
180 
130 
180 
180 


pW AUUND SD 
SUNIIUM AAD 


185 


1956 
1957 
1956 
1957 
1956 


1956 
1956 
1957 
1957 


1957 


> MU Qtr nN 


1956 
1956 
1956 
1957 


tN 
“I nv oO 


+ 
Ass 


38 
35 
45 


_ 
aw 


oe 
° 


5904 1957 + 42 


6121t 1956 29 17 + 58 
ki 1957 19 00 14 73 
me 1957 26 18 77 
“? 1937 26 7 14 57 
Mean 66 5 


+In NGC 6121, stars 1, 2 and 3 correspond respectively to Greenstein’s 113, 259 and 
261 (Ap. F., 90, 401, 1939, in which a finding chart is given). 
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Taste III (continued) 
t 
360 
223 
310 


1956 
1956 
1950 
1956 


1957 


1957 


1957 


+ 
1?) 


{1956 
< 1956 
| 1958 
(1956 
< 1958 
1955 

1957 

1957 


x 
° 


* 


i i 
nN +> 


—_ 
~ 


* 


* 


moo0ow ad 


NNN & N WN 
o) 

— 

“I Gur ui WwW 


ree 


I 
2 
2 
2 
3 
4 


1956 
1956 
1957 
1957 
1957 


1957 


1956 
1957 
1956 
1957 


1956 
1956 


1957 
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Taste III (continued) 

t n 
195 13 
240 II 
155 13 
220 15 
185 16 
IIo 16 
120 15 
10 

92 

190 

60 


AunovFr DOOWAO wo 


ee 


6 
5 
8 
9 
7 
4 
6 
8 
6 
7 


5 
13 


12 190 


1956 155 — 196 17 
1957 145 — 204 4 
1957 155 — 169 28 
1957 165 —171 3 


Mean — 185 9 9 


1.3. Errors in these velocities. (a) Wave-length errors.—The Royal Observa- 
tory, Cape, have kindly supplied wave-lengths determined by them for the 
same dispersion but with half the slit width. Of the 12 wave-lengths in common, 
7 agree within the quoted Radcliffe errors. The difference for the other 5 lines 
is about twice the quoted Radcliffe error and may be due in part to the higher 
resolution of the Cape spectra. Synthetic wave-lengths were calculated for 
these 5 lines from Rowland revised wave-lengths (11) of the lines in the 
corresponding waveband to the measured line widths. The Rowland wave- 
lengths were weighted according to the equivalent widths derived from Mulders’ 
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Finpinc Cuarts (1).—The NGC number of each cluster is given at the bottom left hand corner of 
each chart. The star numbers are the same as those given in Table III, col. 1. 
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calibration (8). ‘These synthetic wave-lengths show agreement with the 
Radcliffe wave-lengths within the quoted error. 

The error in the velocity of a line in a cluster spectrum produced by wave- 
length errors ranges from +6 to +10kms compared with observed r.m.s. 
error of a single line in a cluster spectrum of +25 kms. The variation of 
wave-length with spectral class is in different senses for different lines and an 
error of 0-3 of a class will produce a velocity error of about 2-3 kms. 
Systematic errors derived from the adopted wave-lengths in the final velocities 
are therefore likely to be satisfactorily small. 

(6) Instrumental errors.—The satisfactory performance of the spectroscope 
at higher dispersions (1) suggests that only errors which increase with decreasing 
camera focal ratio (F) need be considered. The Hartmann shift arising from 
an error of focus (€ microns) together with a difference between the centroids 
of illumination of the collimator by the star and by the comparison arc of a 
fraction p of the collimator diameter will produce an error 12ep F? kms in 
the Radcliffe spectroscope. Uniformity of the stellar illumination of the 
collimator is checked periodically, while the arc illumination was adjusted to 
be uniform visually before each night’s observations, so it 1s thought that p 


is unlikely to exceed o-o5. ‘The error in focusing is of the order of + 15 microns 
Ss. 


Since 
the cluster velocities are in effect determined differentially with respect to 


so that the error from Hartmann shifts can hardly exceed 2-3 km 


standard stars, any systematic error from this cause is likely to be much smaller. 
\ similar argument applies to errors from slit curvature which also increase 
with decreasing camera focal length. It may be shown on the basis of 
Wesselink’s measures at higher dispersions (1) that this error is unlikely to 
exceed o-1 km s and so is quite negligible. 
) Guiding errors.—The slit width used corresponds to 162 kms at Hy. 
his width is a compromise between the need to match the plate resolution 
ind the need to obtain maximum light grasp. The present slit width is 
equivalent to 1”-2 at the Cassegrain focus and so enables the spectroscope to 
accept a large fraction of the stellar image under conditions of good seeing 
which frequently obtain. Considerable errors in velocity can clearly arise 
however unless the star is accurately guided on the slit. ‘T'wo cases may be 
distinguished : 

(i) Stars observed near the zenith. A sample of 16 spectra of G and K type 
stars is considered (from the standard stars used to determine wave-lengths) 
which were observed within +10 of the zenith. The mean value of the 
residuals: General Catalogue — Radcliffe o-5+i:t1kms with an r.m.s. 
error for a single spectrum of +4:5kms. ‘The corresponding figure at 
49 A mm and slit width 0-075 mm for this spectroscope is about + 2-0 kms, 
so that the guiding errors which result from using a wide slit near the zenith 
must be quite small. 

(ii) Stars observed away from the zenith. For these stars the stellar image 
is drawn out by atmospheric dispersion into a small spectrum which will be 
perpendicular to the slit for a star on the meridian. Colours in this atmospheric 
spectrum (on the slit of the Radcliffe spectroscope) are indistinguishable for 
stars fainter than about roth visual magnitude. For fainter stars the observer 
may therefore be expected to guide with the visual wave-lengths of the image 
centred on the slit and with the photographic image systematically displaced 
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from the slit so that a guiding error will result. The spectrum lines will be 
shifted to the blue for objects south of the zenith in the case of the Radcliffe 
telescope and spectroscope in their normal orientation. The effect should 
increase with increasing zenith distance and increasing slit width; dt should 
decrease with increasing size of the seeing image. 

Velocities obtained from integrated spectra are not subject to this error 
because an extended source is used. Consequently one might expect to find 
a difference (A) between the velocities found from the integrated spectra and 
from the spectra of the individual stars of a cluster; this difference should 
increase with the zenith distance ¢ of the cluster. Mean values of ¢ and A are 
given for three clusters in Table IV. It is seen that the differences A are smaller 
than their r.m.s. errors and show a slight change with zenith distance in the 
opposite sense to that expected. It may be concluded therefore that no 
significant systematic errors due to guiding exist in the velocities discussed 
in this paper. 

TABLE IV 


The difference (A) between the cluster velocity obtained from individual stars 
and from the integrated cluster light as a function of zenith distance (€) 


NGC 6397 6752 104 
4 28°S 35°5S 47° S 
A —4+16 —I1+7 +4+6 


1.4. Comparison of cluster velocities with those from other sources.—Mayall 
(loc. cit.) gave weights to the velocities found from the spectra of different 
observers which are inversely proportional to their dispersion. Following this 
procedure, the present spectra are given weight 4-4 on Mayall’s scale on which 
each of his low dispersion spectra (430 A/mm) have unit weight. Those 
Radcliffe spectra which contain only 6 lines or less have been given half this 
weight (2:2). The mean velocities obtained by this simple weighting in fact 
differ little from those which would be obtained by weighting according to the 
square root of the number of lines in each spectrum. 

The weights used show only a weak correlation with the r.m.s. errors 
derived from plate residuals (¢), which perhaps is hardly surprising because 
of the small number of plates available per cluster. If it is assumed that the 
weight should be proportional to o-*, then the Radcliffe spectra have been 
slightly underweighted with respect to Mayall’s spectra by taking only the 
dispersion into account; this presumably is because of the larger number of 
lines in the Radcliffe spectra. Since, however, this effect is difficult to allow 
for satisfactorily, Mayall’s system of weighting has been retained, 

Table V gives a comparison of the velocities for twelve clusters which 
have been observed in common with either Mayall or Joy (3) together with 
their r.m.s. errors and weight on Mayall’s system. ‘The mean difference 
(Radcliffe-Mayall) is +4+8km/s and the average difference without regard 
to sign is 20kms. ‘This may be regarded as satisfactory since, although the 
ditferences for some clusters (NGC 1904 and NGC 6864) are greater than 
would be expected from the internal r.m.s. errors, similar comparisons by 
Mayall suggest that his adopted velocities may be in error by several times 
their internal error (i.e. from +15 to +20km/s). If we discount Joy’s single 
spectrum of a variable star in NGC 5139 (wCen) the mean difference 
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to sign is 11 km,s, which is also satisfactorily small. ‘The final column of 
‘Table V gives the mean adopted velocities with r.m.s. errors based upon the 


(Radcliffe-Joy) equals +3+7 km/s and the average difference without regard 


weighted velocities of the ditferent observers. 


TABLE V 
A comparison of globular cluster velocities (in km/s) 


NGC Radcliffe Mayall Joy R-M R-J Adopted 
(R) (M) (J) velocity 
1851 +317+ 3(18) +2914 8 (8) ane +26 soe + 309 + 12 (26) 
1904 +181 8 (13) +228 10 (6) Eo —47 wd +196 22 (19) 
5139 +231 4 (22 des fo + 5 +230 8 (22) 
5904 42 10 (4) 5 (tg) + 49 3 (42) 
6121 66 5 (18) com - ; = ; 65 2 (26) 
6266 61 (4) 8 17 (9) ean us 7% 6§6¢ G3) 
6656 142 5 (15) § 6 8 (8) + 144+ 4 (34) 
6715 30 3(13) +107 9 (7) ai - e 122 2 (20) 
6864 187 (15) 8 (7) ea nee 198 17 (22) 
7078 98 2(11) 7(15) —138+16 (3) + 3 107+ 8 (29) 
7089 12 10 (7) 7417) -— 4+ 4(18) s+ 24) 
7099 -185 9 (9) 3. 10 (9) neh a 174+ 9 (18) 


2. Determination of velocities with the Newtonian spectroscope 


2.1. Observations.—A single prism spectroscope has been used at the 
Newtonian focus of the Radcliffe reflector to obtain low dispersion (380 A/mm 
at Hy) spectra of 2 clusters (NGC 1261 and NGC 6139), for which only a 


few low weight Cassegrain spectra were available, and also 11 spectra for 
2 clusters (NGC 5927 and NGC 6584) which are too faint for the Cassegrain 
spectroscope. ‘The few remaining clusters which are important for a solar 
motion study proved to be beyond the range of this instrument but the 
programme would in any case have been limited by the time available. In 
this spectroscope the light is collimated by an aluminized parabolic mirror of 
diameter 55mm. ‘The collimated beam is then intercepted by a perforated 
flat placed just behind the slit—an arrangement originally due to Pfund (9). 
Dispersion is produced by a single 60 dense flint prism and the camera consists 
of a detachable Schmidt system with uviol correcting plate and focal length 
45 mm nominally working at F\1. The effective focal ratio must be larger 
however because of substantial vignetting in the system. All optical components 
of this spectroscope except the dense flint prism were made by Messrs. Cox, 
Hargreaves and ‘Thomson. 

‘The H and K lines must be clearly recorded in order to obtain useful 
velocities at this dispersion. At this wave-length the transmission of the 
spectroscope after the slit can hardly exceed 15-20 per cent with the principal 
light loss occurring in the prism: the difficulty of obtaining spectra of faint 
reddened objects will therefore be appreciated. Nevertheless the mean exposure 
times for six clusters showed that this spectroscope is 4:1+0°9 times faster 
than 

A cluster was maintained on the slit by off-set guiding on a nearby star. 


1 


the Cassegrain spectroscope at its lowest dispersion (86 A mm). 


Since the cluster could not be viewed on the slit, a suitably integrated spectrum 
covering most of the dekker width was obtained by placing the slit 8 mm 
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outside the Newtonian focus. A comparison arc was provided from an He 
discharge tube which together with a diffusing screen could be swung in front 
of the slit. A straight slit of width 0-15 mm (projected width 28 4) was used 
and because of the short focal length camera the spectrum lines are therefore 
noticeably curved (radius ~5 cm). Spectra were obtained on disks of 103aO 
film, 0-85 cm in diameter, which were cut out by a special stamp. The film 
was processed in a perforated test tube and was then permanently mounted 
between two microscope slides. 

2:2. Measurement and reduction.—The mounted spectra were measured direct 
and reverse on a Hilger photomeasuring micrometer. ‘The reduction was 
performed using a standard Hartmann curve as for the Cassegrain spectra; 
M.I.'T. wave-lengths (6) were used for the He are and the Lick extra-galactic 
nebulae wave-lengths (2) for the cluster absorption features. ‘Twenty spectra 
of nine clusters with known velocities were used to derive corrections (in km/s) 
which were applied to individual lines to obtain the correct velocity. Since 
the calcium H line is blended with He, two corrections were found for this 
line: one for type F (NGC 1261, 6139 and 6584) and one for type G (NGC 5927); 
no correction for spectral type was needed for other lines. This correction to 
individual lines takes into account line curvature and other systematic errors 
such as may be caused by slit irregularities or differential Hartmann shifts 
between arc and cluster spectra due to non-uniform illumination of the 
collimator. As with the Cassegrain spectra, the velocities are essentially 
measured differentially with respect to those known velocities of objects of 
similar spectra. The accuracy obtainable from different lines shows considerable 
variation so weights were given which were proportional to the inverse square 
of their r.m.s. errors in the standard spectra. These weights are listed in 
Table VI; the difference between Lick and Radcliffe weights for the K line 
may be due to underexposure of the Radcliffe spectra at this wave-length. 
Final heliocentric velocities are given in Table VII. 


Taste VI 
Newtonian spectroscope wave-lengths 


Radcliffe Lick Mt Wilson 
A Wr A Wy A 
3933°8640°34 0°376 3933°28 384 3933°7 
(G) 3968°37+0°24 1:000 H 3968°38 301 3968°5 
(F) 3968-72 H(2)+He(1)* 3968-54 10 ree 
H(1)+He(z) 3969-01 37 
H+He(=) 3969°23 34 itis 
4102:06+0°38 0°376 4101'25 120 4101°7 
4302°76+1°06 0:°085 4393°52 140 4304°4 
4339°'17 0°90 O'l4I 4340°61 36 4340°5 


Wy and W, are the relative weights adopted by Radcliffe and Lick respectively. 
* Numbers in brackets give hypothetical intensities of the unresolved components. 
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Tasre VII 
Radial velocities of globular clusters from Newtonian spectra 
NGC number of cluster. 
Date. 
Exposure time in minutes (f). 
Heliocentric radial velocity in km/s (V). 
R.m.s. error of this velocity in km/s (derived from line residuals in the case of 
individual spectra and from plate residuals in case of mean velocities) (c). 
Assigned weight on Mayall’s system (W). 
NGC Date t : ; W 
1261 1957 10 2 280 : +5 o'sI 
1957 II 2 240 31 1*0O 
Mean (Newtonian) 3 1°8 
Mean (Cassegrain) 6-6 


Mean I 8°4 


185 0°74 
185 o's1 
190 0°74 
230 0°74 
360 0°74 
360 0°74 


Mean 4°5 
1955 8 14 150 10 
1958 8 14 145 12 


Mean (Newtonian) | I 1°9 
Mean (Cassegrain) 29 6-6 


Mean + 20 I 8°5 


1957 38 +114 - os! 
1955 + 210 . 0°93 


5 ; 
1955 5 5 +245 5 0°70 
1958 6 ‘ +184 0°93 
1955 6 + 42 y 0°89 


Mean +160 4°0 


2.3. Errors in these velocities.—I\t is hoped to publish a fuller description 


of the spectroscope and its performance later; consequently only the results 


of 


tests of various sources of error are summarized below. 


(a) Flexure and thermal stability.—A null result with a r.m.s. error of 
s km/s was found for flexure alone. The mean shift of the A 4046 Hg line 


m long exposures with respect to the laboratory value was 

gt25 Km)S. 
(6) Emulston shifts—For differential shifts occurring over short ranges 
mm), no difference was found between 103a0 Cassegrain plate and 
‘wtonian film. ‘The r.m.s. error corresponds to +25 km/s. 


(c) surement,— The r.m.s. consistency of measuring a 
ngle spectrum is +37 kms. 


The r.m.s. error of a single Radclitfe spectrum derived from plate to plate 


variation is +58kms, which is little larger than that derived from the 
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agreement between the lines within a spectrum (Table VII, col. 5). For the 
20 denser standard spectra the error is +50 km/s. Corresponding errors for 
other low dispersion spectra are: Lick (globular clusters (§)) +33 km/s; 
Mt Wilson (extragalactic nebulae (2)) best quality +39 km/s, poorer quality 
+52km/s; Lick (extragalactic nebulae (2)) +81 km/s. 

The errors in the present spectra may be considered to be largely random 
and the spectra may therefore be weighted according to line weights: spectra 
on which all listed lines were used are given unit weight on Mayall’s system. 
‘The adopted weights are given in the final column of Table VII. 


3. Velocities and distances for seventy globular clusters 

It has been suggested that in the present case weights should be assigned 
according to the dispersion rather than by using the r.m.s. errors. Apart 
from invalidation by systematic errors (which should be small for the Radcliffe 
velocities), r.m.s. errors afford a poor indication of the accuracy if the number 
of spectra is small. Application of the x? test shows that if the r.m.s. error 
derived from a sample of four spectra be o, then there is a 20 per cent chance 
that the true value of this error in the parent population exceeds 20 and a 
7 per cent chance that it exceeds 30. The r.m.s. errors of Mayall’s velocities 
of weight four (four spectra) range from +5 to +53 km/s; weighting by the 
standard method according to o~* is thus inadmissible. 

The probable errors have in general been assigned according to the weight 
on the following scale: 


Weights 4- 6 (low dispersion only) P.E. +20 km/s 
8-10 15 
10-20 10 
20-30 7 
above 30 5 


These values have however served only as a general guide: each cluster 
has been considered individually taking into account the agreement between 
different observers. 

Tasie VIII 


Velocities and distances for seventy globular clusters 


NGC V W P.E. R (kpc) 

104 — 24 (22) a 4°8 5°8 6°9 
288 — 47 (22) | 9°6 12°6 16°6 
362 +221 (24) y | 9°1 12°6 17°4 
1261 + 46 ( 8) 15 19 29 44 
1851 +309 (26) 8 11‘o 16°6 25 
1904 +196 (19) 15 9°6 13°2 18-2 
2298 + 64 ( 5) 20 30 46 
2419 + 14 ( 5) 20 69 87 
2808 +101 (26) 7 . : II's 
3201 +493 (35) 5 : ‘ 6°3 
4147 +191 ( 4) 20 we 
4372 + 66 (13) 10 . . g'I 
4590 —116 ( 4) 20 3° 15°8 
4833 +204 (13) 10 . % 72 
5024 —112 (10) 15 . 25 
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Tasie VIII (continued) 
IGC 4 W PE. 
13° + 24 (22) a 
27 57 (54) 
5256 . ; (13) 
5634 ( 4) 
5694 ( 4) 
5324 5! ( 4) 
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Table VIII gives the radial velocity in km/s (V), weight (W) and probable 
error (P.E.) for seventy clusters. The heliocentric distances of these clusters 
in kiloparsecs (R) found in a previous discussion (4) are also given; the second 
figure in this column is the estimated distance, and the first and third are the 
limits corresponding to the probable error of the moduli. 
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THE BRIGHTNESS DISTRIBUTION 
WITHIN THE RADIO SOURCES CYGNUS A (19N4A) 
AND CASSIOPEIA A (23N5A) 


R. C. Fennison and V. Latham 
(Received 1958 December 8) 


Summary 


Measurements are described of the amplitude and phase of the Fourier 
transform of the distribution of brightness at 127 Mc/s within the intense 
radio sources Cygnus A (19N4A) and Cassiopeia A (23N5A) along an axis 
in the East-West direction. ‘The brightness distribution is computed from 
the observed reaaings and in the case of Cygnus A is shown to confirm the 
existence of two major emitting regions of almost equal intensity. The 
x component, leading in right ascension, is slightly brighter than the B 
component. Measurements are also described of the determination of the 
position of the major axis of the source. It is evident that the diameter 
of the Cygnus source in a position angle close to 187° must be relatively 
very small. The measurements on Cassiopeia indicate the presence of a 
faint spur near position angle go together with probably slight limb 
brightening of the main body of the source. 





1. Introduction.—The existence of asymmetry in the radio source Cygnus A 
(1g9N4A) was shown in the initial measurements of Hanbury Brown, Jennison 
and Das Gupta (1). ‘These measurements indicated that, to a first order, the 
shape of the source was elliptical, approximately two minutes ten seconds of 
arc in diameter in position angle 113° and thirty-five seconds of arc in position 
angle 177. ‘The major axis appeared to lie between position angle go° and 
position angle 120. ‘hese dimensions were not compatible with those of the 
photographic object observed by Baade and Minkowski (2) within the position 
rectangle given by Smith (3). The photographic object, interpreted by Baade 
and Minkowski as a collision between two late-type spiral nebulae, is elliptical 
in form, measuring 18” x 30” with the major axis in position angle 150°. The 
nuclei of the two galaxies are separated by 2” in position angle 115°. The 
distance of this object, obtained from the red shift of a neighbouring galaxy, 
is approximately 100 M parsecs*. ‘This distance has been confirmed in the 
radio frequency spectrum by the hydrogen line absorption measurements of 
Lilley and McClain (4). 

A more detailed study of the structure of the object at radio frequencies 
was performed by Jennison and Das Gupta (5, 6) who showed that the source 


consisted of two major components, each approximately 50” x 30”, separated 


along the direction of their major axis by 1’ 25” between their centres. On 
the assumption that the position and identification of the source were correct, 
this model straddled the photographic object and required the generation of 
the radio frequency radiation to take place at a considerable distance from the 
main mass of the nebula. The interferometer system used by Jennison and 


* Using a value for the Hubble constant of 180 km/sec/Mpsc. 
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Das Gupta was not suitable for the measurement of the phase of the Fourier 
components and the model given was not unique, though other solutions in 
general increased the overall size of the radio object. It was evident that a 
more refined technique was required for work on this problem. 

The intense radio source Cassiopeia A (23N5A) has been identified (2) as 
an unusual emission nebulosity within the Galaxy. The structure of the 
Cassiopeia source which appears on photographic plates taken at the 200-inch 
telescope, consists of a number of small filamentary structures and condensations 
distributed over an area five minutes of arc in angular diameter. To the north 
of the area there is a brighter arch of nebulosity of bluer colour. The diffuse 
filaments exhibit large internal random motions and are extremely transient; 
their structure has been observed to change over a period of a few years. The 
filaments also exhibit signs of an ordered expansion which indicates an origin 
280 years ago if the expansion velocity is still constant (8). 

The angular diameter of the Cassiopeia A source was first determined at 
radio frequencies by Smith (g) in position angle go°, and by Hanbury Brown, 
Jennison and Das Gupta (1) in three position angles inclined to each other 
at 120°. These measurements showed that the object was approximately 
radially symmetrical and had a diameter of the order of five minutes of arc. 

Further measurements by Jennison and Das Gupta in position angle 113° 
extended the resolution to baselines of over 3.000 wave-lengths (6). It was 
found that the correlation function of the radiation from the Cassiopeia source 
fell off rapidly up to a spacing of goo wave-lengths, thereafter apparently 
exhibiting a tail which approached zero in the neighbourhood of 1 500 
wave-lengths. No subsidiary maxima were observed but, as the correlation 
function measured by the equipment corresponded to the square of the 
Fourier’ transform obtained by conventional interferometer systems, the 
equipment was not sensitive to small readings and it was fundamentally incapable 
of determining the phase associated with such maxima. In order to obtain 
a more detailed picture of the radio frequency distribution of brightness of 
these objects, a high resolution, phase-sensitive interferometer was constructed 
and used in a series of measurements on the source at a frequency of 127 Mc/s 
during the period 1955-1956. 

2. The equipment and method.—A three-station interferometer system was 
developed in which most of the errors usually associated with the equipment 
and aerial siting were nullified, and both the amplitude and phase of the Fourier 
components could be measured concurrently over large spacings of the aerial 
systems. The equipment and the method of measurement have been more 
fully discussed in a previous paper (7). 

The measurements were made by spacing three aerial systems A, B and C 
and their associated relay transmitters in a line along which the readings were 
required. The signals received from the three stations were then combined 
in pairs by means of three switched interferometer systems in which the 
arguments of the fringes could be controlled by electrically-driven phase-rotating 
condensers. By spacing the aerials at equal distances the measurement of 
amplitude could be made absolute whilst relative measurements could be 
obtained at other spacings. ,The measurement of phase was relative to a datum 
obtained from measurements at short spacings and involved a cumulative error 
which increased with aerial separation. 
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Under the illumination of a complex source, the phase angle given by 
subtracting the sum of the arguments of the fringes formed between 4 and B 
and between B and C from the arguments of the fringes formed between A 
and C is given by the expression 


$=€ ac (f4p +£ pc) 


where £,,, €n¢ and €,,, represent the phase of the Fourier transform at spacings 
AB, BC and AC respectively. 

This result is independent of the phase errors and the angular velocity of 
the phase rotation. It enables an estimate to be made of the relative phase 
contributions due to the source structure only, and is applicable when the 
aerials are spaced by arbitrary distances. In the general case, when the transform 
is complex but the aerial spacings are unequal, the phase change over any 
extension of the baseline may be determined and normalized against those 
points where the phase is absolutely known or where the phase change is 
constant. By this means the whole phase function may be mapped. 

The measurement of the amplitude of the Fourier transform of the intensity 
distribution across the source may be made concurrently with the measurement 
of phase. We may form the function 


[AB] « LAC 
|BC| x |A|? 


where |AB|, |BC| and |AC| are the moduli of the fringe systems between 
channels A and B, B and C, and A and C respectively, and |A[? is the total 
power due to the source recorded in channel A. It has been shown (7) that 
the function 7» may also be written as 


Kf (ab) x f (ac) 
f(bc) 
where f(ab), f(dc) and f(ac) are the amplitudes of the transform for the baselines 


AB, BC and AC, and K is a constant of the equipment which may be checked 
regularly. 





3. Measurements on Cygnus (19N4A) 


(a) The East-West distribution.—A series of measurements was undertaken 
in position angle go in which the phase and amplitude of the Fourier components 
were determined at a number of spacings extending to the region of the second 
maximum of the transform. ‘Topographical features prevented the taking of 
readings in the immediate vicinity of the first minimum, though these readings 
would have been of considerable value in the analysis. At each spacing of the 
aerial systems the angular velocity of the rotating phase shifters was set to a 
suitable value for measurements to be made of the fringe pattern on the 
recording meters. ‘The relative values of the natural fringe frequency and the 
phase shifter frequency were noted, since the sign of this parameter determined 
the interpretation of the measured phase shifts as leading or lagging and hence 
the interpretation of the relative position of the asymmetrical component of 
the source. 
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The measured values of the amplitudes and phases of the complex visibility 
function corresponding to the Fourier transform of the brightness distribution 
are listed in Table I. 

TaBLe I 
Aerial spacing . 
peecnsat irs Amplitude Phase 
22 0°98 + 0°02 -4 
534 0°67 + 0°05 +2" 
828 0°51 +0°07 -—5° 

1060 0°32+0°05 +22 + 

1540 0°24 +0°03 + 166° + 

2160 0°46 +0°04 +190°+ 


Each of the readings tabulated was the mean of the observations obtained 
on several successive days, usually about ten in number. In addition each 
daily record yielded a number of separate measurements of the phase function 
from the individual fringe displacements, whilst the amplitude was also deduced 
independently from several points on the envelope of each fringe pattern. 

It will be seen that the amplitude of the second maximum of the transform 
of the source is comparable to, but slightly less than, the value previously 
obtained by Jennison and Das Gupta (6). The phase function undergoes a 


— 


° 
© 





© ° 
> 





v 
A . 4000 $000 


Fic. 1.—Cygnus A (19N5A), amplitude of visibility function at 127 Mc/s in position angle 90°. 
The vertical lines represent readings taken with the phase sensitive interferometer. The crosses 
are readings by Smith and the circles are readings by Fennison and Das Gupta; the readings 
between 4.000 and 5000 wave-lengths were zero in the latter measurements of the function p*, 
on conversion to p the errors are as shown. The two curves refer to the alternative brightness 
distributions in Fig. 3. 


sudden change in the region of the first minimum, in agreement with the 
assumptions made by the above workers in interpreting their results, but there 
is a detectable sinusoidal component contributing a small degree of asymmetry 
to the source. 

(6) Analysis of East-West distribution.—The readings given in Table I 
may be transformed to give the brightness distribution of the source smoothed 
by the convolution with the Fourier transform of the maximum aperture of 
2160 wave-lengths; only the major components of the source are therefore 
resolved but the solution within these limits, and within the errors, is unique. 
The readings obtained for the amplitudes are plotted in Fig. 1, together with 
the readings obtained by Jennison and Das Gupta and by Smith. The phases 
are plotted in Fig. 2. Superimposed on the amplitude and phase plots are 





178 R. C. Fennison and V. Latham Vol. 119 


the curves corresponding to the transforms of the distributions shown in Fig. 3. 
It will be observed that both these distributions transform within the error 
of the readings and it is not possible to distinguish between them without 
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Fic. 2.—Phase of the visibility function at 127 Mc/s, in position angle 90° across the radio 
source Cygnus A (19N4A). The vertical lines refer to the observed readings whilst the two curves 
correspond to the two distributions in Fig. 3. 


exceeding the baselines used. ‘The two distributions are similar in that the 
integrated brightness of corresponding components is the same, but in one 
distribution the diameters of the components are the same whilst the brightness 
of one of the components is increased, and in the other distribution the 
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FIG. 3.—Alternative primary brightness distributions at 127 Mc/s in position angle 90° 
across the radio source Cygnus A (19N4A). 


brightness of the two components are the same but the asymmetry is obtained 
by controlling their diameters. The separation between the centres of 
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radiation, 82 seconds of arc, is accurate to a few seconds of arc. The ratio of 
the intensities of the two components has an associated error of the order 
of 10 per cent. The major component, which we shall designate Cygnus Az, 
leads in right ascension. Though the distributions in Fig. 3 indicate only the 
primary regions of emission and cannot show the finer details, the contribution 
at 127 Mc's from the central reading of the system is evidently very small and 
cannot exceed a few per cent of the brightness of the « and 8 components. 
The errors in the readings combined with the smoothing imposed by the 
restriction of aperture do not enable the exact extent of this radiation, if it 
exists, to be determined from the present measurements. 

(c) The determination of the position angle of the major axis—The measure- 
ments of Hanbury Brown, Jennison and Das Gupta showed that the major axis 
of the Cygnus source had a position angle between go° and 120°. Measurements 
of this position angle by varying the azimuth of the baseline of a high resolution 
interferometer are liable to the large errors usually associated with changes of 
baseline in a two element system. In order to establish the position angle with 
greater accuracy three stations were used in the following ways. 

The intensity distribution measured by an interferometer system is that 
projected on to the tangent of a great circle on the celestial sphere which has 
the baseline as its diameter and which also passes through the source. As the 
source moves across the sky, the observed distribution will change. If the 
source is not radially symmetrical the observed fringe amplitude will be a 
maximum when the major axis is orthogonal to the plane containing the baseline 
and the source. The method used is therefore to observe the time at which the 
fringe amplitude reaches a maximum; the position angle of the major axis 
may then be found from the geometry of the situation. 

The main aerial system was phased to look in the approximate direction 
at which the maximum amplitude was expected and a second aerial system 
consisting of a simple dipole and reflector was sited at a distance of 4 km and 
an azimuth of 177°. A third aerial system, also consisting of a dipole and 
reflector, was sited 50 m from the first station at an azimuth of approximately 
315°. The phases of the signals from both outstations were artificially rotated 
to give suitable fringes on recorders multiplying their outputs with that of the 
home station. The small aerial sited close to the home station served as a 
monitor to remove the effect of the polar diagram of the large aerial system. 
The envelope of the fringes formed by mixing this channel with the signals 
from the home station depended only upon the polar diagram of the large 
aerial, whilst the envelope of the fringes formed by the large aerial and the 
other outstation depended primarily on both the polar diagram of the large 
aerial and the brightness distribution of the source. Suitable delays were 
included in the equipment to compensate for the distance of the furthermost 
outstation. 

The results of the observations are plotted in Fig. 4, in which the ordinate 
scale is the observed fringe amplitude after correction for aerial gain. Each 
point represents the mean of five observations. The measured amplitude 
corresponds to a position angle of the major axis of 97°, though this is subject 
to a slight systematic error arising from uncertainty in the correct setting of 
the delay lines. The diameter of the source in position angle 187° will be less 
than that observed by Jennison and Das Gupta in position angle 177°. Further 
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measurements at baselines in excess of 3000 wave-lengths in position 187° 
would be of considerable value. 
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Fic. 4.—Cygnus A (19N4A), 127 Me/s. Variation of fringe visibility with position of source 
/ f J 


the determination of the position angle of the major axis. Azimuth of the baseline of the 
rferometer is 177 The scale in the centre of the diagram gives the position angle corresponding 


the maximum fringe visibility 


4. Measurement on Cassiopeia (23.N5A).—The measurements made with 
the equipment were confined to the distribution across the source along an 
axis in the East-West direction. ‘The measurements were made concurrently 
with measurements on the Cygnus A source and with the time and manpower 
available it was not possible to take a fully representative set of readings. 
Considerable care was taken to establish the position of the first minimum, at 
which the fringe visibility was reduced below the noise level of the equipment. 
‘This minimum was found to occur at 840 wave-lengths; topographical difficulties 
prevented the measurement of visibility in the region of goo wave-lengths 
though a reading at this aerial separation was desirable to establish the position 
of the minimum with greater accuracy. 

‘The measurement of the phase function clearly showed the existence of the 
second and third maxima of the transform of the brightness distribution of the 
source. ‘These maxima could not have been detected from the amplitude 
measurements alone unless many more readings had been taken. Both the 


amplitude and phase functions showed a depression at short aerial spacings. 


This depression indicated an asymmetrical extension of the object seen as a 
projection into position angle go’. 

‘The values obtained for the amplitude and phase of the complex visibility 
function corresponding to the Fourier transform of the brightness distribution 
of the Cassiopeia source are listed in Table II. The magnitude of the errors 
in the phase function at large aerial spacings is due to the cumulative effect 
of the deviations from zero at short spacings (7). The readings for both 
amplitude and phase at 224 wave-lengths were the mean of a number of 
observations on slightly different baselines. In general each of the readings 
tabulated was the mean of the observations obtained on several successive 
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days, usually about ten in number. In addition, each daily record yielded a 
number of separate measurements of the phase function from the individual 
fringe displacements, whilst the amplitude was also deduced independently 
from several points on the envelope of the fringe pattern. 


TaB_e II 


Aerial spacing ‘ 
(wave-lengths) Amplitude Phase 
0°78 +0°05 ° 
0°36 +0°05 
0°045 +0°03 47° 
oo )6t 0°01 ashe 
0028 +0°01 280° + 32° 
0°05 +oO°O! 135°+43° 
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The observed amplitude of the visibility function is plotted in Fig. 5, in 
which the vertical lines represent the errors associated with each point. The 
broken curves drawn on this plot correspond to simple rectangular slit and 
disk models. They do not pass through the two points at large aerial spacings, 
but these are satisfied by slightly tapering the slit distribution or limb brightening 
the disk. 
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Fic. 5.—The amplitude of the visibility function of the Cassiopeia A (23N5A) source along 
the East-West axis. The dotted curve on the diagram represents the transformation of a simple 
disk source; the dashed curve represents a simple slit whilst the full curve is the transformation 
of the mean of the two distributions in Fig. 8. 


The phase of the visibility function is plotted in Fig. 6. From this function 
it is clearly seen that the points at 1540 and 2150 wave-lengths lie in the 
second and third maxima of the visibility function. The shaded portion in 
Fig. 6 encloses the range of brightness distributions lying between the two 
illustrated in Fig. 8. 

Fig. 7 shows the residual to the phase function after subtraction of the 
symmetrical component. The asymmetry apparent in the neighbourhood of 
500 wave-lengths will be seen to coincide with a depression in the amplitude 
function (Fig. 5) in this region. The general slant of the residual phase function 
is due to the shift of the centre of radiation to the main body of the nebula as 
the spur in position angle go” is resolved with extension of the baseline. 

The brightness distribution within the radio source in position angle go° 
is shown in Fig. 8. The two diagrams in this figure embrace the limits of 
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brightness distribution permitted by the experimental results. The most 
noticeable feature of the distributions is the faint spur on the left of the diagrams. 
This is seen as a projection into position angle go” and the true position angle 
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Fic. 6.—The phase of the visibility function of the Cassiopeia A (23N5A) source along the 
East-West axis. The shaded area embraces the region contained between the transformations of 
the tevo models shown in Fig. 8. The dotted line represents the phase of a symmetrical source. 
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Fic. 7.—The residual to the phase function after subtraction of the symme trical component. 


cannot be determined without further measurements taken with a baseline of 
different orientation. A few unpublished measurements, made on short 
baselines at Jodrell Bank in the course of testing the phase sensitive interfero- 
meter, failed to show a significant contribution in the North-South direction, 
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and it is probable that the spur actually lies close to the position angle go”. 
If this is not so, an allowance must be made for the projected length of the spur 
into this position angle; the true length will be somewhat greater than that 
indicated by the present measurements. 

The main body of the nebula appears to be reasonably symmetrical within 
the error of the present readings, but the readings imply that there may be 
slight limb brightening over the distribution otherwise given by that of a 
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Fic. 8.—Limiting brightness distribution for the Cassiopeia A (23N5A) source 
along the East-West axis at a frequency of 127 Mc/s. 


uniformly iiluminated disk. The presence of the asymmetrical spur does not 
appreciably affect the position quoted by Smith (3) for the right ascension of 
the object as measured at baselines of less than 200 wave-lengths. The phase 
function shown in Fig. 6 may be used to correct the position at larger apertures. 
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Summary 


The similarity between the spectra of cosmic rays and of the relativistic 
electrons responsible for cosmic radio waves suggests that the electrons may 
derive from the collisions of primary cosmic rays with the interstellar gas. 
The relation to be expected from this process is examined, and it is shown 
that the form of the electron energy spectrum is determined primarily by the 
different losses which prevail for different parts of the electron energy 
spectrum. 

It is found that the observed spectra of various different types of radio 
source agree with the theoretical spectra, and that it is possible to obtain 
information on various physical conditions in the radio sources. 





Introduction.—\It is generally agreed that non-thermal cosmic radio noise is 
emitted by high energy electrons accelerated in weak magnetic fields. This is 
the familiar synchrotron mechanism. From the observed variation of radio 
emission with frequency, in which the power in unit frequency interval is 
proportional to »-**, one may deduce directly that the number of electrons in 
unit energy interval, which we shall call 7(£), is proportional to E-“+®, The 
similarity of this spectral variation to that of cosmic rays, the number of which 
in unit energy interval is proportional to E~“+»”, suggests that the electron flux 
is related in some way to that of cosmic rays. Since electrons lose energy much 
faster than protons, no original electrons (that is, those injected with the protons) 
would be left in a steady system. They will be continually injected, though, as 
the decay products from the collisions of cosmic rays with the interstellar gas; 
it is the energy spectrum to be expected from this process which will be examined 
in this paper. 

We cannot expect the indices p and s to be identical, chiefly because the 
electrons’ rate of loss of energy varies very much withenergy. It is found that the 


effective spectral index s increases slowly with energy since the spectral index 


of the observed radio spectrum increases slowly with frequency. ‘The radio 
spectra of several sources are now well known and it appears that they accord 
well with the assumption that the electrons are produced by primary cosmic 
rays. Small deviations may be due to the effects of acceleration processes but, 
provided that these are unimportant, the radio spectra may be used to explore 
physical conditions in the various kinds of radio source. 





No. 2, 1959 Cosmic radio emission and cosmic rays 


Characteristics of synchrotron radiation.—Let 
a(E)dE= KE“+9 dE 
be the number of electrons of energy E to E+dE. The radio emission is then 
J, where 
K 4+/3é ae as , , : 
J,=—. \ 3¢ —3f_) Hite, y-e2 ergs sec~! (c/s)~}. 
2 mec \4nrmic® 
H is the magnetic field in gauss. J is an integral of the order of 2 which varies 
slightly with s. An electron emits most of its radiation near its critical frequency 
v,, Where 
v.= 6°3 x 10'8 HE? c/s. 


E is the electron’s energy in ergs. To account accurately for any given radio 
emission, therefore, we must know the field, because it affects both the energy 
of the electrons responsible for a particular frequency (through y,) and the 
number required (through J ,). 

Energy losses.—The determining factor in the general form of the electron 
energy spectrum is rate of loss of energy by various processes, since these 
processes show marked, and different, variations with energy. Consequently 
we shall first examine these losses and then their effect on an electron energy 
spectrum where the source of electrons at energy E is simply proportional to 
the proton number density, E~“*»”, 

Electrons may lose energy by (i) ionizing the surrounding medium; 
(ii) bremsstrahlung; and (iii) synchrotron radiation in weak magnetic fields. 

The first is predominant at very low energies, decreases quickly as the energy 
increases, and at energies greater than 10°eV is nearly constant. The very 
slow increase in the rate of loss is proportional to log E only. The energy is 
lost by small steps in many encounters with the surrounding medium, so that the 
rate of the energy loss may be treated mathematically as continuous. 
3remsstrahlung is more important at energies greater than the so-called critical 
energy (3°82 x 10°eV in hydrogen) but, unlike ionization losses, it occurs in 
large discrete events. 

In any accurate treatment one must use only the probability that an electron 
will lose a certain amount of energy in a particular length of path. The 
theoretical expression I shall use is valid in neutral hydrogen at energies much 
greater than the critical (3-8 x 1o°eV). The effect of using such an expression 
in the energy region just above 3-8 x 10° eV is to overestimate the loss by no 
more than a factor of two. Synchrotron radiation loss is, mathematically speaking, 
continuous, and is proportional to the square of the energy. In the usual weak 
interstellar magnetic fields, it is smaller than bremsstrahlung except at very high 
energies. ‘The conditions under which it supersedes bremsstrahlung will be 
discussed in detail in another section. 

Positive and negative electrons may annihilate each other; this process has 
a negligible cross-section at the energies under consideration*. 

The effect of energy losses on the electron energy spectrum. Qualitative 
discussion.—In Figs. 1 and 2, the various rates of loss of energy (in ergs per 
radiation length) have been plotted as functions of energy. Even though 
bremsstrahlung losses are, strictly speaking, not continuous (i.e. “‘ rate of loss ”’ 


* Heitler, The Quantum Theory of Radiation, O.U.P., 1944, p. 230. 
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Fic. 1.—Rate of energy loss (ergs per radiation length) plotted as a function of energy (ergs). 
(6) is the loss due to bremsstrahlung. 
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is a meaningless term), yet, on the average, an electron loses its own energy in 
one radiation length. This is a length of path which depends on the medium and 
its density. In neutral hydrogen, it is 69 gmcm~*. The curve R represents 
the average loss by bremsstrahlung. J is the continuous loss by ionizing 
collisions, H is continuous loss by synchrotron radiation for various values of 
H?/N, where N is the number density in cm~*, The total loss curve, therefore, 
is fairly flat in the energy region (a), linear in (b) and parabolic in (c). 


A. 








Log Frequency 


Fic. 3.—Hypothetical spectral indices: (i) acceleration never predominant; (ii) acceleration 
predominant for energies within a range varying by a factor of nearly 10 to 1. 


Suppose now that there is a source of electrons akin to the cosmic rays in 
spectrum, say of index p. Let the source be OQ,E~“+”, In a steady state, the 
number of electrons leaving a particular energy interval is the same as the number 
gained. Therefore 


a m(E) {rate of energy loss by one electron} + Q,E-"*+” =o. 


The rate of energy loss is the sum of a constant, a linear and a quadratic term, 
call ita+bE+ck*. Then 


+p -1 
(B)= 2 "(j++ eb) : 


Compare this with the electron energy cee is: which has been set 
proportional to 
E-A+9, 


14 
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The effective index s of the flux of electrons whose energy is in region (a), 
where the constant term of the loss curve predominates, is p—1; the index 
is p, the same as in the source for electrons whose energy is in region (6), 
where the linear term of the loss curve predominates; the index at the highest 
energies, region (c), is p+1. ‘Then the curve of log radio power against 
log frequency for a radio source due to such electrons would look like Fig. 3 (i). 

The effect of acceleration.—Electrons will be accelerated in exactly the same 
way as other charged cosmic-ray particles, but the effect may be negligible 
compared to the effect of the losses. ‘There are many theories on accelerating 
mechanisms (see, for example, Parker 1958, for a resumé), all of which predict 
the same linear dependence of rate of ‘energy gain on energy. It is sufficient 
for the present therefore to define a coefficient of acceleration « such that «£ is 
the energy gained per radiation length by an electron of energy E. So far as 
we know from the observed index of primary cosmic rays and assumed collision 
cross-sections, « does not exceed unity. This it must do if electrons are to be 
accelerated independently, that is, if the rate of acceleration «EF is to exceed the 
approximate rate of bremsstrahlung loss, FE. For a purely hypothetical case, 
where «= 2, the loss curve and the equivalent radio spectrum have been plotted 
(Fig. 2 and Fig. 3(iii)). The accuracy of present flux measurements is not 
sufficient to determine whether any radio spectra are curves of this sort or not. 

The production of electrons from cosmic rays.—In detail, the discussion above 
has not concerned itself with the precise nature of the source of electrons. It has 
merely assumed that they are secondary particles of the primary flux of heavy 
cosmic rays. We know very little of the cross-sections for meson production 
and especially their variation with the energy of the incoming particle and the 
outgoing electron. The following attempt to devise a source function is therefore 
at best only an approximation. 

We know that an electron of energy E is the result of the collision of a proton 
of some higher energy E’ with the nucleus of an atom of interstellar hydrogen. 
‘The available energy will be distributed among the residual nuclei, mesons and 
neutrinos. ‘The electrons are, incidentally, positive since no 7~ mesons are 
produced in the p+/ collisions in the interstellar medium. 

We shall assume first that EF is proportional to E’ at energies too low for 
multiple mesons and strange particles, i.e. 


E’ = aE. 
At very much higher energies, one may use Fermi’s statistical theory (1950). 
Electrons of energy E are produced with a multiplicity s from a collision of 


energy E’. 


2\ 4 
E=cE'*4, c¢=0'049 (=) 
2 


where 


2 


s=bE'™, b= 8( 


) “ae 


Secondly we assume that the number of collisions is proportional to 
number of cosmic-ray protons, i.e. 


P(E’) = P,E’-8+»), 
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Thirdly, that at the lower energies, where E’=aE, the proportion of 
meson-producing nuclear collisions increases with energy. As an approximation 


to the curve reported by Riddiford (1957), let us say that 


{2} 


is the probability that a proton of energy E’ produces a meson in one radiation 
length, where o is the total cross-section gmcm~*, K is (69gmcm~*)~', 7 is 
a constant of the order of }. This expression is zero, if E’ <E,=300 MeV, 
the threshold for meson production. 

If E’>10' eV the source function used will be 


; b /E\~281+2p) 
pe Sg 5. (1) 
. ¢°te 
There are energies between 10!2 and 10!eV about which little can be said. 
If E’<10!' eV the source function will be 


KoPya{t “ (=) “ (aE)-1+0), (2) 


0 


The equation governing the electron flux is therefore 


on(E x a . , 
y= 2- J ir) Aa(- =.) }W(v)do+ 4 m(E){ey—2E+ BE (3) 
where OQ is the source, either (1) or (2), t is the distance along the electron’s path 
measured in radiation lengths, %(v) is the probability that a fraction v of the 
electron’s energy will be lost by bremsstrahlung in one radiation length. By 
definition, its average value is unity. €, is the energy lost by ionization per radiation 
length. It will be assumed that it will be constant and of the order of the critical 
energy. aE represents the possible energy gained by acceleration per radiation 
length. BE? is the energy lost by synchrotron radiation per radiation length. 
If E is in ergs, the magnetic field H is in gauss, and the number density of the gas 
is N cm~’, then the coefficient B is 3 x 10! H?/N. 

Solutions of the equation.—To simplify the problem first of all, only stationary 
solutions have been looked for; that is, d7/dt has been set equal to zero. 
Secondly the equation has been solved in different energy regions where the 
appropriate approximations may be made. 


E> 10%eV. 

Synchrotron radiation predominates as the loss mechanism at these electron 
energies even in fairly weak fields. ‘The energy E’ of the parent nuclei is 
sufficiently large for Fermi’s statistical theory to be applicable (1) 


. «“ePe 6 /E\-et® I 
a(R) = Sed oe: (; 
P 9-3 Oe ee (4) 
BE 
1o'eV >E>s5 x 10% eV. 
The continuous losses by ionization and by synchrotron radiation have been 
neglected, and the lower energy electron source has been used (2). 
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n(E - (FE) poe. ee : 
= ENE) *a—Ae) (P=) P+ 7h) 
(aE)? aE\~" A(p)—ap 
+ ge —(1-(F) a —)} () 
A(p)—ap\ \Eo) A(p+7)—a(p+7)/J 


1 
A(s)= | (1 —v*)b(v) dv 
0 


and A’(s) is its derivative with respect to s. A(s) is tabulated by Rossi (1952), 
p. 296. 5, is the root of A(s)—as=o. 

If «, the acceleration coefficient, is small, sy is large and the first term of (5) 
varies too rapidly with energy to be important. ‘The electron spectrum is 
governed entirely by the second term, that is, by the proton flux. Physically 
this means that the electrons are kept up by their source in the decayed mesons 
from cosmic-ray collisions against energy losses predominantly by bremsstrahlung. 

If, on the other hand, « is of the order of unity, the root 5) is less than p. Now 
it is the first term of the solution (5) which is important. Physically, this means 
that acceleration is overcoming radiation losses and the electron flux can grow 
independently of its source. 
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where 


3x 108/a<E<5 x 10%, 


lonization losses are now predominant. ‘The spectrum is flatter and joins 
up with the second term of the previous solution when E =, 


» mor, : I ak\-" 1 
n(E) = staal (ak) P = 2S —_ a ° 

€9 Pp E, ptr 
If « is large, ionization losses do not overtake acceleration until lower energies 
are reached, 1.e. aE = €. 


E< 3x 105/a eV. 


The source is negligible at these energies and 7(E) x (ionization loss) is 
constant. 

Non-relativistic energies are of little interest so far as radio emission is 
concerned, for the intensity drops off very quickly with decreasing energy. In 
addition, the simple formulae relating spectral index to the index of the electrons 
no longer hold. 

Applications 

1. The Galaxy.—We may now examine the relation between the spectral 
indices of cosmic rays and galactic radio waves using the following information. 

(i) At the Earth cosmic rays are observed having an energy density of the 
order of 1 eV/cm®* and an index p which varies slightly with energy. 


E Pp 
10° e\ 1:2 (Popov 1957) 
101! 1°5 
10! 20 


(ii) The background radio emission from the Galaxy is due to an electron 
flux of index s=1-0+0°2 (Adgie and Smith 1956). 

(iii) The emissivities at 80 Mc/s of the halo and disk are 1-8 x 10° and 
3°2 x 10° w.s.~! (c/s)! respectively (Baldwin 1955). 
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(iv) The estimated magnetic field in the spiral arm is 6 x 10~* gauss (Fermi 
and Chandrasekhar 1953). 

(v) On the average the number density N of the interstellar gas is roughly 
unity; hence H?/N is roughly 2 x 10-". 

From (ii) and (iv) we may calculate that it is the electrons of energy 10° eV that 
are chiefly responsible for 80 Mc/s emission, and that their energy density must be 
roughly 0-035 eV/cm? in the disk and arms. The same energy density can account 
for the halo, if the field there is 10~-* gauss. Provided of course that the Galaxy 
is in a steady state, the electrons of energy 10° eV in the galactic field should obey 
equation (5), since their energy losses are mostly by bremsstrahlung. The field 
is not large enough for the synchrotron losses to be predominant. 

If we neglect acceleration and hence the first term of equation (5) we have 


7(E) akoP,(aE) +») { A(p)—ap (=) 
 Alp)=—ap A(p+7)—a(p+r)\ Ey } 


as the solution for the electron flux. ‘The second term in the brackets is small 
and the bracket is ~1. 

Now the constant a was introduced quite arbitrarily to represent the efficiency 
of energy transfer from the colliding proton (E’) to its daughter electron of 
energy E. Because the electron is the result of three distinct decay processes 
its share of the original energy E’ may be quite small, 10 per cent or less, 








pt+p-pt+n+nt 
>p+p+n7 
ntopt+y 
pt—>Br+v. 
If then E’ =aE, ais large; the index p is that appropriate to protons of energy 
ten to a hundred times the electrons’ (F), i.e. 


V2<p<i'5. 


The observed index is 1-0 + 0°2, which is not far from the expected value. The 
difference, if significant, may mean that acceleration should not have been 
ignored. 

2. Cassiopeia A.—The radio emission from Cassiopeia A has been shown to 
come from a galactic nebulosity identified by Minkowski (1958) as a supernova 
256 years old. It is therefore not possible to assume that a stationary solution 
for 7(E) applies. Nevertheless, the slope of its radio spectrum, —0°8+0°05, 
suggests that the index s of the electrons is 1-6+0-1, a value not significantly 
different from the value of p for galactic cosmic rays. 

3. Cygnus A.—The radiation is from two colliding galaxies whose radio 
spectrum exhibits a sharp break at about 350 Mc/s. Above this frequency, 
the spectrum has a slope of —1-02+0-1 (corresponding to s=2:04+0-2) and 
below, the slope is —0-66+0°6 (corresponding to s=1°32+0°12). Any 
observed steepening of the radio spectrum at the high frequencies must 
indicate a real change in the index s of the electrons for the only confusion 
which might arise is thermal emission which of course varies in the opposite 
way with frequency. On the other hand the flattening at the low frequency 
and around 10 Mc/s may be due to absorption in H 11 regions through which the 
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radiation has passed, and not due to a real flattening of the electrons’ spectrum, 
which is expected as ionization losses take over at the low energies corresponding 
to 10 Mc/s. 

The break in the spectrum at the high frequency end is interesting. It 
should occur roughly at the critical frequency of those electrons whose synchrotron 
losses are just overtaking the bremsstrahlung losses. If their energy is FE, the 
former is 3 x 10'?(H?/N)E®? ergs per radiation length, and the latter is E in the 
same units. ‘Their critical frequency v, is given by 


v= 6:29 x 10 °HE? c/s, (6) 


E = 3 » ron ee ergs. 


If E is eliminated between (6) and (7) we obtain 
H =8-9 x 10-3 N23y,-1, 

\ is the number density of the hydrogen gas in the source and may be of 
the order of unity. Putting v.=350 Mc/s, we find H = 1-3 x 10~° gauss. 

It is possible to derive an estimate for the magnetic field quite independently 
by using the following data. 

(i) s=2. 

(ii) ‘The power at 80 Mc/s is P,) = 1-87 x 107 watts c/s~! ster~". 

(iii) Distance 93-5 Mpc. 

(iv) Volume 3 x 10% cm?, 

If it is assumed that the electron density and magnetic fields are the most 
important components of the total stored energy in the galaxies and that this 
total energy has its minimum value consistent with the power emitted, then 
the magnetic field is H = 5 x 10~° gauss. ‘The two estimates are reasonably close. 

4. Other sources.—The spectral indices of ‘‘point’’ radio sources taken as 
a whole vary between rather narrow limits: —o-5 to — 2-0. The corresponding 
electron index s lies between 1 and 4, and therefore limits on p may be set 
immediately: p< 5. On closer inspection, however, it appears probable that 
the real limits are much stricter. I have shown that all radio spectra may be 
expected to be steeper at high frequencies. Whitfield has observed that the 
mean index of galactic sources is less than that of the identified extra-galactic 
sources and that both are less than the mean index of unidentified sources*. 
Following this observation, Ryle has suggested (in the Bakerian Lecture to the 
Royal Society, 1958) that most radio sources are extremely distant and that the 
steeper spectra expected at high frequencies are brought into view by the red 
shift. His suggestion is supported by Dewhirst’s identification of a number 
of sources whose distances are so great, and consequently whose radio powers 
are so enormous, that the origin of their energy is difficult to understand if the 
better part is wasted in bremsstrahlung and ionization losses. Synchrotron 
losses therefore, must predominate, and hence it is the steep part of the spectrum 
where s=p+1 which is observed. If so, stricter limits may be set on , i.e. p < 3. 


* The spectral index —o-28 of Taurus A or the Crab Nebula is exceptionally low. The 
corresponding electron flux index is s=o-56. But in this source, a remnant of a supernova 
explosion 900 years ago, no steady state obtains. A radiation length is equivalent to 5 x 10°/N 
years. ‘The number density cannot be greater than 107; therefore, one radiation length is 


5 x 10* years, a value very much greater than the age of the nebula. 
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It is interesting that the index of the local flux of cosmic rays also lies 
within these limits, but it is surprising that they may be generalized to such 
widely differing regions of space such as most radio sources. This conclusion 
is based of course on the assumption of a stationary flux of cosmic rays which may 
be unjustifiable. 
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PESIOME JKOKJIAJIOB, B TIEPEBOJIE HA PYCCKHH A3SbIK 


HYTAIHA: CPABHEHHE TEOPHH C HABJIIOIJEHHAMU 
X. JIonehpuc 


JlaioTtcA TeopeTH¥eCKHe 3HaueHHs JIA OCHOBHOTO, NONyrOMOBOrO HW MOJIYMeCAYHOTO UeHOB 
HYTalMH, OCHOBAHHbIe Ha MOJIeIAX Gemsin: MOC Pomla nH TOUcUHOH. OHH CpaBHHBaJINCcb co 
3HaYCHHAMH, MOUYYCHHEIMH “PeyepoRLimM H3 HaOOAeHHH. Cornacne nonyunsoch yume, 4em 
BO BCeX TIpebIyuIMx Cyuanx. LlonymMecH4HbIe WieHbI HMeIOT MOUTH OJ(MHaKOBbIe 3Ha4eHHA 
win O6eHxX MOjesIeH WH CormacyloTcH c HaOsnojyeHHem. Besmunnbt O6OHX KOMNOHEHTOB Mosty- 
rOMOBOH COCTaBIIALOMIeH JIyulle BCerO COrslacyloTCA CO CxemMOil TOYeUHOH MOyemH. HabsmoyeH- 
Hoe 19-JieTHee H3SMCHCHHe B HaKJIOHCHHH JIy4ule BCero CormacyeTcA C TOUeCUHOM MOebIO, a 
H3McHeHHe B JOIroTe—c MOjlebIO Poula; HeKOTOpOe pasHOrsiacHe C HaGJHOJE€HHAMM JOJDKHO 
OcTaBaTLCcA H JIA TOGO Apyrok MojecIH. 


PEKOMBHHAIHWA HUSITIYYMEHHEM BOJIOPOJIHbIX HOHOB 
M. Cumon 


YuHTbIBad NepBble TPH WieHa AaCHMMNTOTHYECKOrO pasiOx*KeHHA MHO*KHTeIA Kpamepca-l'ayuta 
MpoOBeeHO BbLIUHCeHHe CKOPOCTH PeKOMOHMHAaLHH HW CpewHelt KHHETHYeCKONM 9HEPprHH peKOMOHH- 
HPYIOULHXCA BICKTPOHOB. 


PEINEHHE YPABHEHHM KACKAJIHOrO 3AXBATA JUIA BOTOPOTA 
M. Cumon 


BbruncienvA mpoBefensi Wi T= 2,5 x 10°, § x 10°, 1 104 nH 210" °K, yuntTbipam Geckon- 
euHoe YHCJIO YPOBHeH, OJHakKO NpHHHMaA bni—[(2/4 1)/n*]bn. TlokasaHo, u4To pe3ynbtTaTel, 
nosyueHHpie belikepom n Mensesom, ourmGouHbr npHOnHsuTesbHO Ha 5%, Wia Cnyuaa A w Ha 
20%, qua Cryuasn B. 

YolyuuleHHbili MeTO, BLIUHCICHHH NOuy4uaeTCH, CCIM BBECTH KaCKajIHylO MaTPHIly, 3JIEMeHTbI 
KOTOpOH Cn’n QalOT NOMHY!O BEPOATHOCTL Mepexosa n’--n C YUeTOM BCEX BOSMOMKHbIX KaCKaHbIX 
nyTeH. 


HABJIIOMEHHA COJIHEVHOHM TPAHYJIAWMH C AOPOCTATA C SKHTIAKEM 
Jl. Baexyaaa, Jl. JTowxepem u A. Jloacgrwc 


Ilo-BHAMMOMY CYIIeCTByeT NpakTHYeCKHH Mpesel paspewieHMio, AOCTHXKHMOMY TIpH dpotor- 
padupopannn Cosma c Gemsin. TI brranch NOYYNTh 3HAaUHTeJIbHO JIyUlee paspellleHHe, AaBTOPbI 
MIpOHSBOAHIM (poTorpadHpoBaHHe TIpH MOMOUIH 29-CM pedpakTopa, MOHTHPOBaHHOrO Nos 
KOP3HHOH aspoctata, Ha BbICOoTe ~ 5500 M. Bpisio ocyulecTBJIeHO {Ba NoONeTa; BO BPpeMA BTOporo 
nmoueta, I ampesia 1957 r., Gbim0 CHATO 480 KajIpoB. XOTA yCNOBHA BHMMOCTH Ha Gemsie BO 
BpeMA NoseTa ObLIH OUCHb MWIOXHMH, (poTOrpadun, NONy4eHHble C aIpOCTaTHbIM TeJIECKONOM BO 
BPpeMA BTOPOrO MOeTa, MPeBOCXOAT MO KayecTBYy JIyYYIHe CHHMKH, H3 BCeX NMNOJYYeCHHbIX Ha 
3CMHBIX CTaHIMAX. 

CyomwHand CTpyKTypa dotocdeppi wunoctpupyetca cbotorpadusamu. Jlanbi apToKoppe- 
AWUMOHHBI€ KPHBbIe JUIA JByX (POTOMeTPHYeCKHX pa3pe30B HaHJIyUWIHX CHHMKOB. ‘TouHBIe 
(PYHKUHH MponyckaHHA {UIA TeJlecKoNa B esIOM (OOBeKTHB, OKYJIAP, (boTorpadnyeckan IMYyJIbCHA) 
ObLIM H3MepeHb! Mp MOMOUIM CHHYCOM{@JIbHbIX pellleTOK, a NMOJYYCHHbIe pesyJIbTaTh! O_VIK 
HCHOJIbSOBaHbI {VIA HMCIpaBlieHHA MHKPOdOoTOMeTpH4eCKHX H3MepeHHH KOHTpacta. CpeyHee 
HCIIpaBJIeHHOe SHAYCHHe KOHTPacta MO>K AY PpaHy1aMH MH Me?KIPaHYJIBHbIMH OOJIACTAMH COCTaBIAeT 
40% Hdl A 5300. 
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3BE3],bI TIIABHOM TMOCIIENOBATEJIBHOTCH 
C. Xetizeaepoe u DM. Xoiia 


MnterpupoBanne MpopeeHo WIA OAHOPOAHBIX 3Be3g TuNOB Iu II. Ilo-pugumomy umpnna 
NOOChI, pasqeAOUleH THMbI 3Be3{HOTO HaCeNeHHA, COOTBETCTBYIOIIaA HYJICBOMy BOSpacTy B 
mpeyesaxX riaBHolt MOCIeOBaTebHOCTH, He MpeBbNaeT 0,1 3BE3THOM BENHUMHBI, CCIM TOKO 
HCNOUIb3yeTCA MpaBHJIbHad WikKasla 9dpcdpeKTHBHBIX TemMepartyp. 


BOSPACTbI CYBrHrAHTOB THITIOB I H II 
D. Xoiia 


OBOJTOUMOHHBIe MOC OBaTeIbHOCTH, OUeHb OM3KHe K HAYKHHM YaCTAM 

a) mlocnejOBaTebHOCTel WIapoBbIX CKOMIeHHH, 

é) mocnegqoBatenbHocTH M 67, ObLim NocTpoeHbI Np NMOMOUIM BLIUHCIMTeIbHOM MalIMHbI 
IBM 704. OTH HwKHHE “aCTH NOCcIesOBaTeNbHOCTeM ONpeeNAIOT BOSpacTHI 3Be3y THNOB I u IT. 

Hicnonb3yA MOSyYeHHbIe H3 HaGOAeHHi SBe3HbIeC BCIMUMHbI WIA OM3uexKalMx 3Be3q 
Tuna I c OoubuIMMH TPHrOHOMeTPHUeCKHMH MapaslakcaMH MW Hallarad yCOBHe, YTO 3Be3{bI THIMAa 
II qospKubI GbITb crapute 3Be39 THNa I, aBTOP MPHXOHT K CJIe€AYIOUIMM BbIBOWaM: 

1. Ilepemenupie Tuna ckonsIeHHii He MOryT GbITA Apue My = +0, 6. 

2. Bos3pact TanakTuku 2ospKeH Opts Gombe 107° mer. 

IIpupeyeno takoKe coo6paxkeHue B NOMb3y TOrO, UTO NepBOHAYAIbHOe COepoKaHHe TresIHA B 
3Be3 jax THA II 6110 HHU3SKHM. 


S3AMETKA O ITTEPEMEHHbBIX THITA RR JIMPBbI 
T. Kunman 


Ilynpcaunonnpi Kputepuit P\/p=Q ucnonb30BaH JIA NomyaeHnaA abcOMIOTHBIX BeJIHYHH 
MepeMeHHbIX THMa iXKi.vr Kak (pyHKUKH HX WBeTa MH NepHosla; HcpaBlieHa ommMOKa, JONyueHHaA 
Can HJpKeM B aHasiorH4uHoM BbIBOZe. CormacHo 3TOi MOyemM HaOmO,eHHOe pasiMune B CpeqHHX 
llepHojlax NepemMeHHbIX B M 3 4 w Cen oObACHAeTCA pasIMuMeM B O™,3 abcouIOTHOH BH3yasIbHOH 
BCJIMUMHbI M@*KLy 3Be34aMH 9THX WZByXx cKomeHHH. Tlomyueno rpy6o npuOsmpKeHHOe Cormacne 
M@KIY MPeBbIUHCICHHbIMH aOCOJIOTHbIMH BeJIHUMHaMH HW BeIHYMHaMH, NMOJYYeCHHbIMH 10 
TPHTOHOMETPHYECKHM TlapasLiaKCcaM JIA UeTbIPeX MepeMeHHBIX 3Be3qHOrO NomA Cc P<O4,2, 

IlybcallHOHHadA MOjleIb MpecKasbiBaeT passMune B O",7 MexKiy uedbengaamH Tuma ab c 
@aHOM@JIbHO KOPOTKHM NepHojsioM, OOHapyKeHHBIMH ['anouIKHHbIM B HallpaBJieCHHH raslakTHYecKoro 
ICHTpa, H HOPM@JIbHbIMM eden amu THNOB ab uc B TO xe OGnacTH HeGa. Tak Kak 9TOrO He 
HaO m0 1aeTCA, TO WIM IyJIbCal|MOHHaA MOJ{eJIb, HJIM BeCJIMUMHbI NepHOOB AHOMAJIBHbIX MepeMeHHbIXx 
HeBePpHBbI. 

B cnvicke NpHBOAHTCA 16 3Be3, THNa ab c nepHogzom P<04,36 B oGuem none. TloKka3ano, 
UTO 9TH 3BE3jbI eKaT OsMoKe K rallakTH4eCKO MIOCcKOCTH, 4eM JIerxKaa ObI Cry4aitHaa BBIGOpO- 
4uHaA Ppylila HOPMaJIbHbIX MepeMeHHbIX THMa ¢c TaKOro >Ke NepHosa, B CormacHH Cc OOHapy>KeHHbIM 
Ilerpe OTHOCHTeJIbHO MCHbUIMM COJIHC4HbIM JBHOKeHHeM JIA 9THX 3Be30. HeoOxoqMMbI [aJib- 
HelilmMe HaOOJeHHA JIA NOAOOHbIX 3Be30, Pacnos107KeHHbIX Kak B HalipaBJIeHHH raslaKTHYecKoro 
IlC¢HTpa, Tak H B o6ulem mMovte. 


CHUCTEMA ¢ ®EHHKCA 
I’. Xetizman 


Jly4esble CKOpOCTH H3MepeHbI 110 71 HOBO cneKTporpamme ¢ Phe, nomy4eHHod Ha o6cepBaTopuHn 
Maynt Crpomio. Bprancnena cnekTpockonmueckas opOuta riaBHol 3Be3qbI M M3 AMd>depenun- 
Q@IbHbIX H3MepeHHii ompeseneHa nouyammiutyfa K, cnyTHuKa. AGcomoTHBIe pasmepbI OGoMx 
KOMIIOHCHTOB YCT@HOBJICHbI H3 CONOCTaBJICHHA CIICKTPOCKOMMY4CCKHX op6uT H cboToMeTpHuecKHx 
HaOQmoyeHHit Xorra. OnpefeneH cnekTpanbHbIi Kaacc cnyTHHKa. Hekotoppie TpyqHOCTH 
TeOPHH CHCTCMbI M10-BHAMMOMY Y/laJ1OCb PaspeLllinTb. 
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3AMEUAHHA O HATH IOOKHBIX TASOBbIX TYMAHHOCTAX 
JI. C. Dean 


Coo6miaiorca pesvibTaTb! HaOmMOoWeHHIt CNeKTPOB, JIVUCBLIX CKOPOCTeil, BHYTPCHHHX JIBHK- 
eHHii H COpMbI [VIA NATH 1OOKHBIX OObEKTOB MGC 3918, 6302, 6326, IC 4776 4 Anon «= 16" 13",5 
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IWAPOBbIE CKOILIEHHA. I. JIVUWEBbIE CKOPOCTH IOXKHBbIX IMAPOBbIX 
CKOINIEHHA 


1. Kunman 


HbIX IaPOBbIX CKOILICHHHM HW OOCyxKWalOTCA MX OWIHOKH. 

eCHOBCKOrO CHeKTpockona (jucnepcHA 86 A/MM) nosly4eHbI 

M YHCIe WIA 1§ CKOMICHHH B HHTerpasIbHOM cBeTe, CIIeK- 

CKOMICHHIT H CHeKTporpaMMbI OOOMX THMOB JWIA Tpex 

(OILICHHA paHee BoBce He CbLIO M3BeCTHO cKOpocTteH. IIpu 

) crieK Kora ¢ ucHepcneit 380 A’MM ObILIM NoOlyUeHbI CleKTporpaMMbl 
UIA JIBYX H3 KOTOPbIX MMCJIHCb Cle¢KTPOorpaMMbl, MostyueHHbie c Kaccer- 
<OTIOM (PaBHCHHe HaUIMX JTY4eBLIX CKOpocTeii c AaHHbIMH Meifanua u 
eT 1acue. IIpusogurca taGnMua Jy4ueBbIX CKopocteit u 
bIX CKOILICHHH. 


B THPEQEIAX HCTOUHHKOB 
I9N4A) Hl KACCHOIMEA A (23NsA 


iu B. Jlamam 


Hu (pasbl MpeoOpasopanna ype pacipesenenna 
HMHTCHCHBHBIX HCTOUHHKOB payHoustyuenna JleGenp 
OCH B HallpaBJieHHH BocTOK-3anay. Pacnpeyzen- 
{aeMbIM OTCUcTaM. IloKa3aHO, 4TO B CJIy4ae HCTOU- 
‘CTBOBaHHeM ByX OONbIIHX M3Ty4uaIoUMXx oGsacTeil 
OMMOHEHT, Oepe7KaloulMii MO MpAMOMY BOCKXOsK]ICHHIO, 
LIBAIOTCH TaK7Ke H3MCPCHHA, CBASaHHbIe C OMpeeseHHem 
OueBH HO {WMamMeTp HcTOUHHKa Jle6eb B MOSHIMOHHOM 
OTHOCHTCJIBHO OUCH HeoobWHM. H3amepenna jwiA 

ro cuaGoro BBICTyNa BONM3H NOSHIMOHHOrO yriia 90°, 

) YBCAINYCHHA APKOCTH K KpatO OCHOBHOrO HCTOUHHKA. 


CBASb KOCMHUYECKOPO PATMOHSJIYYEHHA C OIIEKTPOHHbIM 
KOMITOHEHTOM KOCMHUECKHX JIYUEH 


laumep 


ICCKMX JTV4CH H PCJIATMBHCTCKHX IJICKTPOHOB, OOYCJIOB- 
‘; HABOAHT Ha MbICJIb, YTO 9JIEKTPOHbI MOryT MOABJIATBCA 
© KOMIIOHCHTAa KOCMHYCCKHX JIY4YeCH C Me@?K3BeC3HbIM 

MO7KHO OYKHaTb TPH 9TOM TIporecce, MH MOKa3aHoO, U4TO 

‘TPOHOB OlpeseaeTca mpexnjte BCero pa3/IM4YHbIMH 
bIX UACTAX BHEPreTHUCCKOLO CIeCKTpa IICKTPOHOB. 

Pbl HCTOUHHKOB pajIHOH31y4YeHHA pa3JIMUHbIX THMOB 
bIMH CII€KTPpaMH HM UTO HMe€eTCAH BO3MO?*KHOCTL NOJIY4aTb 
OBHAX B 9THX HCTOUHHKAaX payMOn3sly4ueHHA. 
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